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The  National  Bureau  of  Standards  haa  been  engaged  in  the  heat 
transfer  analysis  of  underground  installations  for  the  past  several 
years.  This  report  covers  a  part  of  the  Bureau's  activities  related  to 
the  computer  simulation  of  the  thermal  environment  for  prototype  shelters. 

The  computer  was  used  basically  to  simulate  energy  balance  in 
the  shelter  living  space  and  to  analyze  heat  conduction  from  the  shelter 
walls  (including  ceiling  and  floor)  to  the  surrounding  earth. 

For  the  heat  conduction  analysis,  finite  difference  techniques 
were  employed;  using  a  three  dimensional  model  in  some  cases  and  a  one 
dimensional  model  for  the  remainder. 

Digital  computer  programs  were  developed  and  applied  to  seven  different 
prototype  shelters  for  which  temperature  and  humidity  records  with  simulated 
occupants  were  available  as  a  result  of  studies  by  the  National  Bureau  of 
Standards  and  by  the  University  of  Florida.  In  the  seven  shelters  used 
for  the  investigation,  twelve  different  operating  conditions  were  analyzed. 

Of  these  twelve  conditions,  fen  were  under  sumaer  operation  and  two  under 
moderate  winter  conditions. 

Generally  the  agreement  between  the  computed  and  observed  thermal 
environment  on  these  prototype  shelters  was  surprisingly  good,  in  spite  of 
the  fact  that  numerous  simplifications  were  Involved  in  describing  the  complex 
shelter  heat  transfer  system  for  computer  sna lysis.  TWO  Inherent  uncertainties 
exist,  vhlch  Influence  the  final  reliability  of  the  calculations.  The  first 
involves  the  description  of  the  actual  complex  system  by  mathematical  language 
(or  operat'  nl  uncertainty).  The  second  is  related  to  the  accuracy  of  input 
data  used  for  the  calculations,  (or  data  uncertainty).  Often,  these  are 
interrelated. 
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CUT  ALONG  THIS  LINE 


In  this  study,  the  following  four  different  computer  models  were 
studied  and  two  of  them  were  extensively  utilized  for  the  comparison 
of  calculated  thermal  environment  with  the  observed  data  in  the  proto¬ 
type  shelters. 

M-(l)  Three-dimensional  rectangular  model  with  composite  walls 
and  with  separate  initial  temperature  patterns  normal  to 
the  six  bounding  surfaces. 

M-(2)  Three-dimensional  rectangular  model  with  homogeneous  heat 
conduction  medium  having  separate  initial  temperature 
patterns  normal  to  the  six  boundary  surfaces. 

M-Q)  One-dimensional  compound  model  for  six  composite  wall 
systems. 

M-  (4)  One-dimensional  compound  heat  conduction  model,  same  as 
M-(3)  except  that  the  roof  region  was  assumed  adiabatic. 

To  simulate  the  initial  earth  temperature  distribution,  the  follow¬ 
ing  three  modes  were  employed: 

I-(l)  Earth  temperature  gradients  normal  to  the  six  bounding 
surfaces. 

I-(2)  Earth  temperature  gradient  normal  only  to  the  ground 
surface. 

1-0)  Initial  e  rth  temperature  constant  around  the  shelter. 

One  of  the  factors  not  well  established  for  calculating  the  shelter 
heat  transfer  is  the  heat  exchange  between  the  shelter  air  and  the  inner 
surfaces,  between  the  occupants  and  the  surfaces,  and  among  the  surfaces. 


The  analysis  of  simultaneous  exchange  for  radiative  and  convective 
energy  among  occupants,  air  and  Inner  surfaces  of  a  shelter  Is 
very  complex,  and  it  requires  the  solution  of  a  set  of  integral 
equations  which  are  difficult  to  solve  for  even  very  simple  geo¬ 
metries.  Therefore  conventional  combined  heat  transfer  coefficients 
for  radiation  and  convection  were  used  in  the  analysis.  Several 
numerical  values  and  combination  of  these  combined  coefficients 
were  assigned  to  the  six  interior  surfaces  to  study  the  overall 
effect  on  shelter  thermal  environment. 

Findings 

1)  A  soil  analysis  of  the  earth  around  most  of  the  prototype  shelters 
indicated  that  the  thermal  dlffusivlty  and  thermal  conductivity 
were  in  the  neighborhood  of  0.02  ft^/hr  and  0.75  Btv/hr,  (ft)**, 

°F/ ft  respectively.  These  values,  in  turn,  seem  to  result  in 

a  good  agreement  between  the  calculated  and  the  observed  earth 
temperature  change  surrounding  prototype  shelters. 

2)  The  following  combined  heat  transfer  coefficients  at  the  shelter 
inner  surfaces  produced  satisfactory  simulation  of  the  shelter 
summer  environment  for  most  of  the  prototype  shelters. 

1.0  Btu/hr.  (ft“),  (°  F)  for  vertical  walls 
1.5  Btu/hr.  (ft2),  (°F)  for  the  ceiling 
0.5  Btu/hr.  (ft2),  (°F)  for  the  floor 


D 


Although  there  may  be  some  other  values  and  other  combinations  of 
these  values  that  might  have  resulted  in  a  slightly  better  simulation 
than  those  ucd  in  this  analysis,  these  three  values  can  be  considered 
representative  design  heat  transfer  coefficients  in  the  underground 
cavities. 

For  larger  shelters,  the  one -dimensional  and  compound  model  (M-0) ) 
will  probably  be  adequate  for  calculating  the  shelter  thermal 
environment.  The  complicated  three-dimensional  model,  therefore, 
may  not  be  required  for  the  calculation  simulating  the  14-day 
occupancy  of  many  large  community  shelters.  For  small  sneiters 
(such  as  family  shelters  similar  to  the  NBS  shelter),  however,  it 
la  recommended  that  the  three-dimensional  model  be  uaed  for  the 


accurate  calculation. 
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NOMENCLATURE 


Unless  otherwise  defined  in  the  text,  the  symbols  used  In  this  report  are 
summarized,  as  follows: 


Symbols 

Dimensions 

a 

shelter  dimension  (half  of  the  inside  length) 

ft 

b 

shelter  dimension  (half  of  the  inside  width) 

ft 

C 

P 

specific  heat  of  moist  air 

Btu/lb, °F 

c 

shelter  height 

ft 

d 

shelter  depth,  distance  between  the  eavth  surface 
and  the  ceiling  of  the  shelter 

ft 

G 

flow  rate  of  ventilation  air 

CFM 

h 

g 

ground  surface  heat  transfer  coefficient 

Btu/hr, ft8., °P 

hK°rh 

k 

inner  surface  heat  transfer  coefficient  of  K  th 
exposure 

thermal  conductivity  of  solid 

Btu/hr, ft* , °F 
Btu/hr, ft, °F 

4 

Lewis  Relation  -  _JL_ 

ac 

P 

dimensionless 

Pvs 

saturated  vapor  pressure  of  water  at  temp,  t 

inches  Hg 

Pv 

vapor  pressure  of  water  in  the  air  at  dew  point 
temperature 

inches  Hg 

PB 

barometric  pressure 

inches,  Hg 

%S 

sensible  heat  released  by  the  ventilation  air 

Btu/hr 

qvl 

latent  heat  released  by  the  ventilation  air 

Btu/hr 

qgs 

sensible  heat  generated  in  the  shelter  by 
simulated  occupants 

Btu/hr 

«gl 

latent  heat  generated  in  the  shelter  by 
simulated  occupants 

Btu/hr 

sansibla  haat  released  by  the  shelter  inner  surface 
of  K  th  expoaura 

Btu/hr 

latent  heat  released  by  the  shelter  inner  surface 
cf  K  th  exposure 

Btu/hr 

sensible  heat  generated  in  the  shelter  by  things 
other  than  aimulatad  occupants 

Btu/hr 

* 
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IV 


Symbols 


? 


qml 

latent  heat  generated  in  the  shelter  by  things  . 
other  than  simulated  occupants 

Btu/hr  ; 

J 

1 

qsuk 

SK 

solar  radiation  intensity  at  the  earth ‘iat&faso 

inner  surface  area  of  K  th  exposure 

Btu/kr,fta  | 

:  ft*  j 

t 

temperature 

-  °f  1 

w 

a 

humidity  ratio 

lb  of  water  vapor/ 
lb  dry  air 

W 

8 

humidity  ratio  of  air  saturated  by  water  vapor 

lb  of  water  vapor/j 
lb  of  dry  air 

\  (K" 

1,6)  coordinate  system  used  for  shelter  heat  transfer 

ft 

! 

*K,0  <K 

*  1,6)  coordinates  of  the  system  boundaries 

ft 

*K 

thermal  diffusivity  of  earth  in  the  region  sur¬ 
rounding  K  th  exposure 

ft3 /hr 

Ti 

albedo  of  earth  surface 

dimensionless 

3 

time  coordinate 

hr 

\ 

latent  heat  of  vaporisation  of  water 

Btu/lb 

a 

water  vapor  transfer  coefficient  (lb/ (hr) (ft?) (lb/lb  dry  air) 

P 

density  of  moist  air 

lb/ ft3 

Axk 

finite  difference  length  along 

ft 

A9 

finite  difference  time 

hr 

£ 

summation  symbol 

Subscripts 

Unless  otherwise  stated,  the  following  rule*  of  subscripting 
will  apply  to  all  of  the  variables, 
a  shelter  space  properties 

o  outdoor  air  properties 


v  ventilation  *tr  properties 

c  concrete  property 

8  ground  surface  properties 

«  deep  underground 

K  inneraurface  exposure  index 

K  *  i  **  North 

2  *  South 

3  *  East 

4  *  West 

5  »  Floor 

6  ■»  Roof 

In  scene  cases,  the  subscript  W  is  used  to  denote  the  wall  properties  instead 
of  K  being  1,  2,  3,  and  4* 

The  subscripts  R  and  F  are  employed  in  the  same  manner,  denoting,  respec¬ 
tively,  the  properties  pertaining  to  roof  and  floor  regions. 

S  sensible  heat  property 

s  saturated  air  property 

L  latent  heat  property 


captions  Foa  figures 


Fig.  1.  Schematic  diagram  of  the  matrix  used  for  computer  program  M-(l) . 

Fig.  2.  Schematic  diagram  of  the  heat  conduction  region  used  for  computer 

program  M-(2). 

Fig.  3.  Inner  surface  heat  transfer  coefficient  frequency 
distribution. 

Fig.  4.  Inner  surface  heat  transfer  coeflicient  vs  temperature  difference 
between  the  air  and  inner  surface  of  the  shelter. 

Fig.  5.  Comparison  of  the  calculated  and  observed  shelter  air  tempera¬ 
tures  and  relative  humidities  of  test  1  for  NBS  family  shelter 
(computer  program  M-(l)). 

Fig.  6.  Comparison  of  the  calculated  and  observed  shelter  inner  surface 
temperatures  for  test  1  of  NBS  family  shelter  (computer  program 
M-(l)). 

Fig.  7.  Comparison  of  the  calculated  and  observed  shelter  air  tempera¬ 
tures  and  relative  humidities  for  test  2  of  NBS  family  shelter 
(computer  program  M-(l)). 

Fig.  8.  Comparison  of  the  calculated  and  observed  shelter  air  tempera¬ 
tures  and  relative  humidities  for  test  3  of  NBS  family  shelter 
(computer  program  M-(l)). 

Fig.  9.  Comparison  of  the  calculated  and  observed  shelter  air  tempera¬ 
tures  and  relative  humidities  for  test  3  of  NBS  family  shelter 
(computer  program  M-(2)). 

Fig.  10.  Comparison  of  the  calculated  and  observed  air  temperatures  and 
relative  humidities  for  teat  4  of  NBS  family  aha Iter  (computer 
program  M- (2)). 
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Fig.  11.  Comparison  of  the  calculated  and  observed  shelter  air  tempera* 
tures  and  relative  humidities  for  test  5  of  NBS  family  shelter 
(computer  program  M~(2)). 

Fig.  12.  Comparison  of  the  calculated  and  observed  shelter  air  tempera¬ 
tures  for  summer  condition  test  for  Summerlin  shelter  (computer 
program  M-(2)). 

Fig.  13.  Corparison  of  the  calculated  and  observed  shelter  air  relative 
humidities  for  summer  condition  test  for  Summerlin  shelter 
(computer  program  H-(2)). 

Fig.  14.  Comparison  of  the  calculated  and  observed  earth  temperatures  sur¬ 
rounding  Summerlin  shelter  during  summer  test  conditions  (computer 
program  M-(2)), 

Fig.  15.  Comparison  of  the  calculated  and  observed  shelter  air  tempera¬ 
tures  and  relative  humidities  during  moderate  weather  condition 
test  for  Summerlin  shelter  (computer  program  M-(2)). 

Fig.  16,  Comparison  of  the  calculated  and  observed  shelter*  air  tempera¬ 
tures  for  Broylea  shelter  (computer  program  M- (2)). 

Fig.  17.  Comparison  of  the  calculated  and  observed  shelter  relative 
humidities  for  Broyles  shelter  (computer  program  M«(2)). 

Fig.  IS.  Comparison  of  the  calculated  and  observed  shelter  a tr  tempera - 
tures  for  Kapler  shelter  (computer  program  (2)). 

Fig.  19.  Comparison  of  the  calculated  and  observed  shelter  relative 
humidities  for  Kapler  shelter  (computer  program  M-(2)). 

Fig.  20.  Comparison  of  the  calculated  and  observed  shelter  air  tempera¬ 
tures  and  relative  humidities  for  Read lug  shelter  (computer 

program  M**(2))« 

vin 
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CAPTIOUS  FOR  FIGURES  (cont'd) 

Fig.  21.  Observed  moisture  balance  of  Reeding  shelter. 

Fig.  22.  Comparison  of  the  calculated  and  observed  earth  temperatures 
outside  Reading  shelter  vails  (computer  program  M-(2)). 

Fig.  23.  Comparison  of  the  calculated  and  observed  earth  temperatures 

outside  Reading  shelter  ceiling  and  floor  (computer  program  M-(2)). 

Fig.  24.  Comparison  of  the  calculated  and  observed  shelter  air  tempera¬ 
tures  and  effective  temperatures  of  Ft.  Belvoir  200-man  shelter 
(computer  program  M-(4)). 

Fig.  25.  Comparison  of  the  calculated  and  observed  shelter  air  tempera¬ 
tures  and  effective  temperatures  of  Ft.  Belvoir  1000-man  shelter 
(computer  program  M-(4)). 


COMPARISON  OP  DIGITAL  COMPUTER  SIMULATIONS  OP  THERMAL  ENVIRONMENT 
IN  OCCUPIED  UNDERGROUND  PROTECTIVE  STRUCTURES  WITH  OBSERVED  CONDITIONS 


T.  Kuauda  and  P>  R.  Achenbach 


1 .  INTRODUCTION 


This  report  compares  digital  computer  calculations  of  thermal  environ 
ment  for  underground  protective  shelters  with  observed  conditions  of  tem¬ 
perature  and  humidity.  Several  digital  computer  programs  have  been 
developed  by  the  National  Bureau  of  Standards  for  the  purpose  of  simulating 
the  heat  transfer  of  underground  structures.  These  computer  programs  were 
applied  to  7  shelters,  whose  thermal  environment  under  simulated  conditions 
of  occupation  had  been  observed  experimentally.  Since  the  thermal  environ¬ 
ment  in  underground  protective  structures  may  become  extremely  rnfavOCable , 
particularly  during  the  summer  occupancy  period,  for  lar-je  areas  of  the 
United  States,  the  majority  of  prototype  shelters  mentioned  herein  were 
tested  under  summer  climatic  conditions.  Of  the  7  shelters  whose  thermal 
environments  were  calculated  and  compared  with  experimental  observations, 

12  different  test  conditions  mere  Included,  2  of  which  were  under  moderate 
winter  conditions. 

Analytical  and  experimental  studies  of  various  shelters  have  shown 
thst  ths  temperature  and  humidity  within  the  occupied  underground  shelter 
depend  on  many  parameters,  which  may  be  classified  as  follows i 


e.  Slse  end  shape. 

b.  Physical  end  thermal  properties  of  construction  materiel 


2. 


Site  characteristics. 


a.  Physical  and  tharmal  properties  of  earth  surrounding  the 
shelter. 

c.  Thickness  of  earth  cover. 

d.  Type  of  earth  surface  and  landscape. 

e.  Neighboring  buildings  and  installations. 

3.  Climatic  factors. 

a.  Earth  temperature. 

b.  Pay chrome trie  condition  of  outdoor  air. 

c.  Solar  radiation. 

d.  Precipitation. 

4.  Operational  characteristics. 

a.  Ventilation  rate. 

b.  Psychrometric  condition  of  ventilation  air 

c.  Density  of  Occupancy. 

d.  Activity  of  Occupancy. 

e.  Heat  and  moisture  release  by  equipment  in  shelter. 

f.  Emergency  condition  such  as  sealed  up  or  surface  fire  conditions. 
Testing  of  underground  structures  to  cover  even  a  smell  portion  of  all 

of  the  possible  combinations  of  the  above  parameters  is  a  formidable  and 
expensive  task.  However,  the  number  of  tests  could  be  dreatlcally  decreased, 
and  the'  efficiency  of  testing  improved,  If  Che  effect  of  various  parameters 
In  the  thermal  environment  of  a  shelter  could  be  predicted  by  computation. 

The  mathematical  formulation  of  such  a  computation  should  taka  into  account 
a  majority  of  the  important  parameters  so  the  sensitivity  of  the  overall 
thermal  environment  to  the  several  parameters  could  be  studied  individually 
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or  simultaneously  with  others.  The  mathematical  procedures  should  be 
simple  enough  so  the  computation  time  (or  computer  cost)  would  be  reason 
able.  Finally,  and  most  Important,  the  computed  results  should  be 
reliable . 

Some  previous  computations  of  the  thermal  environment  in  shelters 
have  been  reported  which  take  into  account  simulated  human  metabolism, 
ventilation  effects,  and  heat  transfer  in  the  earth  [1,2,3],  but  there 
have  been  very  few  actual  comparisons  between  calculations  and  observa¬ 
tions  for  a  given  system  over  a  substantial  period  of  time.  It  is  for¬ 
tunate  that  the  observed  results  of  six  of  the  seven  shelters  covered  by 
this  study  were  so  well  documented  by  reports  of  the  University  of 
Forlda  [4] ,  thus  making  possible  comprehensive  comparisons  between  the 

calculated  thermal  environment  and  the  observed  results. 

2.  BASIC  HEAT  TRANSFER  RELATIONS 

Since  the  details  of  the  numerical  technique  employed  for  the  heat 
transfer  of  underground  protective  structures  have  been  reported  pre¬ 
viously  [l],  only  basic  mathematical  formulations  employed  for  all  of 
our  computer  programs  are  given  here: 

2.i.  Shiltti.iii  hilt  bfliaci- 
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This  solar  heat  model  is  essentially  the  same  as  the  sol-air  tempera¬ 
ture  concept  and  ignores  the  direct  radiation  heat  exchange  between  earth 
surface  and  sky,  but  it  does  include  the  effective  radiation  heat  exchange 

between  the  earth  surface  and  ambient  air  by  adjusting  the  value  of  n  . 

8 

According  to  the  last  equation,  however,  the  earth  surface  temperature 
never  becomes  lower  than  the  ait  temperature,  which  is  not  always  the  case. 

2.?.  Psychrometric  calculations. 

Taking  advantage  of  the  large  memory  of  the  high  speed  digital  computer 
of  the  National  Bureau  of  Standards,  all  the  psychrometric  cei  ulatlcns  were 
performed  using  the  thermodynamic  properties  of  moict  air  published  by  Goff 
and  Gratch. 

The  thermodynamic  properties  of  dry  air  and  those  of  saturated  air  at 
one  standard  atmospheric  pressure,  such  as  the  following,  were  tabulated  as 
temperature  functions  by  Goff  and  Gratch  f 5] . 

w„  a  humidity  ratio  of  the  saturated  moist  air  (lb/lb  of  dry  air). 

a  enthalpy  of  dry  air  (Btu/lb  of  dry  air). 
hg  3  enthalpy  of  the  saturated  moist  air  (Btu/lb  of  dry  air), 
hy  a  enthalpy  of  the  wat*r  (Btu/lb  of  water). 

*  volume  of  the  dry  air  (cu  ft/lb  of  dry  air). 

vg  -  volume  of  the  aeturated  moist  air  (cu  ft/lb  of  dry  air). 

*  factor*  related  to  the  relative  humidity  and  degree 
(diamna ton leas)  of  saturation, 

Theae  properties  were  read  into  the  computer  for  the  temperature  range 
frtmi  30  *T  to  120  *F  et  every  one  degree  increment,  except  that  the  humidity 
ratio,  W„,  was  programmed  for  the  temperature  range  from  -20  *F  to  120  *F 
et  every  one  degree  Increment. 


The  following  psychrometric  symbols  and  formulas  are  used  to 
derive  the  desired  properties  from  given  sets  of  properties,  such  as  dry- 
and  wet-bulb  temperatures:  ■  x 


P v  ■  partial  water  vapor  pressure  in  moist  air  (in.  Hr). 

Pp  •  barometric  pressure  (in.  Hg). 

Pvs  ■  partial  water  vapor  pressure  in  saturated  air  (in.  Hg.). 
cp  •  relative  humidity,  as  a  fraction 
p  ■  degree  of  saturation. 

W  *  humidity  ratio  of  moist  air  (lb/lb  of  dry  air). 

V  *  volume  of  moist  air  (cu  ft/lb  of  dry  air). 

Wg*  *  W8  evaluated  at  the  thermodymmic  wet-bulb  temperature 
(lb/lb  of  dry  air). 

hy*  «  evaluated  at  the  thermodynamic  wet -bulb  temperature 
(Btu/lb  of  water), 

h8*  ■  h8  evaluated  at  the  thermodynamic  wet-bulb  temperature 
(Btu/lb  of  dry  air). 


P 

W  -  0.622  ~p  2.7-1 

v 
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s-'  rCi 


h  ■  tk(b8  -  ht)  +  hfi  2,7-6 

h  +  (\i9*  -  V)h^*  -  hs*  2,7-7 

The  last  formula  represents  the  thermodynamic  wet -bulb  temperature  relation 
for  given  h  and  W  of  moist  air. 

Hie  psychrometric  calculations  of  thermal  environment  begin,  usually, 
with  input  data  for  dry-  and  vet-bulb  temperature  of  moist  air.  These  two 
temperature  data  are  sufficient  to  describe  the  complete  thermodynamic 
state  of  the  moist  air  at  one  standard  atmosphere .  The  actual  calculations 
involving  heat-  and  mass-balance  within  a  given  thermal  system  are,  however, 
more  readily  performed  with  the  dry-bulb  temperature  and  the  humidity  ratio, 
as  seen  in  equations  of  shelter  heat  balance.  The  method  used  during  this 
study  to  obtain  the  humidity  ratio  of  moist  air  from  its  dry-  and  wet-bulb 
data  is  described  as  follows : 

The  thermodynamic  wet-bulb  temperature  relation  2.7-7  can  be 
expanded  by  the  use  of  enthalpy  expression  2.7-6. 

~~~~T  Chs )  -  ha(t)]  +  ha(t)  +  [Ws*(t')  -  W]hw(t')  -  hs*(t')  2.7-8 

w8(t) 

In  the  above  expression,  Ws(t),  ha(t),  and  ha(t)  are  thermodynamic 
properties  at  dry-bulb  temperature  t,  while  Wg*(t#),  h8*(t/),  and  hw*(t#) 


are  thermodynamic  properties  evaluated  at  wet-bulb  temperature  , 

Rearranging  the  terms  in  equation  2.7-8,  the  humidity  ratio  W  for 

dry-  ?-A  wet -bulb  temperatures  t  and  t'  can  be  expressed  as 

h,«(e')  -  h.(t)  -  h»»(r')w,«(t') 

*  h.(t>  -  h,(t>  -  **w  Z-' 

The  calculation  of  the  thermodynamic  wet-bulb  temperature  from  t 

given  dry-bulb  temperature  and  humid ity  ratio  is  also  possible  from 
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2,7-8  by  an  iterative  technique.  The  Iterative  technique  found  successful 
during  the  course  of  this  investigation  Is  the  inverse  interpolation 
formula  [6]  of  Newton  applied  to  Goff  and  Gratch  tables. 

2.8.  Effective  temperature  of  shelter  air. 

For  some  of  the  shelter  analyses,  effective  temperatures  have  been 
calculated  from  dry-  and  wet-bulb  temperatures  assuming  an  air,  itlbclty  of 
less  than  20  fpm.  The  effective  temperature  chart  of  the  ASHVB*  has  bean 
stored  In  the  computer  memory  In  a  tabular  format,  and  a  table -searching 
end  lnterpolative  subroutine  used  to  calculate  ther  effective  temperature. 


3.  DESCRIPTION  OF  COMPUTER  PROGRAMS 

3.1.  Computer  mode la. 

Basically  four  different  computer  programs  have  been  developed  during 
this  study.  All  the  programs,  however,  essentially  employ  the  time-iteration 
technique  for  solving  transient  heat  conduction  equations  and  they  are 
designated  as  follows: 

M-fl)  Three-dimensional  model  with  composite  walls  and  with  separate 
initial  temperature  pattern  normal  to  the  six  boundary  surfaces. 

M-(2)  Three-dimensional  model  with  homogeneous  heat  conduction 

medium  having  separate  initial  temperature  patterns  normal  to 
the  six  boundary  surfaces. 

K-(3)  One -dimensional  compound  heat  conduction  model  for  composite 
regions  treating  the  six  exposures  separately. 

M-(4)  One -dimensional  compound  heat  conduction  model  for  six 
composite  regions  assuming  sa  adiabatic  roof. 

Each  program  has  advantages  and  disadvantages,  at  discussed  in  the 


following  pages. 

*  ASHVE  Guide  1950  Chapter  6 
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Thin  program  has  bean  described  In  reference  [2],  and  was 
used  to  evaluate  the  MBS  family  shelter,  The  earth  temperature  field  sur¬ 
rounding  the  shelter  was  divided  Into  6  blocks,  such  as  shown  in  figure  1 . 
This  block  system  enabled  the  program  to  account  £cr  situation*  in  which 
some  well  region (s)  may  he  considerably  different  from  the  othera  in  heat 
transfer  properties  and  earth  tempera ture. 

The  initial  temperature  -in  each  block  was  programmed,  only  in  the 
direction  normal  to  the  wall  surface,  however,  because  the  temperature 
profiles  were  usually  known  only  along  those  directions. 

Temperature  calculations  for  the  entire  earth  region  surrounding  the 
shelter  were  made  with  time-iterative  techniques  on  three-dimensional 
finite  difference  equations.  The  concrete  wall  (including  floor  and  celling) 
temperatures  were  calculated  by  one-dimensional  finite  difference  equations 
separately  for  each  wall,  assuming  that  lateral  temperature  variation  on 
the  interior  surface  for  a  given  wall  could  be  neglected.  Each  inner 
wall-surface  temperature  was  determined  by  surface  heat  balance  equation 
2.2,  including  vapor  condensation  but  excluding  condensate  re-evaporation. 

The  shelter  psychranetric  condition,  dry-bulb  tempera ture,  dew-point  tem¬ 
perature,  and  relative  humidity  were  then  evaluated,  baaed  upon  the  total 
heat  balance  equation  2.1. 

M-(2).  In  this  model,  thermal  properties  or  heat  tranafar  charac¬ 
teristics  around  the  shelter  air  apace  ell  were  assumed  homogeneous.  In 


other  words,  no  distinction  in  thermal  properties  w as  made  from  the  concrete 


wall  to  the  soil,  or  from  one  wall  region  to  the  other,  aa  In  Model  M-(l) 
The  earth  temperature  initialisation  was  parforwad  only  in  the 


direction  normal  to  the  earth  surface.  The  finite  difference  scheme 


employ «d  for  the  three-dimensional  time-iteration  solution  of  the  hone 
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conduction  equation  is  shown  in  figure  2 .  The  homogeneity  assumption  of 
the  heat  conduction  equation  mentioned  makes  it  possible  to  analyse  only 
a  one-quarter  eegment  of  the  shelter,  because  of  the  symmetric  nature  of 
the  entire  system.  The  details  of  this  program  have  been  described  in 
reference  [1].  This  program  was  applied  to  the  NBS  6-man  family  shelter, 
the  Summerlin  shelter,  the  Broyles  shelter,  the  Napier  shelter,  and  the 
Reading  shelter,  for  the  purpose  of  comparing  the  calculated  shelter 
thermal  environment  with  the  experimentally  observed  data. 

M-(3).  The  three-dimensional  effect  on  earth  heat  conduction 
around  the  shelter  becomes  less  and  less  significant  as  the  shelter  size 
increases.  A  one -dimensional  compound  system  was  developed  primarily 
for  large  shelters,  where  a.  major  portion  of  the  heat  flow  is  always 
normal  to  the  shelter  walls,  ceiling,  and  floor.  This  program  is  identi¬ 
cal  with  M-(l) ,  except  that  the  corner  region  of  earth  and  concrete  is 
ignored,  and  the  one -dimensional  finite  difference  equation  mas  used  for 
earth  temperature  determinations .  Comparison  with  the  two  previous 
programs  indicates  that  the  earth  temperature  computation  scheme  was 
drastically  simplified  in  the  model.  The  program  was  applied  to  the 
Summerlin  and  Reading  shelters . 

M-(4).  This  program  was  a  modification  of  M-(3)»  to  simulate 
basement  type  shelters  such  as  the  two  at  Ft.  Belvoir.  The 
heat  transfer  in  the  ceiling  region  of  basesient  shelters  will  be  much 
less  significant  than  that  in  wall  end  floor  regions.  Thus,  in  this 
program,  the  outer  face  of  the  ceiling  layer  was  assused  adiabatic. 


O 


The  calculation  of  the  psychrometric  condition  of  shelter  air  for 
this  program  was  designed  so  that  dry-  and  vet -bulb  and  dew-point  tempera¬ 
ture,  as  well  as  shelter  effective  temperature,  can  be  computed  with  or 
without  the  air  conditioning  system  turned  on.  The  details  of  the  air 
conditioning  calculations  used  in  this  program  are  described  in  the 
appendix.  This  program  was  applied  to  200-man  and  1000-man  shelters  of 
Ft.  Belvoir,  which  were  tested  by  the  University  of  Florida. 

3.2.  Initialization  of  earth  temperature  surrounding  a  protective 
shelter. 

Accurate  heat  transfer  calculations  for  the  early  part  of  an 
occupancy  period  for  these  shelters  are  extremely  difficult,  because  of 
the  uncertainty  about  the  earth  temperature  distribution  around  the 
shelter.  Particularly  for  the  shelters  that  had  been  Installed  wich  only 
a  shallow  earth  cover,  the  temperature  variation  from  roof  region  to  the 
floor  region,  from  one  wall  region  to  another,  and  from  comer  to  flat 
surface  region  was  quite  appre  table,  and  very  complex  to  approximate 
mathematically. 

In  order  to  simulate  this  complex  and  three-dimensional  pattern  of 
the  initial  earth  temperature  around  the  shelter,  assuming  that  such 
three-dimensional  patterns  are  important,  the  earth  temperature  program 
becomes  highly,  and  perhaps  unnecessarily,  complicated.  Therefore,  during 
the  study,  e  simplification  was  made  by  selecting  three  initial  tem¬ 
perature  patterns,  as  follows: 

-1,  Earth  temperature  gradients  are  always  in  a  direction 

normal  to  four  walls,  floor,  and  roof  surfaces  -  block  mode. 
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1-2.  Predominant  earth  temperature  gradient  is  always  along  a 
vertical  path,  from  the  surface  downward  -  vertical  mode. 

1-3.  Earth  temperature  is  virtually  constant  all  around  the 
shelter  -  constant  mode. 

These  three  modes  of  the  initial  earth  temperature  patterns  have 
been  employed  during  the  calculation,  and  presented  in  this  report. 

However,  It  is  important  to  remember  that  the  detailed  earth  tem¬ 
perature  profile  data  surrounding  shelters  seldom  will  be  available, 
even  when  it  is  Important,  for  the  majority  of  actual  shelters.  For 
most  of  the  thermal  environment  calculations  of  underground  structures, 
it  is  usually  assumed  that  the  earth  has  a  single  value  of  temperature. 
This  uncertainty  in  the  Initial  earth  temperatures  is  one  of  the 
sources  of  error  in  predicting  the  thermal  environment  of  shelters, 
particularly  for  the  early  period  of  shelter  occupancy. 

4.  DESCRIPTION  OF  PROTOTYPE  SHELTERS 

Brief  descriptions  of  the  prototype  shelters  used  for  this  analysis 
and  their  characteristics  are  given.  If  the  mathematical  simulation  of 
shelter  thermal  environment  is  to  be  most  effective,  it  is  necessary  to 
secure  accurate  information  regarding  many  parameters  related  to  struc¬ 
tural,  site,  climatic,  and  operational  characteristic*  of  shelters. 

NBg  shelter. 


The  fsailly  else  shelter  in  Washington,  D.  C,,  tested  by  the 
National  Bureau  of  Standard#,  was  constructed  according  to  Bulletin  MP-15 
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of  the  Office  of  Civil  Defense  [3]  with  tone  smell  modifications.  It  was 
a  concrete  wall  shelter  with  external  dimensions  of  12  feet  long,  9  feet 
4  Inches  wide,  and  7  feet  6  inches  high  (not  Including  hatch),  placed  in 
an  excavation  and  covered  with  2  feet  3  inches  of  earth.  The  interior 
dimensions  of  the  shelter  area  were  10  feet  8  inches  long,  8  feet  wide, 
and  6  feet  6  inches  high.  A  2-foot -wide  hatchway  was  Installed  on  the 
north  side  of  the  occupancy  area  and  separated  from  the  8-  by  8-foot  living 
space  by  an  8-inch-thick  concrete  shielding  wall.  The  earth  surface  over 
the  shelter  was  grass  covered  and  partially  in  the  shade  of  neighboring 
trees.  The  soil  around  the  shelter  was  mostly  loam  and  clay;  its  density 
averaged  about  109  lb/ft3  ,  and  its  moisture  content  averaged  about  15 
percent  (dry  weight  basis). 

Ventilation  air  was  controlled  in  a  neighboring  equipment  house  to 
specified  conditions  ami  ducted  into  the  shelter  through  an  inlet  at  the 
mid-height  of  the  shielding  wall  facing  the  occupied  area,  and  the  shelter 
air  was  exhausted  through  an  outlet  located  on  the  opposite  wall  of  the 
hatchway . 

Six  simulated  occupants  (SXMOC)  were  carefully  designed  and  con¬ 
structed  to  produce  sensible  end  ietent  portions  of  metabolic  heat  as 
functions  of  shelter  area  temperature  (deteiU  of  SXMOC  ere  found  in 
reference  [3]).  Four  SXMOC's  bad  a  nominal  haat  output  of  400  Btv/hr, 
and  tha  other  two  had  haat  output*  of  200  and  600  Btu/hr,  reepectivalp . 

Comp let*  psychromatrlc  observations  were  mada  at  the  ventilation 
air  inlat,  tha  axhauat  air  outlet,  and  at  tha  5-foot  lave!  above the 


geometrical  center  of  the  shelter  floor.  Temperature  observations 
obtained  during  this  test  included  earth  temperatures  around  the 
shelter  to  a  maximum  distance  of  4  feet  away  from  the  exterior  surface  of 

the  concrete  wall  for  all  four  walls  and  the  floor.  In  the  shelter  roof  t 

I 

region,  earth  temperature  was  studied  with  four  thermocouples  located  j 

I 

6  inches,  12  Inches,  18  inches,  and  24  Inches  from  the  external  surface  I 

t 

of  the  shelter  ceiling. 

Five  tests  were  conducted  with  variations  in  test  duration,  venti¬ 
lation  rate,  number  of  occupant?,  ventilation  air  condition,  and  sur¬ 
rounding  earth  temperature.  Test  1  was  conducted  with  no  occupants,  with 
42  cfm  of  ventilation  air.  Tests  2,  3,  and  4  were  each  conducted  with 
six  SDfOC's,  and  with  ventilation  air  rates  of  0,  18,  and  42  cfm,  respec¬ 


tively.  Test  5  was  undertaken  to  simulate  '/Inter  conditions  of  occupancy, 
ventilation  air,  and  earth  temperature,  using  six  SZMOC'a  and  18  cfm  venti-  j 

•  '  •  •  j" 

■  i 

’at ion  air.  j 

4.2.  Summerlin  shelter  f4l. 

This  shelter  was  a  welded  steel  structure  located  entirely  below  the  < 


finished  grade  line,  with  a  30- inch  earth  cover  over  the  roof ,  in  a  rural  j' 
■  ■  •'  J 

area  of  Gainesville,  Fla.,  and  thermally  isolated  from  other  be  Id logs,  j 

The  occupancy  ares  dimension*  were  23  feet  7  Inches  long,  and  7  fast  8  j 

.  '  .  V.  •  t 

Inches  wide.  The  shelter  roof  was  arched  over  the  wsll,  sod  the  maximum  j 

ceiling  height  at  the  middle  of  the  arch  wee  7  feet  2 h  inches.  A  4-  by 
3-foot  hatchway  located  st  one  end  of  the  shelter.  The  earth  around 


the  shelter  wes  e  mixture  of  sand  sad  loam,  end  the  earth  cower  was  bare 


at  the  time  of.  the  shelter  environmental  testa.  The  moisture  content  (dry 


weight  basis)  of  t.H*  soil  samples  am  lysed  ranged  from  9  percent  to  19 
percent;  ,  their  dry  density  ranged  from  103  lb/ft*  to  ill  lb/ft*  . 


Ventilation  air  conditions  were  generated  in  an  equipment  trailer 
located  outside  the  shelter,  and  supplied  to  the  shelter  by  a  flexible 
tube  passing  through  the  entrance  hatch.  The  ventilation  air  inlet  was 
near  the  north  end  of  the  shelter,  and  shelter  air  was  exhausted  from  a 
stack  located  at  the  southeast  corner. 

Two  series  of  environmental  tests  were  conducted,  using  18  SIMOC's 
(ref.  3)  of  NBS  400  Btu/hr  type;  one  during  July  1962  with  ventilation 
air  of  700  cfra  and  typical  Florida  summer  outdoor  psychrometric  condi¬ 
tions;  and  the  other  during  April  1963  with  ventilation  rates  of  54  and 
216  cfm  at  August  and  April  psychrometric  conditions .  The  simulated  August 
psychrometric  conditions  used  for  the  tests  were  diurnal  cycles  of  96  ’F 
maximum  dry-bulb  temperature,  80  *F  minimum  dry-bulb  temperature,  and 
/9*8  *P  average  dewpoint  temperature.  The  simulated  April  psychrometric 
conditions  used  for  the  test  were  diurnal  cycles  of  76  *F  maximum  dry- 
bulb  temperature,  61  *F  minimum  dry-bulb  temperature,  and  59.5  *F  average 

dewpoint  temperature. 

4.3.  The  Brovles  shelter [4j  . 

The  Sroyles  Shelter,  at  Gainesville,  Fla.,  wee  so  constructed  that 
a  portion  of  the  shslter  was  below  grade,  and  3  feat  of  earth  was  mounded 
around  tha  above -grade  portion.  The  Interior  of  the  shelter  had  a  ceiling 
height  of  7.33  feet,  end  the  flow  eras  measured  16  b:  7.75  feat,  for  a 
total  crea  of  124  square  feet,  the  floor,  placed  on  e  plastic  «s*abraoe, 
wee  of  waterproof  concrete  reinforced  with  l/2-inch  *t#el  rods  on  12-loeh 
centers.  The  wells  were  hollow  concrete  block*  with  conventional  oerter 
joists,  and  the  cavitte*  were  filled  with  waterproof  cement  as  tike  walla 
were  constructed.  The  roof  we*  e  concrete  slab  reinforced  with  3/g-ineh 
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steal  tods  on  12 -inch  center*.  The  concrete  was  treated  with  a  water¬ 
proofing  material  at  the  time  it  was  mixed.  This  shelter  was  shaded  by 
surrounding  trees,  and  a  heavy  layer  of  sod  and  green  grass  covered  it  and 
the  surrounding  ground .  The  surrounding  earth  waa  a  mixture  of  sand  and 
loam,  whose  dry  density  ranged  from  95  lb/ft3  to  104  lb/ft3  ,  and  moisture 
content  from  1  percent  to  12.5  percent.  ' 

For  most  of  the  thermal  environment  teat  period,  1?  SIUOC's  (ref.  3)  : 
of  NBS  400  Btu/hr  type  were  employed.  The  ventilation  air  was  processed 
outside  the  shelter  to  simulate  diurnal  cyclic  conditions  of  96  3F  maxi¬ 
mum  dry-bulb  temperature,  80  *F  minimum  dry-bulb  temperature^  and  70  *F 
average  dew-point  temperature,  and  was  forced  into  the  shelter  by  a  fan 

■  J 

through  an  air  supply  duct  that  was  laid  along  the  inside  surface  of  the 
longitudinal  shelter  wall,  and  which  had  four  equally  spaced  outlets.  The 
shelter  air  was  exhausted  through  two  outlets  located  at  opposite  «nds  of 
the  shelter. 

In  addition  to  the  regular  psychroiECtric  measurements  of  inlet 


ventilation  air,  shelter  exhaust  air,  and  shelter  air  at  the  geometrical  }■ 

I 

center,  several  ceasurements  were  made  of  the  shelter  inner  wall  surface  i 

I 

temperature  und  surrounding  earth  temperature . 

The  shelter  was  tested  in  four  successive  phases: 

Phase  1:  A  ventilation  rate  of  3  cfm  per  person  was  supplied, 

-f 

and  a  fan-aui-coil  unit  simultaneously  cooled  the  recirculated  air  by  1 

E 

4  gp.i  of  well  water  at  an  inlet  water  temperature  of  71.5  “F,  During  i 

the  first  two  days,  sensible  and  latent  heat  was  removed  from  the  shelter  by 
this  fau-and-coil  unit  at  a  rate  of  6440  Btu/hr.  This  cooling  capacity 
exceeded  the  total  heat  supplied  to  the  SIKDC*s  by  4800  Btu/hr, 


Phage  2 ,  The  we£l-w#tdr  cotl  was  cut  off  and  10  cfa  per 
per a on  of  ventilation  air  was  supplied  for  a  period  of  12  days. 

Phases  3  a  ng  'ft .  Tfegae  two  phases  were  operated  under  the 

^«L»  fm^trnrnm  ~j***L** 

same  conditions  as  phase  except  that  the  ventilation  rata  in  phase  3 
was  6  efm  par  person  (for  4$  hours)  and  in  phase  4  was  3  cfm  per  person 
(for  48  hours). 

4.4,  Saoier  aha Iter  f 4]  , 

The  Napier  shelter  was  a  100-occupant  community  shelter  designed 
for  a  group  of  residents  in  a  subdivision  adjoining  the  city  of 
Gainesville,  Fla.  The  floor  level  was  5  feet  below  grade  and  30  inches 
of  sarth  covered  the  roof.  The  floor  slab  was  wire-mesh  reinforced, 
4-inch-thick  concrete  poured  over  a  waterproof  plastic  membrane.  The 
walls  ware  of  8-inch-thick  hollow  concrete  blocks,  whose  cavities  had 
reinforcing  rods  placed  vertically  through  then  at  selected  intervals 
aod  were  then  filled  with  concrete.  Outside  dimensions  were  20  feet  wide 
and  85  feet  long.  The  interior  floor  area  was  1561  square  ffeet.  Rein¬ 
forced  prestreased  concrete  T-bearas  placed  on  top  of  the  shelter  walls, 
each  in  contact  with  the  beams  parallel  to  it,  formed  its  roof.  The 
surrounding  earth  was  mostly  clay,  whose  dry  density  varied  from  76  lb/ft3 
to  109  lb/ft3,  while  the  moisture  content  ranged  from  3,4  percent  to 
34.2  percent  (dry  weight  basis).  The  earth  surface  over  the  shelter  was 
bare  with  several  patches  of  weeds. 

For  this  test,  100  SIMOC’a  (ref.  3)  of  NBS  400  Btu/hr  type  were 
employed.  The  ventilation  air  was  conditioned  in  the  equipment  trailer 
outside  the  shelter  to  represent  a  diurnal  cycle  of  a  Florida  summer: 

96  *F  maximum  dry-bulb  temperature,  80  *F  minimum  dry-bulb  temperature, 
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and  73  °F  average  dewpoint  temperature .  It  was  forced  into  the  shelter 
at  two  stations  at  the  celling  level  in  the  south  wall,  approximately  8 
feet  from  the  southwest  and  southeast  corners  of  the  shelter.  The 
ventilation  air  was  discharged  in  the  direction  of  the  north  wall;  and 
after  passing  through  the  main  chamber  was  exhausted  through  the  entry 
doorway  and  entry  chamber.  Regular  psychrometric  measurements  were  made 
of  shelter  air,  ventilation  air,  and  exhaust  air.  In  order  to  study  the 
three-dimensional  earth  temperature  profile  during  the  simulated  occu¬ 
pancy  period  of  the  shelter,  numerous  measurements  of  earth  temperature 
were  made  around  the  northwest  and  southwest  corner  regions. 

The  test  was  conducted  in  six  phases; 

Phase  1.  Ventilation  rate  was  maintained  at  3  cfm  per  person 
for  48  hours,  during  which  a  well-water  coil  using  12  gpm  of  72.2  °F  water 
was  In  operation.  The  measured  total  cooling  capacity  of  the  well-water 
coil  during  the  phase  averaged  20,700  Btu/hr,  amounting  to  51,6  percent 
of  the  total  heat  released  by  SIMOC's. 

Phase  2 .  Ventilation  air  rate  was  Increased  to  6  cfm  per 
person  and  the  well  water  coil  was  shut  off  for  this  period  of  97  hours. 

Phase  3.  With  all  other  conditions  being  identical  with  those 
of  phase  2,  30  NBS  SIMOC's  were  replaced  by  one  HASS  SIMOC  of  MRD  [4]  for 
42  hours.  (One  MASS  SIMOC  has  an  adjustable  output  from  1-40  SIMOC's  of 
NBS  400  Btu/hr  type.) 

Phase  4.  The  MASS  SIMOC  was  replaced  by  30  NBS  SIMOC's  for 

51  hours. 


c 


Phase  5.  Ventilation  air  rate  was  increased  to  8.05  cfm  per 
person  with  100  NBS  SIMOC's  for  116  hours. 


Phase  6.  The  total  ventilation  rate  remained  at  5u5  cfm 
and  100  KBS  SXMDC's  were  used  together  with  1  MASS  SBSOC  for  36  hours , 
representing  a  total  of  140  occupants. 

4. 5,  Reading  shelter  [4]. 

This  was  a  community  shelter  located  In  a  park  owned  by  the  city 
of  Reading,  Pa.  It  was  constructed  in  a  hillside  to  take  advantage  of  a 
thick  earth  cover,  and  was  basically  a  rectangular  parallepiped  56  feet 
long,  17  feet  4  inches  wide,  with  a  celling  height  of  7  feet  8  Inches. 

The  ends  of  the  shelter  were  connected  to  separate  tunnel-like  entry 
corridors  so  that  two  right  angle  turns  were  formed  in  each  of  these 
passageways.  The  floor  was  4  feet  6  Inches  below  the  grade  level  that 
existed  prior  to  construction,  and  the  roof  and  all  aidea  had  a  minimum 
earth  covering  of  30  Inches.  All  bearing  walla,  the  roof,  and  the 
floor  were  constructed  of  concrete  reinforced  with  steel  bars.  Water¬ 
proofing  was  applied  to  the  external  surfaces  of  the  shelter  during 
construction.  There  were  many  Internal  partition  walla  constructed  of 
hollow  concrete  blocks  with  mortar  filled  voids,  for  rooms  of  various 
purposes,  such  as  storage,  first  aid,  mechanical  equipment,  and  lavatory. 
The  surrounding  earth  waa  sandy  loam  of  dry  density  between  89  lb/ft3  and 
121  lb/ft°,  and  moisture  content  between  10  percent  anu  22  percent  (dry 
weight  basis).  The  earth  surface  at  the  time  of  testing  was  graaa  under 
snow.  Ventilation  air  artificially  created  In  an  equipment  trailer 
outside  the  shelter  was  carried  to  the  shelter  through  a  duct  connected 
to  e  stack  that  under  normal  shelter  operation  would  be  utilised  as  sn 
exhaust  stack.  A  temporary  distribution  duct  for  ventilation  sir  supply 
was  installed  at  the  celling  level  along  the  west  wall  of  tha  shelter. 
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Three  air  outlets  were  installed  In  the  temporary  duct  at  equally  spaced 
intervals .  The  air  was  exhausted  through  a  stack  located  In  the  center 
of  the  celling  of  the  equipment  room,  which  was  almost  longitudinally  at 
the  opposite  end  from  the  air  supply  duct. 

The  ventilation  air  conditions  used  for  this  test  varied  considerably! 
however,  the  dry-bulb  temperature  was  generally  maintained  between  end 
40  °F,  with  an  average  dewpoint  temperature  of  JQ  °F.  The  simulated 
occupants  used  were  two  MASS  SIMOC's  (ref  4),  which  were  capable  of 
producing  sensible  and  latent  metabolic  heat  equivalent  to  120  sedentary 
adults. 

Psychrometrlc  observation  stations  were  located  at  the  ventilation 
air  inlet,  geometric  center  of  the  shelter,  and  exhaust  air  outlet. 
Thermocouples  were  used  to  measure  interior  surface  temperatures  and 
several  ground  temperatures  extending  downward  45  inches  distance  from  the 
top  of  the  floor  slab  and  horizontally  from  the  outside  surface  of  the  west 
wall.  The  roof  region  ground  temperature  was  also  measured  at  several 
distances  from  the  outside  surface. 

The  test  was  divided  Into  seven  phases,  described  as  follows; 

Phase  1.  Ventilation  at  the  rate  of  150  c£m  was  supplied 
to  the  shelter  in  the  manner  deacribed  previously,  and  the  two  MASS  SXNOC's 
were  adjuated  to  deliver  heat  and  moisture  equivalent  to  a  total  of  50 
sedentary  adulta.  This  phase  continued  from  February  26  to  March  4,  1963. 

Phase  2.  Ventilation  rate  was  than  decreased  to  75  cfm,  or 
1.5  cfm  per  person,  with  operating  conditions  identical  to  those  of 
phase  1.  This  phase  lasted  approximately  5  days. 
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Phage  3.  The  mass  SIMOC ' s  were  readjusted  to  produce  total 
equivalent  metabolic  heat  for  100  occupants.  The  total  ventilation  air 
rate  was  150  cfm  during  this  phase  of  study,  lasting  three  days. 

Phase  4.  During  the  three-day  period  that  followed  phase  3, 
from  March  12  to  14,  shelter  occupancy  was  50  simulated  occupants.  The 


ventilation  air  at  a  rate  of  8J1  cfm  was  supplied  by  the  blower  that 
originally  had  been  installed  in  the  shelter  as  a  part  of  the  permanent 

facilities.  The  shelter’s  own  air  handling  systems  provided  the  distri¬ 
bution  and  exhaust  of  the  shelter  air.  The  purpose  of  this  phase  of 
the  test  was  to  determine  if  the  original  equipment  for  supplying  air, 
and  the  air  distribution  and  venting  systems  were  adequate. 

Phase  5.  During  this  phase,  ventilation  air  at  the  rate  of 
718  cfm  was  supplied  for  50  simulated  occupants.  The  air  distribution 
and  exhaust  systems  were  changed  during  this  phase,  details  of  which  are 
not  important  for  the  purpose  of  this  report. 

Phases  6  &  7.  These  phases  were  conducted  for  sealed-up 
conditions  without  ventilation  for  50  simulated  occupants  for  24  hours, 
and  100  cimulated  occupants  for  58  hours,  respectively. 

4.6.  Ft.  Belvoir  200-man  shelter. 

This  was  an  experimental  200-man  shelter  designed  end  built  by  the 
Protective  Structures  Development  Center  at  Ft.  Belvoir,  Va.  it  waa  a 
two-etory  reinforced  concrete  structure  with  one  story  below  ground  level 
Because  the  building  had  been  occupied  by  personnel  of  the  Protective 
Structures  Development  Center  prior  to  the  environmental  testing,  the 
temperature  and  humidity  had  been  comfort  conditioned.  Only  the  base¬ 
ment  area  was  used  for  the  shelter  environment  test,  while  the  upper 
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floor  area  was  heated  to  simulate  an  adiabatic  roof  condition.  The 
shelter  walls  were  10  inches  thick,  the  floor  slab  6  Inches  thick,  and 
the  ceiling-floor  slab  8  inches  thick.  There  was  a  12 -inch-thick  fill 
of  gravel  between  the  external,  surface  of  concrete  of  the  basement  and 
the  earth  on  all  sides,  except  the  north  end.  The  earth  around  the 
shelter  was  a  mixture  of  sand  and  clay,  its  dry  density  varying  from 
85  lb/ft3  to  120  lb/ft3  and  its  moisture  content  from  9  percent  to  23 
percent  (dry  weight  basis).  The  earth  surface  was  composed  of  sod  on 
all  sides  of  the  shelter,  except  for  a  portion  of  the  north  end  where 
a  bituminous  concrete  driveway  was  located.  Interior  dimensions  were 
37  feet  2  inches  by  37  feet  2  inches,  with  a  celling  height  of  9  feet  6 
inches.  Excluding  a  first  aid  room  and  stairwell,  a  net  usable  floor 
area  of  1032  square  feet  was  available  for  100  simulated  occupants. 

There  were  two  family-type  basement  shelters  that  were  built  for  display 
purpose*  in  the  test  room,  one  of  sand-filled  concrete  block  walls  and 
the  other  a  triangular-shaped  wooden  lean-to  filled  with  sand  and  sand 
bags.  They  iray  have  had  a  considerable  effect  in  absorbing  heat  during 
the  first  few  days  of  the  environmental  test. 

The  temperature-  and  humidity-controlled  ventilation  air  (dry-bulb 
cycles  93  *F  ~  76  *F  and  wet-bulb  78  * F  —  74  *F  were  in  phase  with 
maximum  at  2  p.m.  and  minimum  at  4  a.m.)  was  produced  in  the  equipment 
trailer  parked  outside  the  shelter  and  introduced  into  the  shelter  from 
celling  level  near  the  centar  of  the  occupied  space.  The  exhaust  outlet 
was  also  in  the  celling,  approximately  12  feet  from  the  air  inlet.  The 
thermal  environment  of  the  lover  shelter  space  was  measured  with  100 
SIM0C1*  of  the  MBS  type  (ref.  J). 
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The  test  was  conducted  in  five  phases ,  in  which  the  ventilation 
air  rates  were  3  cfm  per  person  for  the  first  phase  (6  days),  6  efts  per 
person  for  the  second  phase  (2  days),  9  efts  per  parson  for  the  thir4 
phase  (3  days),  18  cfm  per  person  for  the  fourth  phase  (2  days),  and  27 
cfm  per  person  for  the  fifth  phase  (4%  days). 

Psychrometric  measurements  were  made  at  the  ventilation  air 
inlet,  air  exhaust,  east  geometric  center,  and  west  geometric  center  of 
the  shelter.  Temperature  measurements  were  obtained  of  various  inner 
surfaces,  together  with  surrounding  earth  temperatures  to  a  distance  of 
4  feet  from  the  external  surface  of  the  concrete  structure. 

4.7.  Ft.  Belvoir  1000-man  shelter. 

This  shelter  was  built  by  the  Protective  Structures  Development 
Center  at  Ft.  Belvoir,  Va.  It  was  a  two-rtory  reinforced  concrete 
building  with  one  story  below  ground  level.  Prior  to  testing  the  building 
had  been  occupied  by  personnel  of  the  Protective  Structures  Development 
Center,  and  its  thermal  environment  had  been  controlled  at  a  comfort  air 
condition,  by  a  central  heating  and  cooling  plant.  Only  the  basement 
area  was  used  for  test,  while  the  upper  floor  level  was  heated  to  simu¬ 
late  an  adiabatic  roof  condition.  The  usable  floor  area  of  the  basement 
vas  5400  ft?  and  was  given  e  simulated  occupancy  of  11  MASS  SXMOC’s  (ref. 
4).  The  shelter  walls  were  10  inches  thick,  the  floor  6  inches  thick,  and 
the  calling-floor  slab  inches  thick.  The  interior  dimensions  were  74 
feet  4  inches  by  74  feet  4  inches,  with  e  celling  height  of  9  feet  6 
inches.  The  surrounding  earth  was  sandy  clay  and  rangad  in  moisture 
content  from  10  to  20  percent  (dry  weight  bee is),  end  in  dry  density 
from  86  lb/ft9  to  104  lb/ft3  dapanding  upon  tha  dapth,  as  veil  as  upon 
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the  location  around  the  shelter.  During  the  backfilling  and  Bounding 
operations,  a  12 -inch -thick  fill  of  gravel  was  placed  adjacent  to  all 
the  basement  walls. 

The  ventilation  _-ir  was  conditioned  by  the  equipment 
trailer  to  yield  a  loca*.  design  diurnal  cycle  of  dry-bulb  temperature, 

94  “F  maximum,  76  *F  minimum,  and  wet-bulb  temperature,  78  *f  maximum, 

74  ®F  minimum.  The  maximum  and  minimum  conditions  of  dry-  and  wet -bulb 
temperatures  were  both  at  2  p.m.  and  4  a.m.,  respectively.  The  ventila¬ 
tion  air  inlet  and  exhaust  air  outlet  were  both  near  a  corner  of  the 
space  and  separated  from  each  other  by  approximately  12  feat  for  most  of 
the  test  period. 

The  test  was  divided  into  four  major  phases,  excluding  an  initial 
short  period  of  purging.  Phase  1  was  a  sealed-up-condition  test  of 
approximately  8  hours  without  ventilation,  whereas  the  ventilation  ratea 
for  remaining  phaaes  were  varied  from  16.8  cfm  per  person  for  days, 
to  5  cfm  per  person  for  1  day.  The  aheltor  thermal  envlronaient  was 
measured  by  dry-  and  wet-bulb  thermometers  at  the  exhauat  duct,  the 
northeast  area,  the  southwest  area,  and  the  geometric  center  of  the 
shelter,  and  numerous  thermocouple  temperature  readings  wore  taken  for 
many  parts  of  the  wall  surface  and  the  earth  region  extending  outward 
normally  to  the  walls  and  floor  to  a  distance  of  4  feet  from  the 
exterior  surfaces. 

5.  SHELTER  CALCULATIONS 

As  discustad  in  section  2  of  this  report,  the  calculated  shelter 
thermal  environment  la  affected  by  several  paraaaeters,  such  as  thermal 
conductivity  and  thermal  diffuslvlty  of  surrounding  walls  snd  earth. 
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and  heat  transfer  coefficients  along  the  inner  surfaces  of  the  shelter  and 
the  ground  surface.  Based  upon  the  characteristics  of  prptotypa  Shelters 
and  their  test  conditions  described  in  the  previous  section,  table  l  it  pre¬ 
pared  to  summarise  the  parameters  used  for  the  shelter  calculations.  However, 
many  of  the  heat  transfer  parameters  selected  for  the  computation  cannot  be 
too  precise.  The  following  considerations  were  given  for  assigning  numerical 
values  to  the  parameters  listed  in  table  1. 

5.1.  Physical  dimensions  (aism  and  heat  transfer  area). 

The  else  of  the  shelter  is  expressed  in  the  overall  Internal 
dimensions,  which  are  the  overall  external  dimensions  leas  the  thickness 
of  walla,  ceiling,  and  floor.  The  heat  transfer  area  used  for  die  com¬ 
putation  was  calculated  from  the  internal  dimension  only,  thus  ignoring 
the  complex  psttern  of  partitioning  and,  consequently,  the  heat  absorp¬ 
tion  by  partition  walla  or  columns. 

5-2-  Siaa.tf.f.d , .°r-SuP»nt,« 

The  simulated  occupants  (SZMOC's)  used  for  the  experiment  were 
designed  to  simulate  the  heat  output  of  human  bodies  snd  each  genanstad  a 
total  hast  of  400  Btu/hr,  regardless  of  temperature.  The  sensible  heat 
component  wee  regulated  according  to  the  shelter  air  temperature  in  the 
manner  described  in  section  2.  The  number  of  active  SXMOC'm  wee  verted 
during  some  tests  (Kapler  snd  Reading  shelters)  by  turning  the  power 
off  end  on  to  s  pert  of  them.  For  this  reason,  the  area  par  person 
data  shown  in  table  1  for  these  tests  are  not  constant. 

In  order  to  simulate  the  change  in  number  of  SXMOC’s  during  seme 
of  the  teste,  the  KB#  computer  program  was  comstrmstsd  to  taka  « is  number 
of  occupants  as  s  time  varlr.*__  A-  stead  of  a  constant , 


il  cetMfcttXvU. 

b«a  (  * 
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5.3.  Ventilation  air. 

Observed  hourly  data  on  dry-bulb  and  devpoint  temperatures  or 
dry-bulb  and  vet-bulb  temperatures  of  the  supply  air  uuct  were  used  as 
the  input  data,  either  in  the  2-hour  step  or  the  4-hour  step  calculations. 
Since  the  University  of  Plorida  varied  the  ventilation  air  rate  during 
their  teats,  the  computer  program  was  resigned  to  accept  the  ventilation 
rate  as  time  dependent  input  data.  These  variations  of  the  number  of 
SIHOC’s  and  ventilation  rates  during  the  tests  are  the  reason  that  the 
per  capita  ventilation  rate  Indicated  in  table  i  was  not  a  constant  for 
some  of  the  tests . 

The  ventilation  air  temperature  and  humidity  conditions  were 
varied  with  respect  to  time,  closely  following  a  prescribed  diurnal 
cyclic  pattern  in  most  cases.  The  entries  in  table  1  show  approximate 
ranges  of  the  temperature  levels  employed  during  the  test,  whereas  actual 

hourly  vt  uea  were  used  In  the  calculations. 

5.4,  Thermal  properties  of  earth  and  wall. 

For  all  of  the  prototype  shelters,  Soil  samples  vere  taken  from 
several  re presents tiv*  spots  around  the  shelter  at  selected  depths  and 
were  analysed  with  respect  to  soil  classification,  dry  density-  and  water 
content  (dry  weight  basis);  their  typical  characteristics  are  shown  la 
table  1,  The  thermal  conductivities  and  diffusivitias  of  various  types 
of  toils  are  usually  presented  as  functions  o£  moisture  content  [5,6]. 

Such  charts  were  consulted  in  arriving  at  the  values  of  therm*  1  conduc¬ 
tivity*  k,  and  thermal  diffusivity,  or,  in  table  1, 

For  the  HBS  shelter,  the  therms!  diffuaivity  value  of  0.026  ft3 /hr. 
for  tests  l  end  2  was  eat lasted  from  a  phase  angle  shift  of  the  earth 
temperature  cycles  at  two  diffarant  depth*.  The  diffuaivity  values  cf 
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0.022  ft2 /hr  for  teats  3  and  4,  and  9.03  ft? /hr  for  test  5  ware  evaluated 
by  a  numerical  technique  similar  to  that  reported  by  Beck  [7],  The  method 
is  basically  a  reversed  use  of  finite  difference  solutions  of  the  heat 
conduction  equation  for  a  semi-infinite  solid.  For  the  case  of  composite 
wall  models,  the  thermal  diffusivity  and  conductivity  of  concrete  were 
assumed  as  0.036  ft3 /hr  and  l.i  Btu/hr  fta (deg  F/ft),  respectively. 

5.5.  Surface  Heat  ..Transfer  Coefficients 
The  surface  heat  transfer  coefficient  consists  of  a  radiative  por¬ 
tion  and  a  convective  portion.  For  the  heat  transfer  at  the  interior  shel¬ 
ter  surfaces,  the  radiation  portion  plays  a  predominant  role,  since  the  air 
velocity  over  the  shelter  inner  surface  usually  is  very  small.  In  addi¬ 
tion,  the  ..urface  beet  transfer  coefficient  along  any  one  of  the  vertical 
walls  at  a  given  time  would  vary  considerably  from  the  bottom  to  the  top, 
due  to  the  vary.ug  nature,  of  layer  pattern,  as  well  as  that  of  the  radia¬ 
tion  heat  exchange  geometry  and  air  velocity,  and  also  because  of  the 
difference  in  local  temperature  distribution.  However,  it  is  probable 
that  the  local  variation  of  the  surface  heat  transfer  coefficient  along  a 
given  surface  may  be  of  the  same  order  of  magnitude  as  the  convective  heat 
transfer  coefficient  itself.  Several  attempts  were  made  during  this  study  to 
obtain  the  surface  heat  transfer  coefficients  for  the  experimental  obser¬ 
vations  of  the  shelter  wall  heat  conduction.  Since  the  NBS  test  sat  Iter 
was  equipped  with  a  heat  flow  meter  at  the  geometric  center  of  each  of 
all  the  inner  surfaces,  and  since  the  inner  surface  temperatures  and 
shelter  air  temperatures  were  simultaneously  measured,  it  was  possible 
to  calculate  the  heat  transfer  coefficients.  However,  the  accuracy  of 
this  procedure  is  questionable,  because  (a)  the  heat  flow  meter  reading 


at  the  geometric  center  of  an  inner  surface  did  not  necessarily  yield 
an  average  heat  transfer  coefficient  for  the  entire  surface,  because  of 
the  local  variation  of  heat  flow,  and  (b)  where  condensation  of  shelter 
air  moisture  was  taking  place  siaultane/'  ;sly,  the  heat  flow  meter  would 
read  the  total  heat  flux,  which  is  not  proportional  to  the  temperature 
difference  between  the  surface  and  air.  In  fact,  these  difficulties 
were  realized  in  the  MBS  shelter,  since  the  ratio  agreement  between  the 
shelter  total  heat  conduction  estimated  by  the  heat  balance,  and  the 
heat  conduction  based  upon  the  heat  flow  meter  ranged  from  52  percent 
to  108  percent  in  various  tests  [2]. 

Nevertheless,  the  values  of  heat  transfer  coefficients  listed  for 
the  NBS  shelter  tests  1,  2,  3,  and  4,  and  the  Summerlin  shelter  test-  weife 
estimated  from  the  test  3  data  of  the  heat  flow  meter  readings,  adjusted 
by  the  vapor  transfer  due  to  the  difference  of  air  humidity  ratio  between 
the  air  and  the  wet  surface,  and  the  Lewis  relation  of  heat  and  mass  transfer. 

A  low  value  for  the  floor  coefficient  was  observed.  Thia  result 
was  to  be  expected  since  the  floor  surface  was  typically  colder  than  the 
air  immediately  above  it  and  the  downward  connection  heat  transfer  rate  would 
be  very  low  under  these  conditions.  Moreover  the  effective  heat  transfer 
area  of  the  floor  was  considerably  reduced  by  the  presence  of  simulated, 
occupants.  The  same  value*  for  the  surface  coefficient  were  applied  to 
the  teat  5  condition  of  the  NBS  shelter,  but  the  agreement  between  the 
observed  and  calculated  inner  surface  temperatures  was  rather  poor 
when  thia  low  value  of  the  surfac*  coefficient  was  used. 
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Several  other  combinations  of  the  surface  heat  transfer  coefficients 
for  KBS  test  5  were  tried  and  the  results  are  summarised  in  table  2,  which 
will  be  elaborated  later  in  this  section.  For  the  many  other  prototype 
shelters,  heat  flow  meter  readings  were  not  included,  or  were  not  usable, 
so  that  the  detailed  heat  transfer  analysis  on  the  surface  by  surface 
basis  was  discarded.  However,  the  average  heat  transfer  coefficients 
for  the  entire  inner  surface  were  obtained  from  the  sensible  heat 
balance  calculation  for  the  Summerlin  shelter  during  the  moderate  weather 
condition  test,  the  KBS  test  3»  the  NBS  test  5,  the  Ft.  Belvior  1000-man 
shelter,  and  the  Kapler  shelter. 

As  indicated  in  figures  3  and  4,  the  result  of  this  analysis  on 
the  overall  sensible  inner  surface  heat  transfer  coefficient  shows  a 
greatly  fluctuating  pattern.  The  accuracy  of  the  calculated  values  of  the 
surface  heat  transfer  coefficient  is  inherently  related  to  the  accuracy  of 
measuring  air  and  surface  temperatures.  Since  radiation  would  usually  be 
present  and  wetted  surfaces  are  sometimes  involved^  the  error  in  temperature 
measurement  could  easily  be  a  signigicant  part  of  the  observed  temperature 
difference  for  differences  of  2  °F  or  less.  A  significant  trend  can  be 
observed,  however,  from  Figure  5  and  4  that  the  overall  inner  surface 
sensible  beat  transfer  coefficient  seems  to  Increase  as  the  temperature 

difference  between  the  air  and  inner  surface  decreases.  ’  The  majority  of  the! 

i 

sensible  heat  transfer  coefficients  are  in  the  neighborhood  of  1.0,  except 
for  the  Napier  shelter.  Nevertheless,  the  iuner  surface  heat  transfer 
coefficients  selected  for  the  Froyles  and  Napier  shelters  of  table  1  are 
somewhat  arbitrary,  whereas  those  given  to  the  Reading  and  Ft.  Belvior 
shelters  are  estimated  from  ref,  8,  which  was  net  available  prior  to  | 

the  last  three  shelter  calculations.  The  possible  effects  of  the  various  - 


Table  2.  EFFECT  OF  INNER  SURFACE  HEAT  CONDUCTANCE 
FOR  N3S  SHELTER  TEST  5 
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sirfsee  heat  transfer  coefficients  upon  calculations  of  the  shelter 
thermal  environment  can  be  seen  in  table  2,  which  illustrates  a  result 
of  sensitivity  analyses  of  inner  surface  heat  transfer  coefficients  using 
computer  program  M-(2)  on  the  NBS  shelter. 

In  table  2,  the  calculated  shelter  surface  temperatures,  shelter 
air  temperatures,  and  shelter  air  relative  humidities  are  compared  with 
the  observed  data  at  the  end  of  the  seventh  and  fourteenth  day  on  test 
5  of  the  NBS  shelter.  The  symbol  WG  in  table  2  is  zero  when  solar  heat 
input  to  the  ground  surface  above  the  shelter  is  Ignored,  while  it  is 
one  when  the  entire  solar  heat  effect  is  assumed  absorbed  by  the  ground 
surface.  If  they  are  compared  at  an  identical  condition,  the  calculated 
shelter  air  temperature  at  the  fourteenth  day  became  approximately  3  F 
higher  when  solar  heat  input  was  considered  than  when  it  was  ignored. 

For  the  fourteenth  day  results,  the  interior  surface  temperature 
agreement  between  the  observed  and  calculated  is  better  with  rather  high 
values  of  surface  conductances,  while  better  air  temperature  agreement  is 
obtained  with  relatively  low  values  of  surface  conductance.  The  combina¬ 
tion  of  h  =  1.21,  h„  =  1.4*.  and  h_  =  0.57  is  obtained  from  reference 
w  R  F 

[lOj  on  the  basis  that  the  simulated  occupants  obstruct  each  wall  from 
seeing  each  other.  As  far  as  the  shelter  air  temperature  is  concerned, 
a  combination  of  =  1,0,  hR  =  1.5,  and  hR  -  0-3  yielded  the  best 

agreement,  which  was  also  true  for  NBS  tests  5  and  A. 

5.6.  Outdoor  conditions. 

Outdoor  air  conditions  during  the  prototype  shelter  test  periods 

were  not  the  same  as  the  ventilation  air  conditions,  as  would  usually 

be  the  case  for  actual  non-air-conditioned  shelters.  The  ventilation  air 

conditions  were  selected  and  programmed  according  Co  certain  cliniactic 

criteria  to  simulate  typical  operating  conditions  of  a  shelter,  regardless 
of  the  actual  climatic  condition  during  the  test  period. 
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The  cutd >or  air  temperatures  used  for  the  calculations  in  this 
report  were  recorded  separately  during  the  test,  averages  of  which  are 
shown  in  table  1. 

The  NBS  tests  included  observations  of  solar  energy  at  the  shelter 
roof  surface,  but  the  tests  conducted  by  the  University  of  Florida  did 
not  have  these  data.  Therefore,  the  inclusion  of  solar  energy  data  for 
the  shelters  tested  by  the  University  of  Florida  was  accomplished  by 
using  data  supplied  by  Flanigan  [ll]. 

5.7.  Initial  earth  temperature. 

As  described  in  Section  3,  the  computer  model  M-(l)  incorporated 
an  earth  temperature  distribution  normal  to  the  surface  in  the  six 
surrounding  blocks.  The  NBS  test  shelter  included  these  observations  and 
they  were  used  for  the  M-(l)  calculation  of  this  shelter.  The  M-(2) 
calculation  provided  for  only  a  vertical  or  depthwise  distribution  of 
initial  earth  temperature.  The  observed  data  on  roof  region  temperatures, 
average  of  the  wall  regions,  and  observed  distribution  of  floor  region 
were  used  to  arrive  at  a  depthwise  distribution  of  the  earth  temperature. 
Average  earth  temperatures  of  the  wall,  roof,  and  floor  regions  were  used 
for  the  calculations  employing  M-(3)  and  M-(4)  models.  Table  1  shows  the 
overall  average  values  of  initial  earth  temperature  which  were  used  for 
the  calculations  in  these  latter  two  models. 

5.8..  Comparison  between  the  computed  and  observed  shelter 
thermal  environments 

Figures  5  through  2*  represent  some  of  the  results  obtained  by 
computer  analysis  of  thermal  environment,  together  with  the  observed 
data.  Figures  5,  6,  7,  and  8  compare  the  calculated  shelter  thermal 
environments  with  the  observed  data  for  ttoats  1,  2,  and  3  on  tha  NBS 


PROPERTIES  INSIDE  UNDERGROUND  SHELTER 
CF*'  VENTILATING  AIR  -  NO  OCCUPANTS 
TEST  NO. I 


CALENDAR  DATE  —  AUGUST  ,  1959 

Fig.  S  Comparison  of  the  calculated  and  observed  shelter  air  temperatures 
and  relative  humidity  of  test  1  for  NBS  family  shelter  (congmter 
program  M-(l)) . 
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relative  humidity  for  teat  3  of  KBS  family  shelter  (computer 
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family  shelter.  The  M- ( 1 )  computer  program  was  used  for  these  compari¬ 
sons.  Computed  shelter  air  temperatures  agreed  with  the  observed  within 
approximately  2  °F  for  the  entire  period  of  shelter  occupancy.  Agreement 
between  computed  and  observed  relative  humidities  was  good  for  tests  1 
and  3,  and  somewhat  poorer  for  test  2  using  the  M-(l)  program.  Figure  7 
represents  a  sealed-up  condition  with  6  occupants,  as  indicated  in  table  1. 
A  very  good  agreement  in  calculated  and  observed  air  temperature  for  the 
first  4-day  period  is  shown  in  figure  7.  However,  the  calculated  relative 
humidity  was  as  much  as  6  percent  lower  than  the  observed  during  the 
second  day  of  test,  and  did  not  quite  attain  the  saturated  condition  that 
was  observed  during  the  test,  even  after  48  hours  of  sealed-up  condition. 

Figures  9,  10,  and  11  compare  temperatures  and  humidities  obtained 
by  computer  program  M-(2)  with  the  observed  conditions  for  NBS  shelter 
tests  3,  4,  and  5.  Figures  8  and  9,  both  for  test  3  condition,  show  that 
computer  program  M-(l),  the  temperature-block  model,  results  in  better 
agreement  with  the  observed  shelter  temper* tuies  for  the  first  two  days 
than  that  by  program  M-(2),  the  homogeneous -earth  model.  This  was 
expected,  because  program  M-(l),  as  explained  before,  is  considerably 
more  elabo  ate  in  accounting  for  a  complex  nature  of  initial  earth  tem¬ 
perature  distribution  around  the  shelter  than  program  M-(2),  and  also 
takes  into  account  the  differences  ii.  thermal  properties  of  the  concrete 
and  the  earth. 

A.  Howavar ,  tha  agreement  between  the  calculated  and  observed  air 
temperature  at  the  5-foot  level  in  the  shelter  was  somewhat  batter  using 
the  M-(2)  program  than  the  M-(l)  program  during  the  last  4  or  5  daya  of 


taat  3. 
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The  agreement  nf  the  computed  shelter  air  temperature  with  the 
observed  condition  for  test  4  of  the  N3S  shelter  became  poorer  toward  the 
end  of  the  test.  The  lowering  of  the  computed  shelter  temperature  during 
the  second  week  of  Che  test  was  caused  primarily  by  a  decrease  In  outdoor 
temperature.  This  temperature  decrease  averaged  approximately  15  degrees, 
beginning  on  October  12  and  continuing  to  the  end  of  the  test  period.  The 
detailed  report  [3]  on  these  tests  indicates  a  steady  and  significant 
decline  of  earth  temperature  during  this  period.  The  curves  on  figure  10 
shows  that  the  observed  shelter  air  temperature  was  not  affected  as  sig¬ 
nificantly  by  this  cool  spell  as  the  computer  model  indicated  that  it 
wou Id  be . 

The  sudden  and  irregular  drop  of  the  calculated  relative  humidity 
for  the  NBS  test  5  condition  shown  in  figure  11  is  not  reflected  by  the 
observed  data.  It  Is  probable  that  the  relative  humidity  in  the  shelter 
was  sustained  at  a  high  level  by  drying  of  the  shelter  walls.  The  mois¬ 
ture  balance  between  the  supply  and  exhaust  air  indicated  that  this 
evaporation  amounted  to  approximately  2  to  2.5  lb/day  during  this  winter 
test  condition  where  extremely  dry  ventilation  air  was  employed.  As 
indicated  earlier,  the  computer  program  developed  for  this  analysis  had 
no  provision  for  taking  this  drying  process  Into  consideration.  Table  2 
Indicates  that  the  calculated  air  temperature  at  the  end  of  14  days  in 
the  NBS  test  shelter  test  5  condition  is  in  better  agreement  with  the 
observed  value  when  solar  heat  effect  was  not  Included  for  the  ground 
surface  heat  exchange  than  when  It  was  Included.  The  ground  surface  was 
quite  wet  during  this  test  period  because  of  a  anow  prior  to  the  test. 
This  Implies  that  during  this  particular  test  period  the  solar  heat  was 
mostly  absorbed  by  evaporation  of  watar  on  or  naar  the  earth's  surfaca 
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am!  not  <  anse  *•  fable  temperature  rise  in  the  ground .  blmi  lar 

moisture  ev.i porat  ion  near  the  earth  s  surface  also  explains  the  small 
earth  temperature  rise  near  the  surface  for  other  test  conditions,  such 
as  the  Summerlin,  Broyles,  and  Napier  shelters,  some  of  which  will  be 
illustrated  later. 

Figure  12  compares  c-: leu  1. ted  program)  and  observed  air 

temperatures  in  the  Summerlin  shelter  for  its  summer  condition  test. 
Although  the  over.' 11  trend  of  the  actual  temperature  is  closely  followed 
by  the  calculation,  the  computed  amplitude  of  the  diurnal  temper'ture 
variation  within  the  shelter  is  considerably  smaller  than  the  observed. 

The  relative  humidity  calculation  for  the  summer  test  condition  of  the 
Summerlin  shelter  is  shown  in  figure  13.  For  the  first  five  days,  the 
calculated  relative  humidity  was  higher  than  the  observed.  This  dis¬ 
crepancy  may  be  due  to  the  fact  that  there  were  cool  regions  in  the 
shelter  interior  surfaces  where  more  condensation  of  water  vapor  was 
taking  place  than  was  determined  by  the  calculation  which  was  based  urv>r 
average  surface  tempera tures  of  each  exposure.  Some  of  the  high  peaks  in 
the  computed  relative  humidity  during  the  last  three  days  were  not  regis¬ 
tered  in  the  observed  record.  Figure  14  compares  the  calculated  earth 
temperature  surrounding  the  Summerlin  shelter  during  the  summer  test 
period  with  the  observed.  The  agreement  between  calculated  and  observed 
values  was  quite  good  for  the  south  and  east  walls.  The  earth  temperature 
under  the  floor  was  not  observed  during  the  test  because  this  was  a  pri¬ 
vately  owned  shelter  and  it  wae  decided  that  the  water-tightness  oi  the 
floor  should  not  be  Jeopard ited  by  making  a  hole  through  it.  The  cal¬ 
culated  roof  region  temperature  was  much  higher  than  the  observed  vs  lues. 
As  Mentioned  before,  this  discrepancy  ia  assured  to  be  ;aused  by  partial 
utilisation  of  solar  energy  by  surface  evaporation  of  ground  Moisture, 
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Fig.  13.  Comparison  of  the  calculated  and  observed  shelter  air  relative  humidities 
for  suttner  condition  test  for  Sucmerlln  Shelter  (computer  program  M- (2)  ). 
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resulting  in  less  temperature  rise  in  the  interior  of  the  earth  than  cal' 
culated.  It  was  found  for  all  other  shelters  for  which  the  ground  surface 
was  not  dry  and  bare,  the  inclusion  of  solar  heat  in  the  calculation 
usually  gave  higher  shelter  temperatures,  as  well  as  higher  roof  region 
temperatures,  than  observed  values. 

Figure  15  compares  the  computed  Sunsaerlin  shelter  thermal  environ¬ 
ment  with  the  observed  condition  for  a  moderate  climatic  condition. 
Although  the  computed  temperature  amplitude  was  again  lower  than  the 
observed,  the  agreement  between  the  computed  and  the  observed  data  is 
very  good  and  was  better  than  for  the  hot  summer  test  condition.  Three 
high  peaks  of  the  observed  shelter  relative  humidity  during  the  last  two 
days  of  the  test  were  not  reproduced  by  the  calculation  and  they  may 
represent  instrumentation  errors ,  Essentially  identical  results  were 
obtained  when  the  calculations  were  repeated  by  M-(3)  or  one-dimensional 
compound  model  on  Summerlin  shelter. 

Figures  IS  and  17  compare  calculated  results  with  observed  tem¬ 
perature  and  relative  humidities  of  the  Broyles  shelter.  As  described  in 
Section  4-C,  this  shelter  was  conditioned  by  a  cooling  coil  using  well 
water  during  the  first  two  days.  In  figure  16,  a  dashed  curve  shows  the 
computed  result  without  consideration  of  the  cooling  coil,  thus  yielding 
much  higher  shelter  temperatures  during  the  first  two  days  of  the  test. 

The  last  2-day  portion  of  figure  16  was  not  simulated  by  the  computer, 
so  the  comparison  between  the  calculated  and  the  observed  shelter  tem¬ 
perature  without  air  conditioning  should  be  made  between  August  i  and  16. 
However,  an  adjustment  was  made  later  to  account  for  the  observed  cooling 
capacity  of  the  coil  by  making  *n<*  negative  the  equation  for 
miscellaneous  heat  load  in  section  2.  The  first  2-day  portion  of  the 
Broyles  shelter  calculation  was  repeated  with  this  adjustment  and  the 
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results  are  shown  by  solid  dots  on  figure  16.  The  adjusted  and  cal¬ 
culated  shelter  air  temperatures  for  this  Initial  2 -day  period  agree 
very  well  with  observed  temperatures  obtained  at  the  geometric  center  of 
the  shelter 

Figure  17  conpares  the  calculated  and  observed  relative  humidi¬ 
ties  in  the  Broyles  shelter.  The  large  disagreement  of  the  computed 
relative  humidities  during  August  2  to  7  cannot  be  explained  adequately 
by  the  available  data.  The  agreement  of  the  computed  relative  humidities 
with  the  observed  values  for  this  test  were  generally  very  poor,  and  the 
same  trend  was  also  observed  In  the  Napier  shelter  comparison,  as  shown 
In  figure  19. 

The  Napier  shelter  temperature  comparison  was  satisfactory  if  the 
solar  heat  Input  to  the  ground  surface  Is  ignored,  as  seen  from  figure  18. 
The  inclusion  of  the  solar  heat  in  the  calculation  caused  the  computed 
shelter  air  temperature  to  be  approximately  3  to  4  degrees  higher  than  the 
observed  temperature. 

By  contrast,  figure  20  shews  excellent  agreement  between  the 
calculated  and  observed  temperature  and  relative  humidity  during  the 
winter  test  condition  for  the  Reading  shelter.  For  this  shelter,  nearly 
80  percent  of  the  total  heat  generated  was  conducted  into  the  surrounding 
earth.  The  moisture  balance  of  the  test  results  Is  shown  In  figure  21, 
Indicating  the  continuous  condensation  of  water  vapor  on  the  Inner  sur¬ 
faces  of  the  shelter  during  the  first  tan  days  of  the  test.  The  venti¬ 
lation  rata  was  either  3  cfm  or  lfc  efn  per  person  during  this  period. 

The  excellent  agreement  between  the  computed  thermal  environment  with 
the  observed  condition  during  the  first  10-day  period  la  a  good  demon¬ 
stration  that  the  computer  simulation  technique  was  valid.  However,  In 


MO 


figure  20,  the  computed  relative  huaidity  it  approximately  10  percent 

■\ 

Inter  dies  the  observed  during  the  lest  three  days.  Soring  this  perlcn, 
the  simulated  occupants  were  ioereaaed  iron  the  equivalent  of  50  to  100. 

In  figure  21,  the  indies ted  moisture  balance  within  the  Reading  shelter 
during  this  period  shows  that  the  ventilation  sir  actually  carried  out 
•or  e  water  vapor  then  was  produced  by  all  of  the  simulated  occupants* 
this  proves  that  the  shelter  wall  was  drying  out  during  the  period.  As 
LoLol*,  ihc  -uonputer  progzm  does  not  have  proviso*  for  si em¬ 
itting  the  drying  process  of  concrete  valla;  therefore,  the  computed 
relative  humidity  was  lower  than  the  obaerved. 

Figures  22  and  23  compere  the  eooputed  and  obaerved  earth  tempers- 
Cures  surrounding  the  Reading  shelter.  Being  in  the  midst  of  winter, 
the  undisturbed  earth  temperature  was  relatively  low  and  the  ground  sur¬ 
face  was  snow  covered  during  this  period.  Two  computer  models  M- (2) 
(symmetrical  end  three-dimensional)  and  M*<5  (one -dimensional  and  compound), 
ware  applied  to  the  Reading  shelter  calculation,  resulting  In  practically 
identical  thermal  environments.  The  ground  temperature  profiles  ware 
computed  end  compared  with  the  obaerved  profiles  for  other  prototype 
shelters  and  the  agreement  between  the  calculated  and  obaerved  results 
were  generally  similar  to  those  shown  in  figures  22  sad  23  for  the 
leading  shelter. 

The  computer  Model  M*<4),  basically  a  one -dimensional  end  easy  awed 
system  with  the  top  surface  of  tho  roof  being  adiabatic,  was  employed  to 
compute  tho  thermal  environments  of  the  ft.  falvior  200-man  end  1000 -man 
basaawnt-type  shelters.  These  results  axe  shown  graphically  in  Figures  24  . 
and  25.  Instead  of  comparing  tho  relative  hundlty,  as  In  eases  of  other 
shelters,  effective  temperatures  man  competed,  In  addition  to  tho  dry-bulb 
tampers turn  comparison,  tho  a amp at id  and  observed  results  agreed  almost 
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Fig.  23  Ca^trUat  of  tb*  caleulotod  and  ib«*rv«<3  oarth  tooporatvra  out  a  id* 
loading  •boitor  coiling  oed  floor  ( ccopwcar  prograo  N-(2)). 


perfectly  for  the  1000-aan  shelter.  Higher  end  less  flttdSuafciKg  dry- 
bulb  end  effective  temperature#  vers  calculated  for  the  200-«ee  shelter 
than  were  observed,  particularly  during  the  first  seven  day s  of  She  test 
period.  By  comparison  it  appeared  that  the  observed  temperatures 
reflected  shelter  diurnal  cycles  corresponding  to  a  higher  ventila¬ 
tion  air  rate  than  the  indicated  3  cfm  Par  person.  It  also  should  be 
realized  that  the  geometrical  relationship  between  the  inlet  for  the 
ventilation  air  and  the  stations  for  room  temperature  measurement  and 
the  mixing  effects  in  the  intervening  space  would  have  an  effect  in  the 
amplitude  of  the  temperature  variations  observed  in  any  of  the  shelter**. 

6.  CONCLUSIONS 

Except  for  a  very  few  cases,  such  as  for  the  Napier  shelter  and 
the  NBS  shelter  test  4  conditions,  the  computer  simulation  based  upon  a 
finite  difference  solution  of  the  heat  conduction  equation  for  the  cal¬ 
culation  of  the  thermal  environment  of  underground  protective  shelters 
has  bean  found  generally  satisfactory.  Good  agreement  between  calculat¬ 
ed  and  observed  thermal  conditions  was  obtained  for  tests  2  and  3  ih  the 
NBS  shelter,  the  Summerlin  shelter,  the  Reading  shelter,  and  the  Ft. 
Belvoir  1000-man  shelter.  Ihis  agreement  coincided  with  well-controlled 
test  conditions,  and  where  Input  parameters  were  Icnovr.  more  accurately 
than  in  other  cases.  However,  for  this  type  sf  analysis  there  are  two 
inherent  uncertainties  which  Influence  the  accuracy  of  the  calculations. 

The  first  uncertainty  stems  from  simplification  of  the  actual 
environmental  system  when  constructing  the  mathematical  model.  The 


complexity  of  the  system  including  construction  details  of  interior 
partitions,  foundations,  wall  supports  and  reinforcing)  the  heter¬ 
ogeneous  temperatures,  thermal  properties,  and  types  of  the  surround¬ 
ing  ssrthj  and  space  variation  of  conditioning  air  and  its  convective 
pattern  la  the  shelter,  makes  an  accurate  mathematical  description 
of  the  system  very  difficult.  Therefore,  the  simulated  computer 
model  is-  e  simplified  version  of  the  actual  system  that  Is  suited  to  the 
mathematical  handling  of  the  problem.  Tfea  magnitude  of  the  error  due  to 
this  uncertainty  has  not  been  evaluated,  but  it  is  a  part  sf  the  differ¬ 
ences  between  calculated  and  observed  temperature  and  humidity  conditions 
shown  in  this  report  for  seven  different  shelters. 

The  second  uncertainty  is  closely  related  to  the  first,  and  con- 

i 

cems  the  reliability  of  input  data.  It  is  almost  Impossible  to  acquire 
complete  three-dimensional  information  on  earth  thermal  properties  and 
temperature,  and  if  such  information  were  available,  it  would  be  difficult 
to  use  as  computer  imput  data.  Therefore,  when  a  homogeneous  heat  con¬ 
duction  model  is  employed  for  a  system  surrounded  by  earth  of  hetero¬ 
geneous  characteristics  and  temperature,  the  choice  of  proper  input  values 
becomes  very  difficult.  The  uncertainty  involved  in  evaluating  the  Interior 
surface  heat  transfer  coefficients  is  equally  aa  difficult,  as  discussed 
in  section  5. 

Thus  the  simplification  of  the  physical  characteristics  of  the  real 
shelter-earth  system  employed  in  creating  an  analytical  model  and  the 
uncertainties  In  the  thermal  properties  of  the  materials  involved, 
became  the  principal  approximations  incorporated  into  the  computer 


-  68  - 


analysis  of  the  thermal  environment  in  shelters.  In  spite  of  these 
uncertainties,  good  agreement  between  computed  and  observed  results  was 
achieved  for  many  of  the  shelters  desclrbed  in  this  report. 

In  addition  to  these  general  comments,  these  specific  conclusions 
are  emphasized,  as  follows: 

(a)  Soil  analysis  of  most  of  the  prototype  shelters  indicat¬ 
ed  that  thermal  diffusivity  and  thermal  conductivity  are  in  the  neighbor¬ 
hood  of  0.02  ft8 /hr  and  0.75  Btu/hr  fta(deg  F/ft),  respectively.  These  values, 
in  turn,  seem  to  result  in  fairly  good  agreement  of  earth  temperature 

change  during  the  test  period  for  most  of  the  shelters. 

(b)  Heat  conductance  at  the  shelter  inner  surfaces  for  most 

) 

of  the  summer  shelter  conditions  may  be  approximated  by  the  following: 

1.0  Btu/hr  ft?  deg  F  for  vertical  walls. 

1.5  Btu/hr  ft?  deg  F  for  ttye  ceiling. 

0.5  Btu/hr  ft?  deg  F  for  the  floor. 

.Although  there  may  be  some  other  values  and  other  combinations  of  these 
values  that  may  have  resulted  in  a  slightly  better  simulation  than  those 
used  in  the  analysis,  these  three  values  can  be  considered  representative 
design  heat  transfer  coefficients  in  the  underground  cavities. 

(c)  For  large  shelters,  the  one-dimensional,  compound  6- 
directicmal  model  M-0)  will  probably  be  adequate  for  calculating  the 
shelter  heat  transfer.  Thus,  the  complicated  5-diaensional  model  may  not  ba 
required,  at  least  for  tha  heat  transfer  calculation  of  leas  than  14-day 

* 

occupancy.  For  small  shelters  (such  as  a  family  shelter  similar  to  the  NBS 

shelter),  however,  it  is  recommended  that  tha  three- dimensional  modal  be 

used  for  the  accurate  calculations. 
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7.  APPENDIX 


Program  M-(4)  employs  a  subroutine  for  air  conditioning  the 
shelter.  This  subroutine  enables  the  computer  program  to  include  the 
heat  and  vapor  absorbing  capacity  of  a  given  cooling  systaa  (a  fan-and- 
coll  unit)  in  the  shelter  heat  balance  equation  of  2.1.  The  row  by  row 
account  of  the  cooling  coll  performance  is  considered  for  the  cross¬ 
counter  flow  circuited  coll  circulating  the  well  water.  The  Input  data 
required  for  this  subroutine  are: 

Shelter  air  dry-bulb  temperature,  t^  °? 

Shelter  air  wet-bulb  temperature,  t>t  °F 

Ventilation  air  dry-bulb  temperature,  t^  °P 

Ventilation  air  wet-bulb  temperature,  Ty,  °P 

inlet  coolant  temperature,  t.rf  F 

Recirculation  air  rate,  CFM^ 

Contact  factor  of  air  conditioning  coil  per  row, 

Thermal  resistance  between  the  solld-alr  inter¬ 


face  and  the  coolant,  6^ 

Estimated  temperature  rise  due  to  fan  heat, 
Coolant  heat  content,  C 


*P,  ft3  hr/Btu 


Rtu/hr,  °P 


The  contact  factor  per  row  ,  mentioned  above,  1«  a  function  of  air 
face  velocity  across  the  given  air  cooling  coll,  ait-slde  heat  transfer 
coefficient  of  the  coll  surface,  end  the  amount  of  total  heat  transfer 
surf act.  The  factor  may  be  estimated  by  the  following  expressions: 
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where f  *  air  aide  heat  transfer  coefficient,  Btu/hr  ft3  °P. 

a 

8^  «  air  side  heat  transfer  surface  of  the  coil  per  row,  ft3 

CFMt  *  CFK^  +  cyM  *  total  air  flow  rate. 

The  thermal  resistance  value  R  may  be  estimated  by  the  following  expression: 

w 


S  , 

_E  ~L 
s„  f 

t  w 
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where  St  *  coolant  side  heat  transfer  surface  of  the  coil  per  row,  ft3 
f  *  heat  transfer  coefficient  of  coolant,  Btu/hr,  ft3  *F.  t 
Dr  ■  sum  of  all  other  heat  resistance  between  the  coolant  side 
surface  and  air  side  surface,  such  as  thermal  resistance 
due  to  tube  wall,  finish  bond,  finish  metal  (depending  upon 
the  effectiveness  of  finish),  and  condensate  film  thickness. 

The  temperature  rise  due  to  fan  heat  can  be  estimated  by  total 
power  input  to  the  fan  and  air  flow  rate,  assuming  that  all  the  fan  heat 
will  be  expended  to  raise  the  air  stream  temperature.  The  value  of  &tf 
usually  does  not  exceed  3  *F. 

The  coolant  heat  Gy  is  the  coolant  mass  flow  rate  multiplied  by 
the  coolant  specific  heat.  In  the  case  of  a  direct  expansion  refrigerant 
coll,  Where  the  coolant  temperature  change  In  the  coll  is  very  small,  G* 
la  considered  infinity,  or  a  very  large  number.  The  details  of  the  cal- 
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culative  procedure  on  Cf  and  Ry»  end  their  design  value  can  he  found  In 
references  5  and  6. 

First,  the  air  inlet  condition  to  the  air  cooling  coil  is  calculated 
by  the  following  equation: 

(CFM)R(ta)  +  (CFMX^) 

**  “  CFMr  +  CFM  7-3 

(CFH)R(Wa)  +  (CFM)(Wy) 

1  (CFM)r  +  CFM 

where  t£,  W^,  Wa,  and  Wv  represent  inlet  air  dry-bulb  temperatures,  inlet 
air  humidity  ratio,  shelter  air  humidity  ratio,  and  ventilation  air 
humidity  ratio,  respectively.  These  humidity  ratio  values  are  calculated 

by  a  separate  subroutine  from  dry-  and  wet-bulb  temperatures,  as 
explained  in  section  2.7. 

The  air  conditioning  coil  is  assumed  to  be  of  multi-row  structure, 
and  the  overall  direction  of  coolant  is  counter  to  that  of  air  flow, 
although  for  individual  rows  t'le  coolant  is  flowing  perpendicular  to  the 
air  stream.  For  simplicity,  it  is  assumed  that  the  coolant  temperature, 
as  well  as  air-side  surface  temperature  of  a  row,  changes  by  a  step  func¬ 
tion.  The  heat  and  vapor  transfer  calculation  of  the  counter-flow  dehuaidl- 
fying  coil  requires  an  iterative  procedure,  because,  except  for  the  direct 
expansion  coll,  ths  outlet  condition  of  air  or  coolant  is  not  prevlouely 
known.  In  this  report,  the  outlet  coolant  temperature,  which  la  in  the 
same  aids  as  the  inlet  air  condition  of  the  coll,  it  first  approximated . 

The  heat  and  vapor  transfer,  end  subsequent  reduction  of  air  tempsrsture , 
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air  hualdity,  and  change  of  the  coolant  fluid,  are  then  calculated 
row  by  row.  After  the  calculation  la  completed,  the  total  net  change 
In  coolant  temperature  la  aubtracted  from  the  outlet  temperature  of 
the  coolant,  yielding  the  calculated  coolant  temperature  at  the  Inlet 
condition.  If  the  calculated  coolant  temperature  at  the  coll  coolant 
Inlet  is  different  from  the  actual,  the  calculation  Is  repeated  with  a 
modified  outlet  coolant  temperature  until  the  calculated  agrees  with  the 
given  Inlet  temperature  of  the  coolant  flow.  The  general  heat  and  vapor 
transfer  relation  used  for  this  calculation  is  described  for  1  row  as 
follows : 

Cf  [l.08  CFM(tt-tsl)  +  4.5  (CVM)aO<tf1-ff,1)]  »  3%.!-^)  7-5 


US£  in  the  above  equation  is  the  humidity  ratio  of  the  air  saturated  at 
the  surface  temperature,  if  *  W.i’  the  second  term  in  the  left 

hand  side  of  the  equation  Is  set  equal  to  zero.  Since  W8i  is  a  compli¬ 
cated  function  of  t8l,  an  iterative  technique  l#  required  to  aoive  tgj 
from  the  above  expression.  After  t,^  is  obtained,  the  leaving  air  condi¬ 
tion,  ti+i,  from  the  l  row,  and  entering  coolant  condition  to  the 

1+1,  are  calculated  by  the  following  relations; 

*i  *  U+l  *  <*1  “  *s)(Cf)  7-6 


Wt  -  Vui  -  0tt  -  W.i)(Cf)  If  Wt  >  Wgt 
•  o  if  \it  *  »tl 


7-7 


— ^(tgi  -  twt)  *  (t*i  -  tvi+l)(Gw) 


7-8 


Since  all  of  the  relation*  are  linear,  calculation*  of  vl+l*  **** 

ty'i+i  are  straightforward. 

When  the  total  nunber  of  row*  1*  S  and  the  proper tie*  In  the  air 
inlet  point  are  specified  by  subscript  1,  the  properties  of  the  air  out¬ 
let  condition  should  be  specified  by  N  +  1.  Hum,  the  calculations  will 
be  iterated,  as  mentioned  before,  until  t*tn+i  becomes  equal  to  t*.  The 
final  reaulta  of  the  air  conditioning  capacity  will  be  expressed  es 
follows : 

(1)  Sensible  cooling  capacity,  Q^s  ■  (l.OSXCIMXt^-tg+i+Afcjf) 

(2)  latent  cooling  cepeclty,  Q^l  *  (4.5)(CW)(Wi-4fj|+i)fc,. 

With  the  value  of  Qgg  end  Q^l  known,  the  shelter  air  condition 

for  the  next  time  period  now  can  be  calculated  by  adding  to 

the  overall  heat  balance  equation  in  (1). 
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Computer  Programs 

Fortran  listings  of  the  computer  progress  developed  for  the 
heat  transfer  studies  in  underground  protective  shelters  during  the 
contract  research  are  attached  in  the  following  pages.  Except  for 
M-3  program,  the  computer  symbols  are  explained  according  to  input 
and  output  sequences  at  the  beginning  of  each  program. 

A  complete  set  of  input  data  and  illustrative  examples  of 


output  format  for  program  M-4  are  also  shown  . 


KPl:  A  matrix  index  marking  the  interface  cf  earth  and  north  wall. 

Kr’.:  A  matrix  index  marking  the  interface  of  earth  and  south  wall 

XF3:  A  matrix  index  marking  the  end  of  earth  block  outside  the  south 
wall. 

KQl:  A  matrix  index  marking  the  interface  of  earth  and  east  wall. 

KQ2:  A  matrix  index  marking  the  interface  of  earth  and  west  wall. 

F.Q3:  A  matrix  index  marking  the  end  of  earth  block  outside  the  west 
wall. 

KR1:  A  matrix  index  marking  the  interface  of  earth  and  shelter  roc:5. 

KR2:  A  matrix  index  marking  the  interface  of  earth  and  shelter  floor. 

KR3:  A  matrix  index  marking  the  end  of  earth  block  below  the  shelter 
floor. 

KP:  Nearest  integer  to  (KPl  +  KP2)/2. 

KQ:  Nearest  integer  to  (KQl  +  KQ2)/2. 

KR:  Nearest  Integer  to  (KR1  +  KR2)/2. 

NN:  Number  of  time  Iol  ments. 

LC:  Number  of  matrix  points  for  shelter  inner-wall. 

WC:  Work  cell  for  determining  type  of  initial  earth  temperature  data. 

WC  »  -1.  :  Depthwise  variation 

•>  o.  ;  Homogeneous  earth  temperature 

-1.  :  Six  directional  variations. 

TGI:  A  constant  Initial  earth  temperature  whan  WC  ■  0. 


TCi(I),  TG2(I) ,  TG3(I),  15G4(X),  1585(1),  TC6(I):  Initial  earth  tamper* 
tare  profiles  in  earth  blocks  outside  of  north,  south,  east, 
west,  roof  and  floor,  respectively.  *F 
TC(L,M)s  Initial  temperature  distribution  of  the  shelter  inner  wall, 
®F. 

2 

S(M):  Shelter  inner  wall  surface  area,  ft  . 

D(M):  Thickness  of  shelter  inner  wall,  ft. 

ZL(M):  Thickness  of  the  earth  block,  ft. 

CK(K):  Thermal  conductivity  of  shelter  wall,  Btu/hr,ft, *F 

CG(M):  Thermal  conductivity  of  earth  block,  Btu/hr, ft, *F 

2 

AC(M)j  Thermal  diffusivity  of  shelter  well,  ft  /hr. 

2 

AS(M);  Thermal  diffusivity  of  earth  block,  ft  /hr. 

H(M):  Heat  transfer  coefficient  at  the  shelter  inner  surface, 

2 

Btu/hr,ft  ,°F. 

HD(M):  Mass  transfer  coefficient  at  the  shelter  inner  surface, 

lb/hr,  ft2 (lb/lb). 

E3(H):  Heat  transfer  coefficient  at  the  external  side  of  the  earth 

2 

block,  Btu/hr,  ft  ,  *F. 

DB(N):  Outdoor  air  dry-bulb  temperature,  *F, 

TV(CN)t  Ventilation  air  dry-bulb  temperature,  *f. 

DPV(N):  Ventilation  air  dew-point  temperature,  *F. 

2 

Q(N):  Solar  radiation-*. Btu/hr, ft  . 

A:  Shelter  length  (internal  dimension),  ft. 

Bt  Shelter  width  (internal  dimension),  ft. 

C:  Shelter  height  (internal  dimension),  ft. 

C-3 


QVL:  Latent  Heat  Exchanged  with  Ventilation  Air,  Btu/hr. 

QGS:  Sensible  Heat  Generated  within  the  shelter,  Btu/hr. 

QGL:  Latent  Heat  Generated  within  the  shelter,  Btu/hr 

QWST:  Heat  Transferred  to  the  shelter  inner  surfaces,  Btu/hr. 

QWLT:  Latent  Heat  Transferred  to  the  shelter  inner  surfaces,  Btu/hr. 

TWCT:  Water  Vapor  Collected  on  the  shelter,  inner  surfaces,  lb/hr. 

TQVS ,  TQVL,  TQGS,  TQGL,  TQWST,  TQWCT,  TTWCT:  Cumulative  values  for 
QVS ,  QVL,  QGS,  QGL,  QWST,  QWLT  and  TWCT,  respectively. 

AT(N)  =  TA  Calculated  shelter  air  dry -bulb  temperature,  °F. 

RH(N)  =  RHA:  Calculated  shelter  air  relative  humidity, 

DP(N)  =  DPA:  Calculated  shelter  air  dew-pcint  temperature,  °F. 

Subscript  M  refers  to  exposures  such  that 
M  =  1:  North  wall  and  its  region 
2:  South  wall  and  its  region 
3;  East  wall  and  its  region 
4:  West  wall  and  its  region 
5:  Roof  and  its  region 
6:  Floor  andjgits  region. 

Subroutine : 

GN :  Computes  the  moisture  saturated  air  vapor  pressure  for  a  given 
temperature  using  Goff  and  Gratch  formula. 

FN :  Computes  the  boundary  temperature  by  a  finite  difference  formula. 

QC :  Computes  the  human  metabolic  heat  by  type  of  occupants  as  function 


C  'UNDERGROUND  FALLOUT  SHELTER  THERMAL  ENVIRONMENT  HEAT  TRANSFER  1 185000C 

C  NATIONAL  BUREAU  OF  STANDARDS  PROJECT-10436  T.KUSUDA,  10.3  11850011' 

DIMENSION  TGI 23. 22. 29) , TCI  10. 6 ) .QSUNl 170, 6) ,DB 1 170) , TV! 170) 

OPVI 170),TG1( 10) , TG2  C 10 ) ,TG31 10) , TG4 1 10) ,TG5( 10)  11850031 

TG6I 10 ) , TGO 1 30 ) , S 1 6 ) ,D ( 6) ,ZL (6) . AG (6 ) . AC 1 6 ) ,CK (6 ) ,CG( 6 ) 1 185004C 
H(6),HD(6),HG(6),QWS(6) »QWL ( 6) , TOWS! 6 ) »TQWL 1 6) ,TWC (6 )  1185005C 


TTWCl 6 ) .QSUNT 16), 
DDI6),EI6), 0(170) 


ATI  170), RHI 170) .DPI  170) , QSUN0I6) 


3 

1 

99 

12 

2 

11 

10 


DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
DIMENSION 
FORMAT! 8F9.4) 

FORMAT ( 1414, 2F3.0) 

READ  1,KP1,KP2?KP3,KQ1,KQ2,KQ3,KR1,KR2,KR3,KP,K0,KR,NN,LC,V;C 
IFIWC)  2,  12,11 
READ  3, TGI 
GO  TO  10 
GO  TO  13 


READ 

READ 

READ 

REAO 

READ 

READ 

REAO 

00 


3, 

I  TGI I  I ) ,  I* 

1 ,  KP 

1),  l 

3, 

ITG4I J ) , J» 

KQ2 , 

KC3) 

3, 

( (TCIL.M), 

L=  1 , 

LC), 

3, 

I  SIM), M«l, 

6), I 

DIM) 

3, 

( CGI M ) , M* l 

,  6 )  , 

I  AC  I 

3, 

(HD(M),M*1 

,6 ) , 

IHGI 

3, 

I  OB  I N ) , N= 1 

,NN ) 

,  I  TV 

55 

M=  l ,  6 

TG2I  I )  ,I=KP2,KP3) ,  (TG3I.J)  ,J=1,K01) 

,  (TG5IK)  ,K*l,KRl) ,  ITG6U)  ,K=KR2,KR3) 

M= 1 , 6 ) 

,  M  =  1 , 6 ) , IZLIM) ,M=1 ,6  ) , ICK(M) ,M=1,6) 

M)  ,M*1,6) , (AGIM) ,M=1,6) , (HIM) ,M=1,6) 

M) ,M=1 ,6) 

(TVIN),N=1,NN) , IDPVIN) ,N=1,NN> , I Q I N ) , N=1 , NN) 1 1850401 

1 18504  1C 


1185006C, 

U85007C 
11851 
11851 
l 185010C 
118501 1( 
11850121 
1 185013C 
1 185014C 
1 185015C 
11851 
1 1 8501 7( 
l 1850181 
1185019C 


51 

52 

50 

c4 

5 

9 

14 


13 


15 

300 


IFO-G(M)  )  50, 51, 50 
00  52  N*1,NN 

CSUN(N,M)-0. 

GO  TO  55 
DO  54  N= 1 , NN 
GSUN l N , M ) =C l N ) 

CONTINUE _ 

READ  3,7a,B,C,ZN,1.ZN2,ZN3.TA.OPaI 

FORMAT! 4F9. 4)  & 

READ  l4,|77TrrZN,tlTTl 
GO  TO  300  ‘ 

READ  3,(TGC(K),K=1,KR3) 

GO  TO  10 
FORMAT  I  4  3HI 
PRINT  15 
PRINT  16 


PB , 8S ,BL ,DGX,DGY ,DGZ , DGX l , DGY1 


- 1  r—i 


Lf  Ell 


UNDERGROUND  FALLOUT  SHELTER  DATA) 


1185042C 

11850431 

1 185044C 

1185045C 

1 1  85046C 

11850471 

1 1850481 

11850491 

11851 

11851 

1185052C 

11850531 

1185054C, 

11850551 

1185056! 

11850571 


16 

FORMAT  I 79H0  ZN 

ABC 

V 

G 

1185058f 

1  ZN 1  ZN2  ZN3 

DT  ) 

1185059! 

30 

FORMAT(F9.0,F7.l,2F8 

.1,F10.0,F12.0,F7.0,3F6.2) 

PRINT  30.ZN, A, B,C,V, 

G,ZN1,ZN2,ZN3»DT 

1 185059! 

PRINT  17 

1185060! 

17 

FORMAT  I 64H0  M 

1  2  3 

4 

5 

1185061! 

1  6) 

1185062! 

PRINT  18, (CK(M),M=1, 

6  ) 

1185063! 

18 

FORMAT  I 10H0  CKIM) 

6F10.3) 

1185064! 

19 

FORMAT ( 10H0  CGI  M) 

6F10.3) 

1185065! 

20 

FORMAT! 10H0  AC! M) 

6F10.3! 

1185066' 

21 

FORMAT! 10H0  AG ( M) 

6F10.3) 

1185067' 

22 

FORMAT! 10H0  DIM) 

6F10.3) 

1185068' 

23 

FORMAT! 10H0  ZLIM) 

6F10. 3) 

1185069 

24 

FORMAT! 10HO  S!M) 

6F 10. 3  ) 

1185070! 

5 

FORMAT ( 10H0  H ! M ) 

6F10. 3) 

1 185071 

<.6 

FORMAT ( 10H0  HG1M) 

6F10.3) 

1185072 

27 

FORMAT! 10H0  HOIM) 

6F10.3) 

1185073 

PRINT19, ( CG ! M ) , M= 1 , 6 

) 

1185074 

PRINT20,  ( AC ( M ) , M= 1 , 6 

I 

1  1850/5 

PRINT21, (AG(M),M=l,fc 

) 

1185076 

C-6 


PRINT22, (DIM) ,  Mil, 6) 
PRINT23,(ZL(M),M*l,6) 

PRINT24, ( SIM) »  M*l,6)  * . . 

PRINT25* (HIM),  M*l,6) 

PRINT26, IHGIM)»M*1,6) 
PRINT27,IHD(M),M*1,6) 

GO  TO  100 

C  GROUND  TEMPERATURE  INITIALIZATION 

100  IFIWC)  101,102*103 

101  00  104  1*1, KP3 

00  104  J*l,KQ3 
00  104  K* 1 ,KR3 

104  TG(I,J,K)*TGO(K) 

GO  TO  129 

102  00  105  I-l,KP3 

00  105  J*1,KQ3 

00  105  K*1,KR3 

105  TGI  I , J,K)*TGI 
GO  TO  129 

103  00  900  1*1, KP3 

00  900  J*1,KQ3 

DO  900  K*1 ,KR3 

IFIKPl-I)  111,106,106 

106  IFtJ-KCl)  111,107,107 

107  IFIKQ2-J)  111,108,108 

108  IFIK-KR1)  111,109,109 

109  IFIKR2-K)  111,110,110 

110  TGI I , J,K )  *TG1{I) 

GO  TO  900 

IFII-KP2)  117,112,112 

112  IFIJ-KC1)  117,113,113 

113  IF1KQ2-J)  117,114,114 

114  IF(K-KRl)  117,115,115 

115  IFIKR2-K)  117,116,116 

116  TG( I , J ,K )  =  TG2< II 
GO  TO  900 

117  1 F ( KQ 1- J )  121,118,118 

118  IFIK-KR1)  121,119,119 

119  IFIKR2-K)  121,120,120 

120  TG( I,J,K)*TG3( J> 

GO  TO  900 

121  IFIJ-KC2)  125,122,122 

122  I F ( K— KR l )  125,123,123 

123  I F ( KR2-K )  125,124,124 

124  TGI I,J,K)*TG4(J) 

GO  TO  900 

125  IFIKRi-K)  127,126,126 

126  TGI  1, J,K)  =  TG5(K) 

GO  TO  900 

127  IFIK-KR2)  900,128,128 

128  TGI  I , J ,K )*TG6 1 K ) 

900  CONTINUE 

1900  FORMAT  I 10F9.5) 

129  TM=0.0 
TQWST=0. 

TQWLT=0. 

TTWCT::0. 

TQVS=0. 

T  Q  V  L  =  0  . 

T  CGS - 0 . 

•C-7 


11850771 
U850781 
11850791 
11850801 
11850811 
11850821 
11850831 
11850841 
11850851 
11850861 
11850871 
11850881 
11850391 
11850901 
11850911 
L185092I 
11850931 
11850941 
11850951 
11850961 
11850971 
11850981 
11850991 
11851001 
1185 10 1( 
11851021 
11851031 
11851041 
11851051 
11851061 
11851071 
11851081 
11851091 
11851101 
118511H 
11851121 
i  185 1 1 31 
11851141 
11851151 
11851161 
11851171 
11851181 
11851191 
11851201 
11851211 
11851221 
1185123: 
11851241 
1185125! 
1165126! 
1  185127' 
11851281 


1185133 

1185134 

1185135 

1185136 

1185137 

1185138 


1 


T0GL=0.  h 

11851391 

TTOSUN-O.  f 

•  ■; 

1 185140( 

00  152  M«l,6 

1185141C 

TQWS(M)=0. 

1185 142( 

TQWLIM)»0. 

■v 

•■v  ;  '  . 

1185 143( 

TTWC(M)»0. 

■  jv  .  ■■■ ' 

1185144C 

152 

QSUNT (M)*0» 

1185 145( 

LLC-LC-1 

:  V  ,  ;  "  3 

11851 

ZLC*LLC 

C  END 

OF  INITIALIZATION 

1185146C 

DO  800  N-1,NN 

1185147( 

29 

FORMAT  ( 10F5. 1 ) 

11851 

301 

00  28  M«l,6 

;*  •  '  ' 

28 

PRINT  29,  (TC(L,M),L»1,LC> . 

11851 

PRINT  29,  ITG( I ,KQ,KR ) , I«1 ,KPl ) 

11851 

PRINT  29,  (TGI  I ,KQ,KR) , I*KP2,KP3) 

11851 

PRINT  29,  (TG(KP,J,KR),J*1,KQ1) 

11851 

PRINT  29,  (TG(KP,J,KR) , J*K02,KQ3) 

11851 

PRINT  29,  (TG(KP,KQ,K),K*1,KR1) 

11851 

PRINT  29,  ITG(KP»KQ,K),K«KR2,KR3) 

11851 

1501  FORMAT ( 16F5.  1 J 

PRINT  1501,  TGI  I,  1, 1),TGU.1,KR3)  ,TG(1,KQ3»1)  ,TG( l ,KQ3,KR3)  ,TG(KP311851 
1, 1,1),TG(KP3, l,KR3),TG(KP3,KQ3,l),TG(KP3,KQ3,KR3)*TGIKPl,KQl,KRl),U85l 
lTGtKPl,KQl,KR2)»TG(KP 1 , KQ2,KR1 ) , TG( KP1 ,KQ2 ,KR2) , TGIKP2 ,KQ1 , KRl ) » TGI 1851 


l(KP2»KQI,KR2),TG{KP2,K02,KRl)»TG(KP2fKCI2,KR2)  11851 

gTM«TM+DTT  11851 

%SV«GN(OPV(N) )  1185148 

WV*0.622*PSV/{ PB— PSV)  1185149 

OTM»0.  11851 

1.j5  DTM-DTM+DT  11851 

DTA»0.2 
TA1*40.0 

IR=1  11851 

405  QVS=*  l »08*G*  ( TA  l— TV  IN)) 

QGS*QG (TA1 ,ZN1,ZN2,ZN3,1.} ♦BS 
QWST-O. 

DO  153  M=*  1 » 6 

DO (M)*D(M}/ZLC  11851 

QWS(M)*H(M)*S(M)*ITA1— TCI l»M) ) 

153  QHST*QHST  +QWS (M) 

GO  TO  (1153,406), IR  11851 

1153  ZX=QGS-QVS-QWST  11851 


IF(ZX)  400,401,402 
402  TAI=TA1+DTA 
ZX2*ABSF(ZX) 

GO  TO  405 
400  T A2=  T A 1-OT A 
ZX  1  =  ABSF ( ZX ) 

TA=ITA1*ZX2+TA2*ZX1)/(ZX1+ZX2> 

I  R  =  2  11851 


FA1=TA 
GO  TO  405 
401  T  A=  TA 1 
406  DDPA=0-2 
OPA 1=40.0 
I  R=  1 

409  PSA^GN ( OP A  1 ) 

WA=0.622*PSA/ ( P8-PSA ) 
QVL=4780»*G*( WA— WV ) 
GML T=0. 

T  W  C  T  ~  0 . 

CO  153  M=  1 , 6 


11851 


11851 

) 185155 
1 18516C 
1185I6C 
118! 

1  1 8  ( 
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157  PSN-GNITCll.M))  116516; 

KS-0.622-PSM/IP8-PSNJ  118916; 

IFIMA-WS)  154,154.155  119*16; 

154  QMLIMJ-O.  118518; 

CO  TO  156  1185164 

155  ONLIN)»1061.-HOIM)«SIM)*INA-NS)  1165164 

•56  QMLT«QWLT*QNLIN)  118516* 

TMCIN)-QNLIH)/I06l.  1169165 

158  TMCT«TMCT*TMCIN)  1185165 

QGL«QGITA,ZN2,ZM2,ZN3,0.)m  1I651TC 

GO  70  ( 1158,413), IK  11891 

1158  ZY-QGL-OVL-QMLT  11851 


IFIZY)  410,411,412 
412  DPAI-0PAI*C0PA 
ZY2-A6SFJZV) 

GO  TO  409 

410  0PA2-0PAI-C0PA 
ZYi-ABSFIZV) 

OPA-  I  OP A 1 *2 Y2+DPA2-Z Y 1 1 / ( Z Y  UZ Y2 } 


16-2  11851 

0PA1-0PA 

GO  TO  409  11851 

411  0PA-0PA1 
413  POA-GNI TA ! 

RHA- 100* -PSA/POA  118517$ 

00  161  6-1,6  1165176 

TQMSm-TQHS(M)*QMS(M}-DT  1185176 

TQMLIN)-TQWLIM).ONLIN)-DT  1169177 

ZM-CDIM)-tQMS(M)*QML(M) )/(CK(M)-S(MI)*TC(2,Ml  1185165 

160  rC(l,M)-ZM  1169166 

161  TTWC<M)aTTMC(M)*TWC(N)«OT  1165177 

TQMST-TQMST,OMST*OT  1165173 

1QNLT-TQMLT*0MLT*0T  1165178 

TTMCT-TTMCT-TNCT-OT  1185179 

TQGS-TQGS+GGS-OT  1185179 

TOGL-TQGL+QGL *0T  1169160 

TQVS-TQVS-OVS-OT  1189160 

TQVL-TQVOQVL-DT  1169161 

C  CONCRETE  NALL  TEHPERATl*R£  DISTRIBUTION  1165602 

00  164  MM, 6  1169603 

6IN|-ACIM)*0T/00IN>—2  1169204 

00  164  L-2.LLC  11691 

164  rClL,N»»EIN|-|TCIL-l.NIMCIL-'  M)*I l./EIM)-?. 1-TCIL.N1 I  1165206 

I* IOTN-DTT )  1405,1164,1164  11651 

C  GROUNO  AND  CONCRETE  80UN0ARY  TEMPERATURES  1165607 

1164  00  165  J-KC1.KQ2  11851  ’* 

00  165  K-RR1.KA2  1165209 


TCIKPI,J,*I  «PNICRm,CG(l},TCILLC,l),TG(KPl-l,J,K),  0011)  , OCX, 0. 11185 1 


165  T6IKP2, J,K  «  -FNICKI 2) ,CGt 2) ,TC ILLC ,2) ,  TGIKP2M  ,J«K) , 
00  165  I-RP1*KP2 
CO  166  K-RR1.XR2 

t3ll,KQl,K>  -PNI CRI3)«C&l9),TCtLLC»3t»TGII , KGl-l»K) , 
146  TGU,KC2,K>  -FNICXI4)  »CG14)  .TCILLC.4I  »TGI  l»K02M»IU  » 
00  167  I-KP1 ,RP2 
00  167  J*R01,R02 

TCI I,J,KR1)  -PNICKI 5), CGI 5) , TCILLC.5) , TGI t , J,RRl-l) , 
167  TGI  t« J.KR2)  -PNICXI 6) ,CGI6) , TCILLC.6) ,T6I (,J,«R2«1 1 , 
TCILC,  H-TGIRPl,KO,RR) 

T£lLC,2)-rCIKP2,R0,RR) 

TCILC, 31-TGIXP, R01, *R» 

TCILC, 4|-TGIKP,fcQ2,Kft1 
TCILC, 5I>TG(KP,K0,XR1) 

TCILC, 6I*TGIRP*K0,KR2)  _  . 

PfclRIT  1900,  I  TCILC, N),N*  1,6) 


001 2) .OCX, 0.111651  . 

1169612 

1169613 

com. OCY. 0.1 11651 
00141,067,0. 111691 

1169216 

1169217) 

001 9), 062,0. 111691 
0016), 062 ,0.1 11691 
11691 
11691 
11691  , 
11691  ! 
11691  !. 
11691 


C  GROUND  TEMPERATURES 
KP4  a  KP3-1 
KQ4  *  KQ3-1 
KR4-^  KR3-1 


DO  f.QO 
DO  600 
00  600 
IFIKPl-I ) 

168  IF(J-KCl) 

169  IF { KQ2-J ) 

170  IF(K-KRl) 

171  IFIKR2-KI 

172  AGG=AG(l) 
IF(I-KPl) 

173  I F ( I—KP2 ) 

174  IF(J-KCi) 

175  IFIKQ2-J) 

176  IF (K— KR1 1 

177  IFUR2-K) 

178  AGG=AG(2> 
IFIKP2-I ) 

179  IF(KQl-J) 
130  IF(K-KRl) 

181  I F { KR2— K ) 

182  AGG«AG( 3 ) 
IF { J— KC l ) 

590  IFII-KP1) 
580  IF  IKP2-I 
•’3  I F  (  J—KC2 ) 

- ->4  IF(K-KRJ) 

185  IFIKR2-K) 

186  AGG=AG ( 4 ) 
IF  (KQ2-J: 

570  IF  II-KP1 
560  I F ( KP2-  I  ) 

187  IF(KRl-K) 

188  AGG= AGt  5 ) 
IF  IK-KR1) 

550  IF  U-KP11 
540  IF  (J-KQ1) 
530  IF  (KQ2-JJ 
520  I F ( KP2-  I  ) 

189  IFIK-KR2) 

190  AGG= AG ( 6 ) 


1=2, KP4 
1=2,  KQ4 
:*2,KR4,, 

173.168.168 

173.169.169 

173.170.170 

173.171.171 

173. 172.172 

500,600,600 

179.174.174 

179.175.175 

179.176.176 

179.177.177 

179.178. 178 

500,600,600 

183.180.180 

183.181.181 
183, 182,18? 

500,590,590 
500,580,580 
500,600,600 
187,  184,184 

187.185.185 

187.186.186 

500,570,570 
500,560,560 
500,600,600 
189, 188, 163 

500,550,550 
500,540,540 
500,530,530 
500,520,520 
500,600,600 
600,190, 190 


I  ”•  {  KR2-K  )  500,521,521 

521  IF  (I-KP1)  500,522,522 

522  iFIJ-KCl)  500,523,523 

523  IFIKQ2-J)  500,524,524 

524  I  r ( KP2-  I  ’  500,600,600 
500  IF i  I-KPl)  502,502,503 

503  IFIKP2-I)  502,502,504 

504  OGX  =  DGX 1 

502  IF(J-KCl)  501,501,505 

505  IFIKQ2-J)  501,501,506 
6  DGY=DGY1 

5<)1  EGX  =  AGG*DT  r/ (  DGX**2  ) 
tGY  - AGG*OTT/ I DGY**2 ) 


EGZ=AGG*CTT/(CGZ**2) 
T 1  =  TG I  I- 1 , J , K  ) 

T  2  —  T  G  1  I  +  1 ,  J  ,  K  } 


1165220 

1185220 

1185220 

1185220 

1185221 

1185222 

1185223 

1185224 

1185225 

1185226 

1185227 

1185228 

1185229 

1185231 

1185232 

1185233 

1185234 

1185235 

1185236 

1185238 

1185239 

1185240 

1185241 


1185243 

1185244 

1185245 

1185246 


1185248 

1185249 


1185251 

1185252 


1185252 
1185252 
1 185252 
1185252 
1165252 
1 185252 
1 1851 
11851 
1  1851 
1185256 
1  18525  7 
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T3*TG! I,J-1,«)  1185 

T4-TGI 1185 
T  5  *  T  G  ( 1 » J  *  K—  1 )  1185 

T6=TG!  I,J,*Ul7  1185. 

XX  =  tGX*!Tl*T2)+£GY*!T3*T4-)  ♦£&£*  i  T  5'  t61  1185. 

YY» I »— 2 . • ! EGX+EGY+EGZ )  1185. 

TG( I ,J,K>*XX+TG( I »  J  »K ) *YY  1165. 

600  CONTINUE  1185. 

C  GROUND  SURFACE  TEMPERATURES  1185; 

DO  191  J* 1 » KQ3  1185; 

00  191  K-i,KR3  1185. 

TGI  1  *  J  »K  )  *FN  (  HG(  1  )  ,CG(1)  ,0b  ( N  )  ,  1  G(  2  ,  J  ,K  )  ,  1 .  ,UT.X  ,  OSUN  i  N,  11  )  1185; 

191  TG(KP3,  J,K)*FN<HG12  ),CG<  2>,D8lN),  TS(KP3-1»J»K)  ,  1.  ,DGX,CSUN{N,2  )  )  1185; 

CO  192  1*1, KP3  1185; 

CO  192  K*L,RK3  1185. 

TGI  I  ,  1,  K  IsFNt  HGI  3  )  ,CGl  3 ) »  DO! N  1  ,‘TGC  I  ,  2,  K) , l. *DGY  , uSUN ! N, 3  )  1  1185; 

192  TGI IfKC3,K)*FN(HG(4) »CG14 ) , Uti I N ) , TG( I , KQ2-1 ,K) , 1 . , DGY,USUN(N,A ) )  1185. 

DO  193  Ix  1  ,RP3  ‘  1185.‘ 

DO  193  J31 *<Q3  1185; 

TGII.J,  1  )*FN(HC(5),CGI5),Dtt(Nj#TGCI,  J,  21 , 1 . ,DG7 .CSUNIN-5 1)  1185; 

193  TGI l,J,KR3)*TGt I » J » KR 3 ) 

1185; 


lies; 


AT l N ) *  T  A 
RH I N ) =  R  HA 
DP  I N 1 =CPA 

WRITE  OUTPUT  TAPE  81 , 29, I  TCI  1 , M 1 , M* i , 6) , TGI  1 , l , 1 1 , TGI  1 , l ,KR3) 


WRITE  OUTPUT 

TAPE  81,1 

900, TM 

195 

FORMAT! 10H0 

UWSIM) 

6F15.0) 

196 

FORMAT  I 10H0 

towsim; 

6F15.01 

19  7 

FORMAT  I 10H0 

QWL(M) 

6F15.01 

198 

FORMAT  I 10HO 

TOWL(M) 

6F15.0I 

199 

FORMAT! 10H0 

TWC(M) 

6F15.0) 

200 

FORMAT! 10H0 

TTW! (Ml 

6F15.0) 

201 

FORMAT! 10H0 

TC< l.MJ 

6F15.1) 

202 

FORMAT! lOOHO 

l 

2 

3 

1  4 

5 

6 

J 

203 

FORMAT ( 100*1 

N  TM 

ovs 

QVL 

l 

QGL 

CWST 

OWL  T 

) 

204 

FORMAT ( 100H 

TWCT 

TQVS 

TOVL 

1  TC.GS 

TQGL 

TUWST 

) 

205 

FORMAT ( 100F 
1TTCSUN 

T  OWL  T 

TTWCT 

J 

TQSUN 

206 

FORMAT ( l I  3, 

1 F  6. 1 ,6F15 

.0) 

207 

FORMAT ! 6F  15 

.2  ) 

1185; 

nos; 

ms; 

U85; 

lies; 

H85; 


1135; 

1185; 

lies; 

U85; 

U85; 

lias; 

H85; 

lies; 

1135; 

0GS1185; 

H85; 

lies; 

ii85;- 

1185; 

11852 

lies;. 


FORMAT (5F  15.0) 

1185; 

PRINT 

203 

1135; 

PRINT 

206,  N , rM,OVS,QVL,OGS,OGL,OWSTfOwt T 

11852 

PRINT 

204 

11352 

PRINT 

207,  T WCT.TQVS.TQVL, TOGS, TOOL, TOwST 

11852 

PRINT 

205 

11 65  2 

PRINT 

208,  TQWLT, TTWCT, TT'OfrfcMt 

11852 

PR  I  NT 

202 

1185  2 

PRINT 

195, ( OWS! M) ,M* 1,6 ) 

11852 

PR  i  NT 

196, ( TQWS(M) ,M=1,6) 

11852 

p  ■>. !  N  T 

197, IQWL(M),M-1,6) 

116  5  2 

R  i  N  T 

198,  (TQWLtMI , H= 1 , 6 ) 

11853 

C-ll 


PRINT  1,99* ITtfCIMf *W*1*6) 

PRINT  iOO, ITffclC(N),M*l,6j 
PRINT  201, ITCH, 

PRINT  20* 

w.  -<», 

80?  DiP  20^ 2 » I DFvTm) •0B4W> . TA^OPA . «MA 
801  D1MEHSW*  A*UiteBU,fTN('TO  ’ 

30»  A9UI-0.C 
AAI2)* SOO, 

&«U|«100» 

BBT2)»50.* 

f*-Q. 

00  90  Nrl,  *N 

tw»tm*ott 

P0  TNtN)*TM 

“u  T"-*T- 

corn  99 

tl»E 


DS«) 


U»«7 
iraS3ia 
1185319 
U85370 
1185320 
'  1485320 

1185320 

1183390 

U85?*e 

11851 

1185390 

11853-90 

1.185390 

1185390 

1185390 
I 185390 
1185391 
1185391 


1X1: 


A  matrix  index  marking  the  boundary  between  the  shelter  air¬ 
space  and  surrounding  earth  (longitudinal  direction) 

1X2:  A  matrix  index  marking  the  end  of  earth,  region  along  the 

longitudinal  direction. 

IY1:  A  matrix  index  marking  the  boundary  between  the  shelter  air 

space  am1  surrounding  earth  (transverse  direction). 

IY2:  A  matrix  index  marking  the  end  of  earth  region  along  the 

transverse  direction. 

IZ1:  A  matrix  index  marking  the  boundary  between  the  shelter  air 

and  roof  earth  region. 

IZ2:  A  matrix  index  marking  the  boundary  between  the  shelter  air 

and  floor  earth  region. 

IZ3:  A  matrix  index  marking  the  end  of  earth  of  shelter  floor 

region. 

CG:  Thermal  conductivity  of  earth,  Btu/hr,ft,  °F. 

n 

AG:  Thermal  diffusivity  of  earth*  ftVhr. 

DT:  Finite  difference  time,  hr 

DX:  Finite  difference  length  along  the  longitudinal  direction,  ft. 

BY:  Finite  difference  length  along  the  transverse  direction,  ft. 

BZ:  Finite  difference  length  along  the  vertical  axis,  ft. 

DTG:  Finite  difference  time  for  computing  the  ground  surface  heat 

exchange,  (modification  of  DT),  hr. 

A,  £,  C:  Internal  dimensions  of  shelter,  ft. 
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Crri:  Shelter  ventilation  air,  cu, ft/rain. 

HV  «  (60)  (specific  heat  of  air)(density  of  air)  *  1.08  for  standard 
air,  Btu/hr,  *F,  CFH. 

WG:  Work  cell  for  ground  surface  temperature  calculation 

if  WG  <  0,  ground  surface  temperature  *  outdoor  air  temperature, 
if  WG  >  0,  ground  surface  temperature  is  computed  by  solar  heat, 
convection  heat  and  conduction  heat. 

RUN:  Run  number , 

2 

HW:  Vertical  wall  surface  heat  transfer  coefficients,  Btu/hr,  ft,  °F. 

2 

HR:  Ceiling  surface  heat  transfer  coefficient,  Btu/hr,  ft,  ®F . 

2 

HF:  Floor  surface  heat  transfer  coefficient,  Btu/hr,  ft,  ®F« 

2 

HG:  Ground  surface  heat  transfer  coefficient,  Btu/hr,  ft,  °F. 

TSW,  TSR,  TSF:  Initial  surface  temperatures  of  wall,  ceiling  and  floor 
respectively,  "F. 

TA:  Shelter  air  temperature,  *F. 

DPA:  Shelter  air  dew  point  temperature,  *F. 

ZN1:  Number  of  600  Btu/hr  occupants. 

ZN2:  Number  of  400  Btu/hr  occupants. 

Zn3:  Number  of  200  Btu/hr  occupants. 

BS:  Sensible  heat  generated  in  the  shelter  by  non-human  source  * 

Btu/hr. 

BL;  Latent  heat  generated  in  the  shelter  by  non-human  source, 

Btu/hr. 

PB:  Barometric  pressure,  in.  hg. 

LI:  If  positive,  initial  ground  temperature  varies  with  depth. 


C-14 


I 


L2:  If  positive,  outdoor  thermal  environment  is  constant. 

L3:  If  positive,  ground  temperature  will  be  printed  out  only  at  the 

end  of  total  time  iteration. 

L4:  If  positive,  ventilation  air  rate  varies  with  time. 

L5:  If  positive,  number  of  600  Btu/hr  occupants  varies  with  time. 

L6:  If  positive,  number  of  400  Btu/hr  occupants  varies  with  time. 

L7:  If  positive,  number  of  200  Btu/hr  occupants  varies  with  time. 

L8:  Strength  of  non-human  heat  source  and  sink  will  vary  with  time. 

NN:  Total  number  of  time  iterations. 

?JPT2:  Work  cell  for  two  different  kinds  of  data  input  formats. 

ZLW:  Lewis  number  for  wall  *  1. 

ZLR:  Lewis  number  for  roof  »  1. 

ZLF:  Lewis  number  for  floor  *  1. 

TG*M :  Constant  initial  earth  temperature,  when  Ll^O,  °F 

TCt«5(K^"  Depthwise  variation  of  initial  earth  temperature  when  L1>0,  °F. 

BB^:  Constant  outdoor  sir  dry-bulb  temperature,  °F,  when  L2^0. 

TV^:  Constant  ventilation  air  dry-bulb  temperature,  °F,  when  L2So. 

QSUN^:  Constant  solar  radiation,  Btu/hr,  £t^,  when  L2^0. 

DPVtf:  Constant  ventilation  air  dew-point  temperature,  °F,  when  L2^Q. 

DB(N) :  Outdoor  air  dry-bulb  temperature  (time  dependent),  °F. 

TV(N):  Ventilation  air  dry-bulb  temperature  (time  dependent),  °F. 
DPV(N);  Ventilation  air  dew-point  temperature  (time  dependent),  °F. 
QSUN(N) :Solar  irradiation  (time  dependent),  Btu/hr  ft^. 

DN(N) :  Ventilation  air  rate  (time  dependent),  cfra, 

ZN1N(N) :Number  of  600  Btu/hr  occupants  (time  dependent). 
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ZN2N(N) 
QN3N(N) 
BSS(N) : 
BLL(N) : 
QX: 

QY: 

QYP.: 

QZF: 

TM: 

STM(N) : 
ST1(N) : 
ST2(N) i 
ST3(N) : 
ST4(N): 
RHA(N) : 
TAA(N) : 
Y(N)  : 


:Number  of  400  Btu/hr  occupants  (time  dependent). 

:Number  of  200  Btu/hr  occupants  (time  dependent) . 

Sensible  heat  due  non-human  heat  source  (time  dependent),  Btu/hr. 

Latent  heat  due  non-human  heat  source  (time  dependent),  Btu/hr. 

Heat  flux  on  wall  normal  to  longitudinal  axis,  Btu/hr  ft  . 

o 

Heat  flux  on  wall  normal  to  transverse  axis.  Btu/hr  ft  . 

o 

Heat  flux  on  ceiling  surface,  Btu/hr  ft  . 

2 

Heat  flux  on  floor  surface,  Btu/hr  ft  . 

Elapsed  time,  hr. 

Mean  surface  temperature,  °F. 

Average  temperature  of  the  surface  normal  to  longitudinal  axis,  °F. 
Average  temperature  of  the  surface  normal  to  transverse  axis,  °F. 
Average  surface  temperature  of  ceiling,  °F. 

Average  surface  temperature  of  floor,  °F. 

Shelter  air  relative  humidity,  %. 

Shelter  air  dry-bulb  temperature,  °F.  -  .  . 

Elapsed  time  =  TM,  hr. 


C  SIMPLIFIED  SYMMETRIC  SHELTER  HEAT  TRANSFER  ANALYSIS  ,T.KUSUCA  118510 

C  RECTANGULAR  PARALLELEPIPED  SHALLOW  UNDERGROUND  118510 

C  S*TOTL  iHELTFR  HEAT  TRANSFER  SURFACE,  SC, FT  118510 

r  CFM-VENTI LATICN  AIR  RATE  ,HV»1.08  SENSIBLE  118510 

CG  -TOTAL  FEAT  GENERATED  INSIDE  THE  SHELTER  U95LO 

C  WG  -WORK  CELL,  IF  WG«0, , CRCUNO  SURFACE-DB(N) , IF  NONZERO, CALCULATE  118510 

C  IF  L1»0  TG0IK5-TGCC  118510 

C  IF  L2-0  OB(N)»DBQ,  TVINJ-TVO,  CSUN(Ni*QSUNO  118510 

C  IF  L3*0  PRINT  TEMPERATURE  AT  EVERY  TIME  INTERVAL  118510 

C  IF  L4*l  READ  CFM  FOR  EACH  TIME  PERIOD  118516 

C  IF  L5-1  READ  ZN1  FOR  EACH  TIME  PERIOC  11851G 

C  IF  L6*l  READ  ZN2  FOR  EACH  TIME  PERIOD  11851D 

C  IF  L7-1  READ  ZN3  FOR  EACH  TIME  PERIOD  U851B 

C  IF  18=1  READ  BS  ANC  BL  FOR  EACH  TIME  PERIOD  I1851B 


DIMENSION  STM  1300), ST i(  300 >, ST2 { 30C ), ST3I 30C > , ST4 1300 ) 

DIMENSION  CPV 1 300 ) , Y { 30C ),EB(2),EC(2!,RMA( 30C ) 

DIMENSION  T(20,20,40),TGO(40),0B( 300) , TV! 300) ,QSUN(300) 
DIMENSION  TAA1300),V(30C:  ,DN(  3CC  ) ,  ZN1NI  3C0 ) ,  Z.M2N  ( 300 )  ,  ZN3N 1 300 ) 


DIMENSION  RSSI300) ,BLL( 300) 

5C0  FORMAT ( 72H  11851B 

l  -  ) 11051B 

READ  5C0  L 1 85 LB 

PRINT  500  11851D 

2C1  FORMAT! 10F7.0)  118510 

202  FORMAT (1017)  118510 

199  REAC  202,1X1, 1X2, IY1,IY2,IZ1, 122, IZ3  113510 

REAC  2Ol,CG,AG,CT,DX,DY,OZ,0TG,A,6,C  118513 

REAC  201,  CFM,HV,WG,RUN,HW,HR,hF,HG  118513 

READ  201,TSW,TSR,TSF,TA,0PA  ,  118513 

REAC  20l,ZNl,ZN2,ZN3,BS,BL,PB 

REAC  202,L1,L2,L3,NN,L4,  L5,L6,L7  118510 

REAC  201,  CPT2 

REAC  20 1 , 2LW , ZL R» ZLF 

IFIL1)  203,203,204  118510 

203  REAC  201.TG00  1 185  1C 

CD  205  K-1.IZ3  118510 

205  TGOI K ) =TGCO  116510 

GC  TO  206  •  118510 

204  REAC  201  TGOIK) ,K  =  1, IZ3)  118510 

206  IFIL2)  2C|7 ,207, 208  118510 

207  REAC  20  L  ,'uBO  ,  TVO,  CSUNO,  CPVU  118510 

DC  209  N- 1 , NN  110510 

OB ( N ) =DBO  118510 

TVINJ-TVG  118510 

OPVtN)«nPVC  11851 

209  GSUN  ( N  )  =  QSllNO  118510 

Gt  TC  210  116510 

208  IF  I  OPT 2  )  401,400,401 


400  READ  402,(CB(N),N=1,NN),(TV(N),N«1,NN),(DPVIN),N*1,NN),ICSUNIN),N- 
11, NN) 

402  FORMAT  I8F9  J) 


GC  TC  210 

401  REAC  20 l , ( CB ( N ) , N-l ,NN) 


REAC  20l,ITV(N),N*l,NN) 

118510 

REAC  201, (CPV(N) ,N-1,NN) 

11851 

IF  (WG)  61,61,60 

118518 

60 

REAC  20!,(CSUN(N)»N*1,NN) 

1185IB 

61 

IF  ( L4 )  63,63,62 

11951B 

62 

REAC  20  1 , ( ON I N ) , N* l , NN) 

1 185  U) 

63 

IF  (L5)  65,65,64 

1 1 85  in 

C-17 


64  READ  70irt2NlNtN),N*l»NN) 

65  IF  516)  67,67,66 

66  NEAC  201?  !ZN2N(M,N«1»NN) 

'BTTFriTl'69r;67,68  ~ 

68  REAC  201, ( ZN3NCN 1 ,N* l,NN 1 

69  RE AC  202 ,  18 

IF ! L8)  210,210,70 

TO  READ  201 , I 3SS (N) ,N*1,NN) 

READ  201, !0LL!N),N=»1,NN) 

210  DTT  30CTK*I,IZ3  '  - 

DC  300  J*l, IV2 

DO  300  T»1VTW  - 

211  IFII-IXI)  212,210,216 

212  IF ( J-im-7TJ72T67716 

213  IFIIZ1-K)  214,216,216 
Z14~IFTR-Tr7)~2T5~,  216, 21  S' 

215  GO  TO  300 

"216  TTIVJ,K)«TGOm - 

3C0  CONTINUE 

PRINT  217, RUN  - 

217  FORMAT ( 50H1  SYMMETRICAL  SHELTER 

pRTfIT-  7rg - 

218  FORMAT ( 65H0  S  CFH 

1  HG  . . . - 

SR*A*B 

"  '  '  SF=SR . - . 

SH*( A+B)«2.«C 

S*SH+SR*5F -  - — 

PRINT  2l9,S,CFM,TA,hV,HW,HG 

‘65 0'  F  C  R  MAT  riO  0  HO  HG  - HH 

l  2LR  ZLF  CX 

PRINT  650  . -  —  - 


HEAT  TRANSFER, RUN 


TA  HV 


HR  -  HF 
OY  D2 


118518 
118518 
118516 
'  118518* 

118518 
118518 
118518 
118516 
118518 
118510 
118510 
118510 
11851C 
118510 
118510 
118515 
118510 
_  uf510 
118510- 
1X8510, 
F3.0  1118510 

rrmo 

HH 

1 118510 
118513 
116513 
118513 
'  n85I3 

TJ1 - 

1 


651  FORMAT! 10F10. 31 

Tin  NT  651,  HG,  HU,  HR,  HF ,  Z  L¥77tR7ZrF7DT,T)Y ,  DT 

219  FORMAT (6F10.3) 

PRTNr  220 - - 

220  FORMAT! 55H0  DB(N)  TV!N1  QSUNtM) 

DC~722  N*T7NN - - - 

222  PRINT  221,  DB ! N 1 ,TV IN) ,QSUN 1 N ) 

221  F0RMAH  JH0.21 - - - 

READ  202,11, J3,KK 

C  END  OF  DATA  INPUT 
KX*IX2-l 

KJOt-TXl-1  -  - - -  - - - 

KY-IY2-1 

^yyrmr-1 - - - 

KZ- 123-1 

' ”  xzi-izi-i  - . — •  - - - - - 

KZ2-IZ2-1 

EX«*G*OT/(DX*CXJ  -  - 

EY«AG*OT/!OY»CY) 

Fr-xwrr/TCTtfcm - - 

Hl«HH«OX/CG 

H2*HH«TJT/CG - - - — - — - - - 

H3  «HR«CZ/CG 

. '  — iff  •HF»cr/cr  - -  - - - - - 

HGU*HG*02/CG 

RtrarTV-rm  1 X  T  E  YTOT - - 

RX«R0-2.»EX»H1 

-  RY«R0-2.»ET»H2  - — - - 


118510" 
118510 
1118510 
_  U05TC 

_ _ 118510 

1T85T5 

118510 
118510 
* "  118510 
118510 
-  It85TU- 

118510 
T18510 
118510 
— n¥5T0 
118510 

rrenro 

118513 

~  nun 

118513 

Tiisn 

1 1-8510 

rnrrro 


C-18 


RZR*R0“2«*EZ*M3 

"  R2F*ftO-2.*ET*W5 -  - 

RG«1 .-2. •(  £Z  )  • 1 1 .  ♦HGU)*CTG/OT 

/  '•  TK-fl.  ■  .  "  '  . 

V,'  END  OF  INITIALIZATION 

PRINT  225  . . . 

223  FGftFAT  1 lOOHO  RO  RX  RY 

1  HI  W2  . H3  . . H4  -  - 

PRINT  229»R0,RX,RY,RZR,RZF,H1,H2,H3,H4 
329  PRINT  228 

22 8  FCRPAT I I05H  CG  AG  KGU 

1  OZ  CT  EX  c»  cl 

PRINT  229, CG, AG«HGU  ,OX , OY.OZ. CT, EX, EY.EZ 

229  FCRPATH0F10.3I - 

234  FCRPAT(20F6.1) 

XX*IXl  .  . . 

YY-IY1 

ZZ-IZ2-IZ1 

AA»CX*XX«2. 

8B*CY*YY«? . 

COCZ*ZZ 

SX«AA*CC  * . .  .  ■* 

SY«88*CC 

SZR*AA*B8 

SZF-SZR 

612  FCRPATC  50H1  QX  QY  QZR 

PRINT  612 

KA«Nft)PA)  -  -  -  "  " 

NH>h(TSU) 

•UR*h(TSR) 

WF-MTSF) 

C  START  OF  TINE  ITERATION 
00  907  N-l.NN 

. .  rF  IL4)  51Y51,5t'  ••*** 

30  CFN  •  DM  NJ 
51  IF  ( L5 )  53,53,52 
52  ZN1  -  ZN1NIN) 

53  IF  <L6)  55,55,54 
54  ZN2  »  ZN2MN) 

55  IFIL7)  71,71,56  - 

56  ZN3  •  ZN3NIN) 

71  I F l L8 )  57,57,72  - - - 

72  6S«8SS(N) 

BL-8LWN)  . . 

57  TP«TH*OT 

. . TMV»-fcWT5TJtrTSm,6Cr2~"  - -  ~ 

602  ZZW106t.«tHA-MNi/(  TA-T$W*O.OOOC l J /O. 243 

- rzn*zrv7zur  ~  — - - * — 

GO  TO  603 

601  ZZWO.  ~ 

603  IF(MA-UR)604,604,605 

605  2ZR*1T)61 .  *  t  NA-WRT7Tr*=TSR  *XT,  COCCI  I  AJ72T3 - - 

ZZR-ZZR/ZLR 

GO  TO  606  - - - - - 

604  ZZR-O. 


118510 

n¥5ie 

118510 

- HH5I0 

RZF  118510 


118510. 

1118510 

118510 

118510 


TN - 1118513 

118513 

* - nesTT 

11851. 

- ntsn 

118513 

-  - 118510 

118510 

- - LJiSTI 

118518 

- nnmo 

118518 

- rrora 

118518 

- TT851H 

118518 

- 1I55T8 

118518 

mm* 

118510 


118513 


J|&06  IFINA-KF >607,607,608 
W60«  ZZF»1061.*(NA-HF)/ITA-T 


SF *0.0000 1  J/0.243 


-  ZZFalZF/ZLF 
GC  TO  609 
607  ZIF*0.  - 


118513 

118513 

~imrr 

TFtfSTT 

t£8513 
_ 118313 

110513 


C-lt 


60 v  TSM*!TSh*SMTSR»SK+TSF#$F)/!S*  +  $R+SFj  115513 

P«*TA 
PR-TA 

?F«TA  •  -  ~  '  '  "•  . .  ~  .  '  ' 

NX*Hl*<  1,+ZZW)  118513 

HY-H2*! l.+ZZW)  118513 

HZR-H3»< l.+ZZR)  118513 

HZF*H4»  t 1 . +ZZF )  118513 

0UU*CSUN(N)*02/CG  118510 


QS*CG(TA , ZN1 1 ZN2 , ZV3,1. 1 +BS 
TAi«TA*l. 

CSl*SG(TAl,ZNl»ZN2tZN3, 1.  )*8S 
V1*CS1-CS 
Y2*CS*TA1-CSI»TA 
CL«CG(TA,ZN1,ZN2,2N3,0. I+BL 
VZV*WtOP-/(Nn 
RX*R0-2  ■  »EX»HX 
RY«R0-2.»EY«HY 
RZR*RO-2.»fZ*HZR 
RZF*RO-2.»EZ»HZF 
7C2  I F ( L 3 ! 23 l » 231 i 845 
645  IffN-NN)  245,23i;21I 


231  PRINT  2 32  118511 

232  F0RPATI50HI  EARTH  TEMPERATURE  ON  Z-X  PLANE  AT  J*1  1118511 

OC  235  K«l,IZ3  118511 

215  PRINT  234,  mi,  1  ,X  )  ,  !«  I,  IX2J  - 118511 

PRINT  236  118511 

236  FORMAT! 50HI  ~  EARTH  TEMPERATURE  ON  Z-X  PUNE  *T  J»JJ  ■'1118511 

00  2 37  K» 1 , 1 Z  3  1 1851 1 

23T  PRINT  234, (Tir.JJ.KJ, 1*1, 1X21  '  "111511 

PRINT  238  116511 

238  FORMAT f 50H1  EARTH  TEMPERATURE  ON  Z-Y  PLANE  AT  1*1  1118511 

00  239  K« 1 , I Z  3  118511 

- rJTTRnmiA.ITU.J.KIVJ-l.IYZl  -  118511 

PRINT  240  1 185  l i 

240  FORMAT (50H1  EARTH  TEMPERATURE  CN  Z-Y  PLANE  AT  I-II  1118511 

00  241  K-I.IZ3  118511 

241  PRINT  234 ,( T II T ,J ,K 1 ,  J«l,  IY21  116511 

PRINT  247  118511 

. .  T42-  FORMAT  1 50Hf  EARTH  TEMPERATURE  ON  X-Y  PLANE  AT  '  K«XK  -  "  1118511 

OC  243  J«1,1Y2  118511 

243  PRINT  234,  IT!  I , J.KK  ) ,  1« l ,  1X2  MS511 

C  EARTH  TEMPERATURE  BOUNDARIES  H8511 

245  00  103  1»1, 1X2  118511 

OC  103  J-1.1Y2  118511 

. . IF!  VG 1  101,101,102  '  “  - -  1 18511 

ICt  T(I,J,il-OB!Ml  118511 

GO  TO  IC3  ■  '  118511 


102  TG  »  0. 

8C2  T(1 ,  J,  l)*2.»£Z*lTi  I  »J ,  2 )  »HGU»D6 !  N )  ♦QUL  >  •OTG/DTtT  1 1,"3»  1  }*RG  - 

If  ITG-Cfl  803,103,103 

8C?  TG*TC*DTC» .  . *  '  ‘  . .  “  -  " 

GO  TQ  802 

IC3  TtT,  J,rZ31«TCmiZ3»  -  -  1 16511 

TGOdl-T!  iX2,lY2,l! 

UNDISTURBED  EARTH  TEMPERATURE 

OC  104  r«2,RZ  1 1851 1 

”1  C 4  TCO ( K 1  * f Z •  C TCSTK-T)  > TGO TK ♦ 11 J  VTG  0 ( XT*  1 1  r-^TTSfTT  ~  - - - - “TmTT 

CO  105  K*2,K2 
DO  105  3«3,IY2 


116511 


105  T  <  IX2,J,K}«TG0(K) 

- DC  1D6  K»2*K7  . . 

00  106  I-1.IY2 
TDF  T(I,  rY2tKJ«TC0(!O 
DC  900  K*?,*'* 

DC  900  J«l,KV 

00  900  I  *  1 1 KX 

P2«TlI*l,J,K)  '  •  - - 

P4»T  (1 ,  J+l  ,K  ) 

‘  P5»TT1 » J,K-1 ) 

P6»T ( I , J,K. 1 ) 

IF(I-t)  1,1,2 

1  P1-P2 

- GO  TO  3  —  ~  - -  '  •  - . 

2  Pl-TU-l.J.K) 

— rTFCJ=iT~¥.4Y5 - - * -  *•  "  - 

4  P3»P4 

GO  TO  6 . ~  '  "  .  . 

5  P3*T ( I , J-l ,K ) 

6  P0*TIT*  J  ,K  J 

ENO  OF  PREPARATION 

- rFrr=m7“T,n;vi - — . 

7  IF(J-IYi)  8,11,11 

8  IFUZ1-K)  9,11,11 

9  IF( K~ 1 22 )  10,  1 1 ,  11 

10  GO  TO  900'  ' 

11  IF (K-lZl 1800, 12,19 

“iziFii-ivn  rvi5,80ff 
13  IF(J-IY1>  14,18,800 

i  14  P0>»£X*{Pl»P2)+EY»lP3*P4U2.«£ZMP5*HZR*PR).P'0»«2R 
'  GO  TO  900 

15  1FU-IYU  16,17,300 

16  P0»<  E*/3.)«(2.»Pl*4.«P2)» 

- -  CTTO  9G0 . ‘ 


£Y<M  P3.P4 1+IEZ/3, }•( 4,«P5*2.*P6>  +  P0*RC 
IEV/7. l*(6.*P3*8«*P41*(EZ/7, )»(8.«P5*6. 


17  P0-(EX/7.1*16.«P1+8.«P21> 

16T*P0«R0 
GO  TO  900 

11  PO*£X.  ( P1+P2 )• I  EY73.  )«(2.*P3*4.«P4)4‘|£2/3,1»(  4  ,»P5*2.'*P61*P0*R0 
GO  TO  JOO 

~rrTFnir-xr  800,26, tot —  — - - -  - - 

20  I F I  I— 1X1)  24,21,  800 

21  IFU-1Y1)  22,23,  SCO  - 

22  P0-EY.(P3.P4)*£Z«tP5.P6)* 

- GO  TO  9C0 

23  P0"(EX/3.)» ( 2,*P l *4, «P2 ) ♦ 


118511 


118511 
118511 
118511 
110511 
118811 
118511 
118511 
118511 
116511 
'118511 
118511 
“118511 
118511 
1 18511 
118511 
118511 
118511 
“118511 
118511 
118511 
118511 
118511 
118511 
118511 
118511 
11851? 
11851' 
118511 
118511 
118511 
•Pi 13511 
1 16511 
113511 
It 85  1 1 
118511 


>2.»£X»IP2*HX^PW).P0»RX 

-IEY/3.  )*I2,«P3*4.»P4JHE2*IP5+P6)*P0*R0 


24  IF ( 

rr"PD» 

GO 

‘  26  IF ( 
2?  I F ( 
2FWi 
GO 

29  PO- 

GC 

*10'  IF! 

31  PO* 

- GO 

32  PC- 
T1)- 


Td  900 
J-l Y1 1900,25,  800 

FX *T PTVPTT. EZ  *TPT.T6J  ♦  7.4TY.T  PFWY  *FB  J  rP0VR¥ — 
TO  900 

I-IXU  77,30,800  - - ““ - 

J-1Y1 >  28,29,800 


■ )  ♦PDiRTF - - — 

TO  900 

ex.  ( p  1.P2 1  vrmrn  «r?;  *T5** .  •  p*j  ♦tei  n .  r*iz.»vst*:i*sr+p(nm 

TO  900 

J-IVU  31,32,800  . 

( EX/3.  > •  1 2.»Pl»-4, »P2  )*EV«I P3*P4l-*IE2/3.  l«l 2«*P5*4. •P61*P0»R0 

rrm - ~ - - - - - ? — - * — 

IEX/7. >*<6.*Pl*8.*P2 )*(€Y/7. !•( 6.«P3«8.*P41 ♦IEZ/7. )M0.*P6+6. 
PifiFD - * - * - - - — 


118511 
118511 
118513 
118511 
118511 
118511 
118511 
118513 
118511 
T18511 
118511 
ITi5l  3 
118511 

rmu 

118511 

118511 

1195U 


nwri 

•P118511 

'  T185U 


c<l 


CC  TO  900  118S11 

800  POEX«  tPl«-P2  )*EY*l  P3«-P4  1*EZ«  (  P5+P6  )*P0«RO  IlfSli 

9C0  T (  1 1  J«K 1»P0 

910  zt-o.  . - . .  "  ”  - -  Tfi?r 

ZN*0 •  11851 

DC  911  K*'Zl,tZ2  '  .  itffsr 

DC  911  J*l,lYl  11851 

ZT»?T*TUX1,J,K>  11851 

911  ZN*ZN*1.  11851 

TSX*ZT/ZN  11851 

ZT«0.  11851 

ZN*0.  11851 

DC  912  K-IZ1.IZ2  11851 

OC  912 - T-T7TX1  11851 

ZT-ZT*T< f.IYl.K)  11851 

‘912  ZNiZNTT.  .  11851 

TSY-ZT/ZN  11851 

-  Z>T*0.  -  11851 

ZT-0.  11851 

oc  9i  j  1*1, ixi  nasi 

00  913  J*l,IVl  11851 

. ~7T»rr*rn,j,izn  "  .  '  uesi 

913  ZN*ZN+1.  11851 

TS7R-ZT/ZN  lie5l 

ZMO.  11851 

'  ZT*0.  -  na3l 

DO  914  1-1, 1X1  11851 

- off—  9n — J«lf  1Y1  •  .  '  11851' 

ZT*ZT+T ( 1 , J, IZ2 )  11851 

T15ZN*ZN+i:  T1851 

TSZF*ZT/ZN  11851 

- -q x«HU«  T  '  .+2ZH)  •fPH— TSTJ . . .  ~  118513* 

QY*HW»  ( l .  +  ZZW 1 *( PM-TSY)  118513* 

- 0Z5iTHr»TT.^ZZrr?rPR-T5ZST  1 18515* 

CZF«*HF*U.»ZZF1*MPF-TSZF)  110513*- 

- TW*rr5X'*5X4Tyy^r/rSJr<5TT~ -  118513* 

T5F*  SZF  118513* 

- TSR.TTZW  -  "  .  -  -  ^15,3, 


T$M*( T  Sta*5W*TSR«SR*TSF«SF )/(SM+SR+SF)  118513* 

~TA*(  1T0T»TFM»  TVTNTVY2WH  *  T  S  W • SW* R  R .T  S  R • S  ftVHF  •  YST.5  F  PA1 1 . 0  8«rF  R^Yl  ♦  1 1 8  5  lT 


lHfc»5N*HR»SR*MF«SF ) 
0L-CC(TA,ZNl,ZH;tZN3,0.  J+Bi: 
UW*H( TSM) 

WR*h(TSRl 

MF*W  t  TSF ) 

■W-Kk/Zltf'  ■'  - 


118513' 
118513’ 
118513* 
1 185131' 
118513* 


HR*MR/ ZtR 

HF*HF/ZLF  . . -  . 

HA*(4780.«CFM«WV*Ql*( 1061 . 1 • ( 6W • Sta •HW*WR»SR *MR SF*HF 1 /0.243) / « 41 1 851 3* 
1760. *CFP*(106l. >*<SW«HU*SR»Hir»ST»HF 1/C. 7431  '  11851T 


HW*HW«ZIW 

"  HR*KP«ZU 


HF«*HF*ZIF 

HS-k(TSP)  -  . 

IF(KA-MS)  610,610.611 

610  tfA.t470a..CFP»WV+CUm78O.«Xf*T~' 

611  PSA  •  WA«P8/IUA«0.622) 

— HHTTNi  *i  wi'irsKretrt  n  j - 

613  FCRMT15F10.2) 

- PI  IWT51TV0T,  QY‘,  0nT707F7TH - 


1185134 
— -Ttim* 

1185131 


1185134 
'  11*5131 


o «n 


STP(N)«TSM 
STi(N)*TSX 
ST2INl«T$r 
ST31N)«TSZR 
$T4(N)*T$ZF 
V(*a«TM 
Y(N)*TH 
907  TAA(N)»TA 
PRINT  950 

950  FCRPAT(  100HI  t*  STM  Stl  $T2  ST3 

i  ST4  RHA  TAA 

OC  951  N*1,NN 

951  PRINT  954,YtL),STN<N),STl(N)  ,S72(NJ ,  ST3J N »  ,  STMN )  ,RHAJ N >  ,  T AA< Nl 
954  FORMAT ( 6F10.2 ) 

E»m*o.‘ 

E8( 2 )"500. 

Ecm-ico. 

EC(2)*40. 

CAU  PLOT (3,NN, Y, STM,NN ,V ,TAA,2, £8,EC ) 

901  CGNTINUE 
CALL  SYSTEM 
GC  TO  199 
ENO 


1185121 

118512' 

"1185121 


C-C3 


Program  M--4  Input  aad  Output  Symbol* 


NN:  Total  number  of  time  iterations. 

II:  Number  of  ventilation  air  inlets  where  temperatures  were  observed, 

JJ:  Number  of  shelter  air  temperature  stations. 

KK:  Total  number  of  shelter  inner  surface  temperature  stations. 

KWC:  If  ^0,  Ft.  Selvoir  1000-man  shelter  analysis. 

If  >0,  Ft.  Belvoir  200-man  shelter  analysis. 

WC:  If  £0,  constant  earth  temperature,  TG^,  is  read  in. 

If  >0,  six  directional  earth  temperature  profiles,  TG(I,M), 
will  be  read  in. 

AI:  If  >0,  air  conditioning  calculations  are  included. 

WC1:  If  —0,  cfta  varies  with  time. 

WC2:  If  —0,  number  of  occupants  varies  with  time. 

TG^:  Constant  initial  earth  temperature,  °F. 

TG(I,M) :Six  directional  initial  earth  temperature  profiles,  *F. 

TC^(M);  Initial  inner  wall  temperatures,  ®F. 

TC(L,K) : Initial  temperature  profile  within  the  shelter  inner  walls,  °F. 
S(K):  Shelter  inner  surface  area,  ft  . 

DD(K) :  Finite  difference  length  for  inner  wall,  ft. 

DG(K) :  Finite  difference  length  for  earth  rsglon,  ft. 

D(F) :  Thickness  of  the  inner  wall,  ft. 

ZL(1C.) :  Thickness  of  the  earth  block,  ft. 

AG(K) :  Thermal  diffusivity  of  earth,  ft^/hr. 

AC(K):  Thermal  diffusivity  of  inner  wall,  ft^/hr. 
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CK(K): 

CG(K): 

H(K): 

HD(K): 

DT: 

DTT: 

GS: 


GL: 

GA: 

FR: 

CF: 

DATT: 

TW1: 

XDB(N) : 
TV(N) : 
DPV(N) : 
EG(K) : 


E(K): 


EC(K) : 


TM: 


Thermal  conductivity  of  inner  wall,  Btu/hr  ft,  °F. 

Thermal  conductivity  of  earth  block,  Btu/hr,  ft,  °F, 

Heat  transfer  coefficient  of  the  inner  surface,  Btu/hr,  ft' ”,  °F. 
Vapor  transfer  coefficient  of  the  inner  surface,  Ib/hr,  ft  (lb/lb)  . 
Finite  difference  time  for  inner  wall  temperature  calculation,  hr. 
Finite  difference  times  for  earth  temperature  calculation,  hr. 

Per  capita  sensible  heat  generation  from  non-human  heat  source, 
Btu/hr,  person. 

Per  capita  latent  heat  generation  from  non-human  heat  source, 

Btu/hr,  person. 

Recirculation  air  rate  for  conditioning,  cu  ft /min. 

Coil  effectiveness  of  the  air  conditioner  (dimensionless) . 

Coil  contact  factor  of  the  air  conditioner  (dimensionless) . 

Coil  leaving  air  temperature  rise  due  to  fan  heat,  °F. 

Entering  temperature  of  the  coolant  to  the  air  conditioner 
coil,  °F . 

Outdoor  air  dry-bulb  temperature,  °F. 

Ventilation  air  dry-bulb  temperature,  °F. 

Ventilation  air  dew-point  temperature,  °F. 

Finite  difference  stability  modulus  for  earth  temperature 
calculation,  which  should  not  exceed  0.5. 

Finite  ’ifference  stability  modulus  for  shelter  inner-wall 
temperature  calculation,  which  should  not  exceed  0.5. 

E(K). 

Elapsed  time  (hr) . 


C-25 


p,VS:  Sensible  heat  carried  out  by  %'entilation  air,  Btu/hr, 

QVL:  Latent  heat  carried  out  by  ventilation  air,  Btu/hr. 

QG5:  Sensible  haat  generated  in  the  shelter,  Btu/hr. 

QGL:  Latent  heat  generated  in  the  shelter. 

QW5T;  Sensible  heat  absorbed  by  shelter  surface,  Btu/hr, 

QW'LT:  Latent  heat  absorbed  by  shelter  surface,  Btu/hr. 

TWCT:  Condensate  collected  in  the  shelter,  lb/hr. 

CAS:  Sensible  heat  absorbed  by  the  air  conditioner,  Btu/hr. 

QAL:  Latent  heat  absorbed  by  the  air  conditioner,  Btu/hr. 

TA:  Shelter  air  drv-bulb  temperature,  °F. 

DPA:  Shelter  air  dew-point  temperature,  °F. 

WRA:  Shelter  air  wet-bulb  temperature,  °F. 

RHA:  Shelter  air  relative  humidity,  %, 

EFS1:  Shelter  air  effective  temp-  mature,  CF. 

TSM:  Average  shelter  inner  surface  temperature,  eF. 

QWS(K):  Sensible  heat  transferred  to  the  Kth  surface..  Btu/hr. 

QWL(K):  Latent  heat  transferred  to  the  Kth  surface,  Btu/hr. 

TWC(K):  Condensate  collected  onto  the  Kth  surface,  lb/hr. 

FLX(K);  Heat  flux  of  the  Kth  surface,  Btu/hr,  ft^. 

TS(K) :  Average  temperature  of  the  Kth  surface,  °F. 

TC(L,K) : Temperature  distribution  In  the  Kth  inner  wall,  rF. 

TG( I, K) : Temperature  distribution  in  the  Kth  block  of  the  earth. 

TIME(N) :Elapsed  time  for  the  observed  date,  hr. 

CFM(N) :  Ventilation  air  rate  used  in  the  experiment  (time  dependent),  cfm. 
ZN(N);  Number  of  occupants  used  in  the  experiments  (time  dependent). 
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XI: 


■  DBV(N)  -  Observed  ventilation  air  dry-bulb  temperature,  *F. 
X2:  "  WBV(N)  -  Observed  ventilation  air  wet-bulb  temperature,  *F. 

DBS:  Observed  shelter  dry-bulb  tenp  rature* 

WBS:  Observed  shelter  wet-bulb  tenders tore. 

STS1:  Shelter  effective  temperature  based  upon  DBS  and  WBS. 

Subscripts: 

MwK  =■  Type  of  wall  exposure. 

1.  North  wall 

2.  South  wall 

3.  East  wall 

4.  West  wall 

5.  Floor 

6.  Roof 


N  *  Time 

L  =  Spatial  matrix  index  for  inner  wall 

I  «*  Spatial  matrix  Index  for  earth  region 


Subroutines : 

{1)  FREP:  This  subroutine  is  designed  to  accept  observed  data  of 
prototype  shelter  and  yield  the  psychrometrically  processed  mean 
shelter  conditions  in  terms  of  temperature,  humidity  and  heat  flux. 
Following  are  the  input  data  specific  to  this  subroutine. 


KWS: 

The 

first 

NWS: 

The 

total 

XWS(I): 

Saturated 

KM; 

The 

first 

The  first  temperature  entry  for  XHA  table. 


C-27 


NHA; 
XHA(I) : 
KHS: 
NHS: 
XHS(I): 
KFS: 
NFS; 
XFS(I) : 

KHW: 

NHW; 

XHW(I): 

KVA; 

NVA; 

XVA(I) ; 

KVS: 

NVS; 

XVS(I) : 

KQS: 

NQS; 

XQS(I); 

EDB(I) : 

EWB(I) : 


The  total  number  of  temperature  entries  for  XHA  table. 

Dry  air  enthalpy  table  of  Goff  and  Gratch,  Btu/lb. 

The  first  temperature  entry  for  moisture  XHS  table,  ®F. 

The  total  number  of  temperature  entries  for  XHS  table. 

Moistiite  saturated  air  enthalpy  table  of  Goff  and  Gratch,  Btu/lb. 
The  fi-st  temperature  entry  for  the  XFS  table. 

The  total  number  of  temperature  entries  for  XFS  table. 

Factors  relating  the  degree  of  saturation  to  the  relative 
humidity. 

The  first  temperature  entry  for  the  XHW  table. 

The  total  number  of  temperature  entries  for  XHW  table. 

Liquid  water  enti.alpy  table  of  Goff  and  Gratch,  Btu/lb, 

The  first  temperature  entry  for  the  XVA  table. 

The  total  number  of  temperature  entries  for  XVA  table. 

Dry  air  volume  table  of  Goff  and  Gratch,  cu  ft/lb. 

The  first  temperature  entry  for  the  XVS  table. 

The  total  number  of  temperature  entries  for  XVS  table. 

Moisture  saturated  air  volume  table  of  Goff  and  Gratch,  cu  ft/ib. 
The  first  temperature  entry  for  the  XQS  table. 

The  total  number  of  temperature  entries  for  XQS  table. 

Per  capita  sensible  heat  table  of  human  body,  /tu/hr. 

Dry-bulb  temperature  entries  to  ASHRAE  effect' ve  temperature 
table,  “F. 

Wet-bulb  temperature  entries  to  AiSHRAE  effe  ;tive  temperature 
table,  °F. 
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(2)  STRCQN:  Subroutine  to  confute  adiabatic  shelter  air  conditions  by 
accepting  ventilation  air  dry-  and  wet-bulb  temperatures  and  ventilation 
air  flow  rate. 

(3)  AIRCgSN:  Subroutine  to  compute  sensible  and  latent  heat  absorption 
by  accepting  the  following  input  data. 

GA;  Recirculation  air  rate,  cfm. 

CFM:  Ventilation  air  rate,  cfm. 

DBS:  Shelter  air  dry-bulb  temperature,  °F. 

WBS;  Shelter  air  wet-bulb  temperature,  ®F. 

DBV;  Ventilation  air  dry-bulb  temperature,  °F. 

WBV:  Ventilation  air  wet-bulb  temperature,  °F. 

TW1:  Inlet  temperature  of  the  air  conditioner  coolant,  ®F, 

E:  Effectiveness  of  the.  air  cooling  coil. 

CF:  Contact  factor  of  the  air  cooling  coil. 

DTA:  Dry-bulb  temperature  rise  of  the  air  conditioner  outlet,  air 

due  to  the  fan  heat,  °F. 

(4)  FN:  Boundary  temperature  calculation  subroutine  by  a  finite 
difference  formula. 

(5)  EFTI;  Subroutine  for  computing  effective  temperature  by  using  the 
input  data  of  dry-  and  wet-bulb  temperatures. 

(6)  DBWBH:  Subroutine  for  computing  enthalpy  of  moist  air  using  dry 
and  wet-bulb  temperatures. 

(7)  DBWBRH:  Subroutine  for  calculating  relative  humidity  for  given 
dry-  and  wet-bulb  temperatures  and  barometric  pressure,  PB  (in.  Hg). 
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(8)  DBWWBH:  Subroutine  to  calculate  wet-bulb  temperatures  and  enthalpy 
of  moist  air  based  upon  given  values  of  the  dry-bulb  temperature  and 
humidity  ratio. 

(9)  PV:  Subroutine  for  calculating  the  vapor  pressure  in  inches  of  Hg 
of  saturated  moist  air  at  given  temperatures, 

(10)  DBWBDP :  Subroutine  for  calculating  the  dew-point  temperature  and 
humidity  ratio  of  moist  air  for  given  dry-  and  wat-bulb  temperatures, 

(11)  DBWPWB :  Subroutine  for  calculating  the  wet-bulb  temperature  «-nd 
humidity  ratio  when  dry-bulb  temperatures  and  dew-point  temperatures 
are  given. 

(12)  TBLU:  Table  look-up  and  linear  interpolatior  subroutine  when 
the  independent  variable  is  incremented  by  unity. 

(13)  XTBLU :  Table  look-up  and  linear  interpretation  subroutine  when 


the  incrementation  of  the  independent  variable  is  non-unity. 


nooooonooooooooonooooooooooonooo 


ONE  DIMENSIONAL  SHELTER  11851003 

SUBROUTINE  TO  BE  USEO  IN  THE  PROGRAM  1 1851004 

PREPIDT.II.JJ.KK.KMC.NN)  I1851M03 

STRCON 1 08, MB, G)  ADIABATIC  SHELTER  AIR  CONOITION  CALCULATION 

AIRCONIGA.CFM.DBS.MBS.DBV.MBV.TWt , E.CF,  QAS.QAL.OTA)  11851A02 

FUNCTION  FN 


EFT  1 1 DB ,M3 »EFX ) 
OBMBHIOB.MB.H) 
DBMBRHI DB.MB.PB.RH) 
DBMMBHI DB ,M?  MB ,H) 
FUNCTION  PVIX.PB) 
OBHBOP ID8.W8.0P.M) 
DBOPMBIDB.DP.MB,  w) 
T8LU  IX.KX.V.Y) 


GIVEN  OB  AND  MB  .COMPUTE  H  11851012 
GIVEN  OB.MB.PB  .COMPUTE  RH  11851014 
GIVEN  OB  AND  M  .COMPUTE  MB  AND  H  11851011 
CALCULATE  VAPOR  PRESSURE  FOR  GIVEN  OB  AND  PB11851018 


OBWBOPIDB.WB.OP.M)  GIVEN  OB  AND  MB  .COMPUTE  DP  AND  M  11351010 

DBOPMBIDB.DP.MB,  W)  GIVEN  DB  AND  DP,  COMPUTE  MB  11851013 

TBLU  IX.KX.V.Y)  KX-FIRST  VARIABLE  11851015 

V-TABLE  11851016 

Y-VALUE  FOR  X  11851017 

XTBLUlX.OX.KX.LXtV.Y) ,  OX- VAR  I  ABLE  INCREMENT  11851005 

KX-FIRST  VARIABLE  11851006 

LX-NUM8ER0F  ENTRIES  11851007 

V  -TABLE  11851008 

Y  * VALUE  OBTAINED  FOR  X  11851009 

1 -NORTH, 2- SOUTH, 3-EAST ,4«MEST,5»FL0QR, 6-ROOF, J-ROOF  MATRIX  11851008 

AI *0  -  NO  AIR  CONDI T IOING.  +-AIR  CONDITIONING 

NC-0-,  CONSTANT  TGO  .^-PROFILE  11851M12 

II  -NUMBER  OF  VENTILATION  TEMP.  DATA  11851004 

JJ  -NUMBER  OF  SHELTER  TEMP.  DATA  ,  11851005 

KK  -NUMBER  OF  MALL  TEMP.  DATA  11851006 

KMC  -♦  200MAN  INPUT  ,0-  -1000  MAN  INPUT  11851007 

MCI—  OR  0  CFM  VARIES 
MC2—  OR  0  ZN  VARIES 
JT-  FIRST  ESTIMATE  OF  CONCRETE  TIME 
OTT*  EARTH  TIME 

OIMENSICN  XMS ( 150I.XHAI I50).XHS( 1 50 ) , XHMI 150 > , XFS 1 150 1 , XVA 1 1 501 . XVI 1851H04 
1SI 150) , QS! 50 ) «CFM! 300) ,ZNi  300 1 » 

2DBV I  300 ! , OPV 1 30C I , M6V ( 300 ) * 08C I 3 , 300 > .M8CI2 , 300 ) , PS  1 10, 300 1 , 

3T VI  3001  »QPC I  300 1 , LA  C 300 . 10  I . LBI 300 ,10),X0BI300),EDB!100),EM8(100), 

4ETI 50.50  I .  TMCI6) , EC  1 6}, EG! 6 1 .TSI6J .TSSI6) 

5  ,  TGI  10,6 J , TCI  5. 6) ,S(6i , 91 6 > , DO! 6 > , ZL I  6 J , DGI 6 > , AGI 6) , AC  I  6 ) 11851004 

6,  CK 1 6), CGI  6) ,H<6! .H0I6I,  OHS  I  6 ) , GMt I C > , TQMS 1 6 > , TQNL 1 6) , TQMI 6 >11851005 

7,  DXI6) »CSI6) r  QM ( 6 ) , FLX I  6 ) ,EI6)  , AS  1 6 > , ES ( 6  I , TCOI 6 ) , T IME I 30U ) 

COMMON  XMS,XHA,XHS,XHM,XFS,KMS,KHA,KFS,KMM,KHS,E0B,EMB.ET,  KET »N1 1851N09 
lfcTIQS,C0l,CO2,U0,TM.GS,GL,2LEMS,CFM,DPS,0BS,M8S,ETSl,ETS2,WCC,CO3,ll851M10 
2MCD,DPA,ZN,QG,0GL»D8V,0PV» DBC , DPC  *  MRC «  S , TA , DT , KMC , NN , X VA , 
3XVS»XVA»KVS,PS»LA,LB.X03»M8V,  MCI.MC2  > DTT 
77  F0RMAT1 IHOIOFIO. 1 ) 

1  FORMAT  I 10F7.0)  "  11851009 


2  FORMAT  I  1017) 

12  FORMAT  I 43HI 

13  FORMAT! 43H0 

14  FORMAT  I 43H0 


UNDERGROUND  FALLOUT  SHELTER  DATA 
INITIAL  EARTH  TEMPERATURES 
INITIAL  CONCRETE  MALL  TEMPERATURE 


11851009 

11851010 

11851057 

11851058 

11851059 


15  FORMAT  I 43HI  TIME 

T6'7l'GMAT'  43K0  081  NI 

VARIABLES 

TVIN) 

DPVINI 

QSUNIN! 

) 

) 

11851060 

11851061 

17  FORMAT I 4F 1 0.  1 ) 

18  FORMAT! 63H  TGI1) 

TGI  2) 

TGI  3  ) 

TGI  4 ) 

T5I  5 ) 

1 1851062 
TGI  1851063 

116) 

19  FORMAT I 63H  TCI  1) 

TC(2) 

TCI  3) 

TCI4) 

TCI5I 

>11851064 
TCI 1851 065 

iib: 

20  i  URMAT16M0.1) 

124  FORMAT 1 43HI  PHYSICAL  OATA  USED  FOR 

COMPUTATION 

) 

>11851066 

11851067 

11851081 

125  FORMAT  1 110H0SIK) 

OIK) 

ZL IV  ) 

CK(K) 

CG1K) 

ACIKi 1851082 
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1)  AG ( K  S  HI K) 

HDfK ) 

)) 1851083 

126  FORMAT I 10F10.4) 

11851084 

64  FORMAT! 16H0  QWS(M) 

6F15.1) 

11851206 

45  FORMAT! 16H0  QHLIM) 

6F15.1) 

1185120? 

6  FORMAT 1 16H0  T  WC ! M  ) 

6F15.1) 

11851208 

67  FORMAT! 16H0  FLXIM) 

6F15.1) 

11851209 

68  FORMAT! 16H0  TS(M> 

- 

6F15.1) 

11851210 

69  FORMAT ! 14H  CONC.TEMP. 

Ill, 

6F15.1: 

11851211 

74  FORMAT! 14H  EARTH  TEMP 

lit, 

6F15.1) 

11651212 

63  FORMAT! 120H-  EXPOSURES 

1 

2 

1  3  4 

5 

6 

) 

203  FORMAT! 100H0 

EG!  1 )  EG!  2  ) 

EGO) 

1  EG ( 4 ) 

EG!  5 ) 

) 

205  FORMAT! 100H0 

EC!  1) 

EC  ( 2 ) 

ECO ) 

1  EGI4) 

tC  ( 5 } 

EC  16 1 

) 

204  FORMAT! 1 OH 

6F15.3) 

210  FORMAT! SOHO  DBV 

DPV 

wav 

CFM 

)  11851.M50 

75  FORMAT! 10F10.2) 

76  FORMAT! 60H0  08S 

DPS 

MBS 

RHS 

ET 

2  TS 

> 

7100  FORMAT! 10H- 

» 

62  FORMAT! 10F10.0) 

1 1891202 

211  FORMAT ( 3F 1 5. 1 ) 

158  FORMAT! 50H-  OBSERVED 

SHELTER  CONDITION 

) 1 1851M53 

159  FORMAT! 50H-  CALCULATED 

SHELTER  CONDITION 

) 11851M55 

5000  FORMAT! 72H1 

1 

*""  61  Format! ioohi  time 

QVS 

QVL 

QGS 

QGl 

) 

2  QWST  QWLT  TWCT 

QAS  QAL 

) 

6JjftP  FORMAT!  50H1  A01ABAT1C 

SHELTER  RESULTS 

1 

TWO  FORMAT (90H0  TIME 

CFM 

ZN 

DBV 

W8V 

1  OBS  MBS  EFT 

) 

7001  FORMAT ( 10F10, 1 ) 

TT599  FORMAT! iOF12. 21 

READ  5000 
PRINT  5000 


READ  2 1  NN 

READ  2  •  1 1  f  J.I » KK, KMC 
P.EAO  l ,  Ut  YAI  .WCI.WC2 
PRINT  158 
IF<WC)  3,3,4 

3  READ  1,  TGO 

00  5  H-1,5 

00  5  I-l  ,10 

*/  T6(r:*)  .TGo - 

GO  TO  T 

4  00  1000  M-l,5 

1000  READ1.1TG! l,M),J«l,lO> 

?  READ l , I  T60i  M)  ,lt-  l  ,6) 

DQ  8  1-1,5 

- - - “DO  TS*  -M-T.fi  ' - ' - 

8  ?C!L,M)«TCO!M) 
i  1  MiO  l ,  { SI  *1,K«  l ,  6 ) 

READ  l , (GDI K ) , K- l , 6 ) 

*"■  ftCAO  t,  (OGtlO  ,A-1,6) 

READ  1, <0(K),K-1,6> 

- “~R¥*TTT7TZinc  F,  K-T.6 ) . * 

READ  1, < AG ( K ) ,K» 1 ,6) 

READ I, UC(KJ,K«i,6> 


H651H33 


11851013 
U851M35 
11851015 
11851016 
1185101? 
1 1851H38 
1 1851K36 
1J851M37 
1 1851M39 
l 1851N40 

U651M41 

11851043 

11851044 

11851045 

11851046 

11851047 

11851048 

11851049 
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READ  li  ICKIK) ,K«1,6) 

READ  1,  (CGIIO  ,K*1,6> 

READ  li (HIK),K~l,6k 
READ  1, CHDIK},K»1,6» 

P8-29.97 

READ  1,0T,DTT,GS,GL 

C0l»1.08 

C02*4.5 

C03-1050. 

U0*0. 

TN»70. 

2LEMS-I. 

KCC«-l„0 
WCD»~1. 0 
CALL  PREP 
TINE  ID  *6. 

DO  9  N*1,NN 
TV(N)*DBViN) 

9  T IME IN>  »TI HE  I N-l I +DTT 
iFUI)  5005*5005,5006 
5006  READ  1,  GA,ER,CF,OATT,TWt 
5005  CONTINUE 
PRINT  12 
PRINT  13 
PRINT  18 
DO  21  1*1,10 

21  PRINT  20, (TGI I,M),M*1,5» 

PRINT  14 
FR.INT  19 

DO  22  L *  I. ,  5 

.2  PRINT  20,  ITCH., Ml, Ml, 61 
PRINT  15 
PRINT  16 
DO  5015  N=  1 , NN 

5015  PRINT  2  0 , XDB ( N) , T V ( N ) ,OPV(N) 

PRINT124 
PRINT125 
D012  7  K*  1 , 6 

12  7  PRINT126,S(K)  ,0 1 K )  .  7.L  » X )  ,CM  K 1  ,f  G  ( K » ,  AC  I K) ,  AGC  Kl  ,  HU)  ,HOU  ) 
"  26  DU  26  K* 1 , 6 
OX ( K } *D0 ( X ) 

AS  ( K  I *AC ( K  J 
28  CS(;U«CK(KI 
OST-OT 

00  801  K» 1  , 6 

86 1  E  Ul  « AC ( K ) *0 T / 1  PD ( tO  *00 ( K » » 

00  802  K*l,5 

802  EG(K»«AG(K»«DTT/(DG(K(*OG(X) } 

ECMAX-E ( U 
JX*1 

00  58  K«2,6 

IF ( 6 (K> -CCNAXI  58,58.59 
59  ECMAX*F <K) 

JK«K 

58  CONTINUE 

IF (ECMAX-0.5  >  91,91,92 

92  DST-0.4»00( JRI»OD( JR )/AC( JK) 

00  93  K-l  ,6 

93  ECK1-AC (Xl»DST/(OOCK»»OD(K) ) 

91  CONTINUE 


11851050 

11851051 

11851052 

11851053 


l 1851*1 5 


U851M31 


11851066 

11851069 

11851070 

11851071 

11351073 

11851074 

11851075 

11851076 

11851077 

11851078 


11851085 

11851086 

11851087 

11851097 

11851100 

11851101 

11851102 

11E51176 

11851184 

11851*16 

11651M17 

11851*18 

1 1  351*20 
11851*21 
11651*22 
11851*23 
11851*24 

1 18511126 


C-33 


11651105 


PRINT  203 

PRINT  204, (EG(K) »K*lf5) 

PRINT  205 

PRINT  204, (EIK),K*1,6) 

TM-0. 

00  202  K=l,6 

201  TSS (K)*TC ( 2» K  > 

202  TS<K)*TCI t ,KI 
00  800  N*l,NN 
TM»TM+DTT 
Xl=D8V( N) 

X2*WBV(N) 

G=CFM(N) 

ZP*DPV( N) 

PSV*PVI  Z<>,PB>  11851107 

HV=0.622*PSV/(P8-PSV) 


dtm=o. 

53  DTM-DTM»OST 
I  R=  1 
OTA=1«0 
T  Al  =  T  A-20. 

30  QVS= 1.08*G« ITAl-TVIN) 1 


11851108 

11851111 

11C51109 

11351112 


IFI ZX)  34,35*33 

33  TA1-TA1+DTA 
ZX2-ABSFIZX) 

GO  TO  30 

34  TA2*TA1-0TA 
ZX1*ABSF<ZX) 

TA*«TAl*ZX2+TA2*ZXl)/tZXl*ZX2) 

1R*2 

TA1-TA 

GO  TO  30 

35  TA*  TAi 
35  D9PA-0.2 

OPAl*DPA-lO. 

IR*l 

37  PSA*PV(0PAl»P8) 


11851120 
11851121 
11851122 
11651123 
11851124 
11851125 
11851126 
11851127 
11851128 
1 1851129 
11851130 
11851131 

11851133 


WA*0.fc22*PSA/ (P8-PSA) 

QVL*4780.»G>»IiM-*eV) 

CALL  DB0PWB(TA,0PA1,WBA,MT> 

5011  CALL ^  A’ RC0nIga*6» TA^ 1 WBA* Xl*X2»TMl»ER»CF*QASf QAL  »0TAA1 
GO  TO  5012 

5010  QhS=0.  . 

QAL=0. 

5012  CONTINUE 
QWLT*0. 

TwC7*0* 

00  42  K»',6 

PSW*=PV<  T5IK)  ,PB1 
WS”0.622»PSW/ I  PB-PSW) 

ZP-WA-  WS 
IF(ZP)  39,39,40 
39  CWL(K)=0. 


GO  TO  41 

40  QWL I K 1 = 1081. *H0( K)*S(K) *2P  . 

41  QWLT=QWLI< OML(K) 

TWC t  K)*QWL ( K ) / 1061 • 


42  TWCT=TWCT+TWC IK) 

CALL  XTBLU(TA,5.?60,9,Q5,0GP) 
QGL=*  (400.-QGP+GL)*ZN(N) 


11851135 

11851136 


11851137 

11651138 

11851134 

11851141 
11851142 
11851143 
11851144 
11851145 
11851146 
11851147 
1 1B51148 
11051149 


GO  TO  I  43,47)  ,IR 

43  ZY=QGL-GVl-CWLT-QAL 

IFIZY!  45,  46,44  _ 

44  OP A 1=DP Al+OOP A 

ZY2=A9SFIZY) 

GO  TQ  37 

,5  OP  *  2=0?  Al-DOP A 
ZY1=A3SF(ZY) 

OPA* (pPAl«ZY2+DPA 2*Z Y 1 ) / ( Z  Y 1+2 Y 7 . 

’”1*1*2  .  ' 

DP A  5 -DP  A 
GU  TO  3 i 


11851151 

11851153 

11351154 

11851155 

11651156 

1185115? 

11851158 

11851159 

11851160 

11351161 

11851162 


WB*WBA 

CALL  EFTH0B,*B,Er511 


DO  50  K*1.6  11851 166 

FLX(K1*(QW$(K)+QWL(K) 1/SiKl  11851167 

ES(K)*ASIKl*OST/tOX(M»OX(iO  1 

TS(K>»2.*ESIKl*(FtX(Kl«DX(K)/CS<K»*TS$(K)*TS<tU*(0-5<M./eSnO-U  11 

49  TC(l,tO*7SIKl  1 1851172 

50  CONTINUE  11851173 

IF(KTEST)  6985,6985,6984 

6984  CONTINUE 

PRINT  8999,  <T S(K)  ,K*1,<.  ' 

6985  CONTINUE 

51  00  52  L  *2 « 4 
00  52  K  =  t ,  6 

52  Tt(L,K) »E(K)*(TC(L-l,K)«-TClL*l,K)M 1 .  /  E  t  tO-2 . 1  *TC  t  L  ,X  >  >  11851178 

DO  8002  K*  1 » 6 

8002  T<r ( K 1*  TC ( 2 »K ) 


IF(QTH-OTT)  53,54,54 

11851179 

54 

00  55  K*  l ,  5 

l 1851180 

TG(l,KJ*FN(CXIK),CG(Ki,TC(4,K),TG(2,XI,00(K),DG(X),0.) 

11851181 

55 

TC($,K)»TGU,K) 

11851182 

TGI  1 ,6) *TC (4,6) 

11851*46 

TC(5,6)»TG(1,61 

1  1851*47 

56 

DO  57  1*2,9 

DO  57  K-1,5 

57 

TG( I  ,K1  »EG(K) *(TGI I-l,X)*TG( I  * 1 , K  » ♦ 1 1./EGIKI-2. 1 • TG 1 I.K) ) 

11851186 

60 

PRINT  61 

T5X*0. 

OU  8001  K*  1 ,6 

11851198 

8001 

T5*«TS*+TS (K) 

IS*’TS*/6. 

PRINT  62  ,  TM,5v5,SiVL,0GS,aGL,0WST,CWLT,T¥CI  .OAS.QAl 

PRINT  159 

PRINT  76 

PRINT  75, TA, OPA,WBA ,RHA, S F S l, TS* 

11851*49 

PRINT  6*3 

11851203 

PRINT  64*<«WS<K>  ,K*l,6) 

11851213 

PRINT  65  ,  t  OWL ( K I  ,K»t,6l 

11851214 

PRINT  66,ITWt(K),K*l,61 

1185121' 

PRINT  6  7 ,  (  FL  X  (  K  1  ,  K*  i  ,6 1 

11851216 

PRINT  68  i  i T S I K )  ,K*I ,6! 

PRINT  7100 

1185121 7 

/l 

OU  70  L  *  2 , 5 

11851219 

70 

PRINT  69,  L, (TC(L.K) ,K-l>61 

PRINT  7 i 00 

11251220 

72 

00  73  i  *  1  >  1  0 

1.851221 

C-36 


11851222 


73  PRINT  7A,1»<TGU»K).K-1*5) 

PRINT  75»TVIN) ,QPV(N) »WBVIN) , CFMIN) 

P"0  CONTINUE 
PRINT  6000 
PRINT  7000 
00  8000  N=1»NN 
X1*03V<N) 

X2=NBVCN) 


GH=CFH( N) 

SOOO  ^Tnfo?r;!;Em?^mN1.ZN,N,.XUX2.DBS,-B5.ET5l 

CALL  SYSTEM 


END 


11851223 


U85122* 


o -3  / 


SUBROUT  I NE  PREP 

0  [MENS  I  ON  XWS(l50),XHA<150)fXHS< 150)  ,XHWC1S0<  ,XFS(  150)  ,XVAU50)  ,XV 
1 S ( 1 50) , QS ( 50 ) ,CFMC300),ZN(300>, 

20BV<300)  ,DPV<  300),  OBC(  3,300  ),OPC(  300) ,  PS  1 10,300.) 

3  , LA (300, 10) , LB (300, 10) , XOB ( 300 ) , WBV( 300 ) ,E0B(100) ,EWB( 100 

4  ) ,  E T ( 50 , 50 ) , WBC ( 2 , 300 ) , S ( 6 ) 

COMMON  XWS,XHA,XHS»XHW,XFS,KWS,KHA,KFS.KHW,KHS,EDB,EW8,ET,  KET »Nl 1851KG9 
1ET,QS,C0L,C02 ,UC!,TW,GS,GL,ZLEWS,CFM,OPS, DBS, WBS, ETS 1 , 6TS2 , WCC , C03 , l 1851M1 0 
2wCD,DPA*ZN,QG» OGL ,  OBV, 0PV,08C,0PC, W8C ,  S ,TA, OT,KHC ,NN, XVA , 

3XVS , KVA ,KVS,PStLA,LB,XDB,WBV,WCl,HC2  , DTT 

26  FORMAT ( 120H0  TIME  OBV  WBV  DBS  WBS 

1  ET  OGT  OVT  QWT  QWS  TSM  HTX ! 

50  FORMAT ( 1017)  11851002 


52  FORMAT ( 7F 1 0. 1 ) 

1  FORMAT <  10*-  7. 0  > 

2  FORMAT ( 1017) 

READ  2, KWS ,NWS 

READ  1, ( XWS ( [ ) ,1*1, NWS) 

READ  2 , KHA , NHA 

RE AO  1 , ( XHA ( I ) , 1*1, NHA) 

READ  2,  KHS , NHS 
'READ  1, (XHS(I  )  *  I  =  1 » NH  S ) 

READ  2, KFS , NFS 

READ  1, (XFS( I ) , 1*1, NFS) 

READ  ?,KHW,NHW 

READ  1, ( XHW ( I ) , I  *  1 , NHW  ) 

READ  2 i KVA , NVA 

RE  AO  1,  l  XV  A(l'),  1*1,  NVA) 

READ  2.KVS.NVS 

READ  1, (XVS(I) ,I*l,NVS> 

READ  2,KQS,NQS 

READ  1, (QS ( I ) ,1=1 , NOS) 

jlEAD  1,(EDB(I),I*1,29) 

READ  IV  (EW8U),  1  =  1  ,30) 

DO  73  1=1,29 

73  READ  1 ,  ( ET  < I , J ) , J= 1 , 30 ) 

IF(WC.l)  4,4,5 

4  READ  1, (CFM(N) ,N*1,NN) 

GO  TO  6 

■'  5  READ  1 ,  CFMO 
00  7  N=l(NN 

7  CFM(N)=CFMO 

6  I F ( WC2 )  8,8,9 

8  READ  l, (ZN(N),N=1, NN 5 
GO  TO  10 

9  REAC  1, 7 NO 
DO  11  N  =  l »  NN 

11  ZN ( N ) =Z  NO 

10  CONTINUE 

28  FORMAT ( 12F10.2! 

_ T  I  ME  =  0. 

ST=S( 1) +S(2 )+St 3)+S( 4)+S(5)+S(6) 
iF(KWC)  61,61,  62 

61  DO  300  N=1 , NN 

r-.  50,  (LAIN,  L),L=  1,10) 

READ  50  » ( LB ( M , L ) v  L= 1 , 10 ) 

GO  TO  500 

6.:  RE  AD  1  ,  I  D 8  7  (  H  )  ,  N=  1 , 240  > 

READ  1,(VBV(  M),N=i,240> 

READ  i,  (D?C!  ! .  ,Ni ,N=!,240) 


11851002 

11851024 

11851011 

11851012 

11851015 

11851016 

11351017 

11851018 

11851019 

11851020 

11851021 

11851022 

11851023 

11851024 

11851025 

11851026 

1185102? 

11351028 

11851029 

11851030 


11851048 


11851C1-3 


11851F05 


1185JF09 
1 1 85 IF  1 0 

ll85ir-12 


St  AO  i  ,  t  M»t  (  i  i  N  i  ,  *■  L  ,  240  ) 

llessF  i  3 

r-*0  1  ,  t  DSC  (  2  ,N!  ,  N-  1  »  240  ) 

l 1651F14 

St  AO  1 ,  1 *BC< 2 ,N) ,N*1 ,24^) 

Kf  AO  i , 50SCI 3,NI ,N> 1,240) 

SEAO  L , ! SOB* N) ,N» 1  ;?40 ) 

1 1  4  5 1 F  1 5 

0063  J»1,10 

1 1 85  IF  1 8 

63 

HEAD  l, (P$( J,N),N=1, 240 

00  600  N*1 ,240 

1 1 851 F 1 9 

600 

CFM{Nl=CFM(N)/2. 

500 

PRINT  26 

00  400  M*1 , NN 

IFIKWC5  401,401,402 

1185101* 

401 

DBVIN)*LAIN,i )/10 

X«DBV(N 1 

WB V 1 N ) *  L A I N , 2 ) /  10 

Y- W8 V I N ) 

CALL  DB*80P< X, Y,Z ,WQ) 

OPV I  N )  =  Z 

1.851010 

X=(l.A(N,  3)  nAIN,5)*LAiN,7)+LAIN,9)  )/40 

OBSM=X 

11851012 

Y=(LAIN,4)  ♦LAIN,6)aLA(N,8)«-LAIN, 10)  >/40 

WBSM* Y 

11951014 

CALL  DB WBOP ( X  ,  Y ,Z ,MQ ) 

OPSM-Z 

11851016 

X*(LB(N,mLBl*,2)*L8(N,3i+LB<N,4)*LBIN,5)4LBIN, 

6)*LBIN, 7) 4L BIN, 8)11851018 

1  fLB(N,10) )/90 

11851019 

1 F ( L 8 1 N , 4 ) )  53,53,54 

1 1.851020 

53 

X=X*9./8. 

54 

T  SM=  X 

XDBIN1»L8IN,9)/10 

GO  TO  70 

11851F07 

402 

X*OBVIN) 

Y-WBVIN) 

CALL  08WB0PIX,Y,Z,W0) 

OPVINl-Z 

X* ( OBC  (  1 , N ) ♦ DBC ( 2 , N ) *DBC I  3 ,  N  >  )  /  3. 

OOSN»X 

Y-IMBCI 1  ,N)*KBCI2,N) 1/2. 

WBSH-Y 

1 1051F23 

CALI.  DBW0DPIX,V,Z,NC) 

OPSN-Z 

1 1 851F2  7 

SUN-O. 

11851F31 

SUN-O. 

HB51F32 

DO  64  J-1,10 

1 1 85  IF.’  3 

I  F ( PS ( J ,N ) )  64,64,66 

1 185  IF  34 

66 

SUN*SUNaP$< J,N) 

1 1851F35 

SUN  > SUN* l . 

11851F36 

64 

CONTINUE 

1 1851F37 

67 

I SN* SUN /SUN 

*1-* 

'  n*4  T  !  NU  c 

UPV3fl»«UP5N 

1 1  Nr, '  Fox 

TlN£«'INE40Tt 

X»DBVIN> 

Y40PVIN) 

Z-NBV(N) 

CALL  TBLUI X , X  MS , XKS , WSS ) 

11851106 

CALL  T6LUIY,KmS,XNS,NSV) 

11851107 

CALL  TBLU(X,KVA,XYA,VAi; 

11851108 

CALL  T8LUI X.KVS.XVS.VSi) 

11851109 

V«VA1*I VSI-VA1)«*SV/NSS 

11851110 

C-J9 


CM  L  Ff?l(XX,ZZ,Sfn 
C  AU.  T  B !.  U  (  y  f  ,  K  W  S  ,  X  w  S  ,  W  S  A  i 
CAM.  XffL,MxXfS.,60,9tQ,i,0GSi 
OGT* i 400, *GS*GL ! »?N( N! 

UGS  *GGS»?N(N) 

wr,L  =  (400. -0GS)*ZMN)/  1053.? 
l)VS  =  Cf:M(N)/V»0.24«60.»(XX"X) 

WVL=CFM<Nl /V»60. »  <  WSA-WSV ) 

UVT=CVS*WVL»1054. 

UWS=QGS -0VS-GS 
OWT  =  QGT -QVT 

nrx=uwr/sr 

RAT  I O-QWT/OGT 

!Hrx,'IT  •?8,T,ME'0BV,N,'k'BVi;'',»D&SH»wes«,SFT,0GT,(}Vr,Qwr.Q«SfTSMf 

00  CONTINUE 

ta=dbc(  i,n 
DPA=DPC (  I  ) 

RETURN 

END 


U85lu2* 

12851113 

1165111* 

11851118 

1185111? 

naaius 
1 1851116 
11851026 

11851027 

21851028 


11851029 


c  SHELTER  AIR  CONDITION  CALCULATION 
SUB  ROUT  INC  STRCCN!OB,WR,2fM) 

0 1  MENS !  CN  XW$! 150)  ,XHA( 150) ,XHS< 150) ,XHW( 150 )  ,  XF S < 1 50 ) , 

IEDB! 100) ,£W8! 100) ,£T(50,50) ,0S(50)  ,CFM<  300) 

COMMON  >  WS . XHA  t  XHS  «  XMW  ??Ir  j  ,KWS,  K  HA  tiers » KHW,  KHS,  CDB ,  E  WB  ,  ET  , 

3KET  .NET  » US  ,COl , C02  .UG , TW , GS  ,GL,  ZLEWS .  Cf  M,  OP S,  DBS , KBS, E TS 1 ,  ET$ 2 


4wCC , CG3 , WC  0 
12  FORMAT (  10H 

ETSl 

F10.2) 

13 

FORMAT  1 10H 

ETS2 

F10.2) 

9 

FORMAT! 10H 

QGL 

F10.2) 

4 

FORMAT! 10H 

DBS 

F 10 . 2  1 

5 

FORMAT! 10H 

9GS 

F 10. 2 i 

14 

FORMAT!  *  OH 

DPS 

F10.2) 

6 

FORMAT! 10H 

THS 

F 10.2  ) 

7 

FORMAT! 10H 

DP 

F10.2) 

11 

FORMAT! 10H 

WBS 

F10.5) 

10 

FORMAT! 10H 

WA 

F10.5) 

8 

FORMAT! 10H 

WV 

F10.5) 

YI.  =  C01»ZFM+UC 
Y2=C01*ZFM*0B+UC»TW 
Y3= Y2/Y 1 

IF! Y3-100.  )  30,30,31 

30  X1=Y3 

3  X2  =  X 1  +  0 . 5 

CALL  XTBLU(Xi»5. ,60»9,GS,QS15 
CALL  XTBLUtX2 ,5, ,60,9,95,052) 
Z1=!QS1+GS)/Y1+Y3-XI 
Z2=!QS2+GS)/Y1+Y3-X2 
IF ( Z2 )  1,1,2 
7  X 1- X2 
GO  TO  3 
l  22- aBSF I Z2  ) 

CBS»Xi+0.5*Zi/Z2/(l.+Zl/Z2» 

CALL  XTBLU(0BS,5, ,60,9 , OS , 9GS I 
GO  TO  32 

31  DBS*U400.+GS+Y2)/IY1  +  14J  ) 

QGS*-14.*!DBS-100.> 

32  CONTINUE 

CALL  TBLUI TW , KWS , XWS « T WS ) 

CALL  DBWBOPiOB.WB.DP.WV) 

0GL-400.-QGS 

ZF«ZlEWS*U0/0.243 

WA* ( ( QGL+GL ) /C03  +  ZF  »TWS+C02+ZFM*WV )/ (C02»£FM+ZF  ) 
IF ( WCC )  19,20,20 

20  PRINT  10, WA 
19  CONTINUE 

1FIWA-TWS)  101,101,102 

101  .-CD)  ,03, 103,  lvic 

103  WA*WV+!QGL+GL)/!C02*ZFH*C03> 

102  CALL  DBWWBH(0BS,WA,W8S,H) 

IF (WCC)  21,22,22 

22  PRINT  1 1 , WBS 

21  CONTINUE 

CALL  OBWBDP(OBS,WBS,OPS,WA) 

IF ( WCC )  23,24,24 

PRINT  14, OPS 
PRINT  10, WA 

23  CONTINUE 

IF ( WCC !  25,16,16 

16  PRINT  4.0BS 


! I  85 1 fc 02 


1 1851E05 
,11851028 


1 1851 

11851E10 
1 1 85  IE  1 1 

1 1851E14 


1 1 85  IE  1 5 

liesieib 

11851E1T 

11951 

11851 

11851 

11551 

1 1 85  IE  l  3 

11851E19 

U851E20 


1 185IE24 


C-41 


PRINT  *>  ,vJOS 
m  I  N  I  ft  ,  T  R  s 
PR  .'.NT  7,  DP 
PRINT  8 ,WV 
PRINT  <J,  COL 
PRINT  JO, WA 
PRINT  1 1 , WBS 
PRINT  14, DPS 
25  CONTINUE 

CALL  EP T1 1  DBS, MBS  ,ET<  1  ) 
IF(WCC)  15,17,1/ 

17  PRINT  12.ETS1 
15  RETURN 
END 


r  Oh 

SUE  ROL  ■  T  ;  Nr  A !  ROL-t  i  GA  »k W  »GP  *  C  F  ’G*»DBG»WBS*TvI  *QAStQAL» NS  OH  ) 

GA  =  AIH  PLOW  THROUGH  CO'L  *  C  H  M 

Gw  =  COOL  AMT  FLOW  i:  «S/KS,.  (artUHC  HEAT)  >  6  TU/HR /DEG- 500*GPM  FOP! 
RW  =  CC‘!L  i Ht RMAL  RESISTANCE  EXCLUDING  AIR  FILM  RESISTANCE 
SP  =  CO  1 L  A 1 17  SIDE  HEAT  TRANSFER  SURFACE  PER  ROW 
CF  =  COI L  CONTACT  FACTOR  PER  ROW 
NRwW-  <  0 1  i.  PD  WG 

DB  S » W  R  S  •  -  C  0  I L  INLET  AIR  DPY-AND  WET-BULB  TEMPERATURES 
T  W  I  =  C  0 1  L  COOLANT  ENTER  I NG  TEMPERATURE 

C AS  -COIL  SENSIRLF  CAPAriTY 


GAL  -COIL  tATFNT 


CAPACITY 


3  I  AS  100.1 
11.851002 

i  -I  s  L  5  f;.(L  > 

H2Q1 1 85 1 004 
U8310G5 
11851006 
3165100? 
11853 008 
11851009 
11851010 
11851011 
1 1851012 


DIMENSION  XW5 f 1 50 1  *  X H A (150! >XHS( 1 50 ) 

» XHW (150) *  XFS ( 1?0! 

.T (  20  3  *W<  2011851013 

1  )  t  T  W  1  2  0  )  .  T  AL  f  2  )  *  TWG  !  2  )  ♦  FGS  1  2  ) 

11851014 

COMMON  XWS>XHA*XHo*XHW,XFS.KWS*KHA.KFS.K;HW.KHS 

11851015 

T ( 1 ! =DBS 

11851016 

CALL  DBWBDPfDRStWRS.DPS.WX ) 

11851017 

W!  1 ) =WX 

11851018 

TWG ( ]  ) =0ES-1  , 

11851019 

TWG ( 2 ) =T WG ! 1  )-l  . 

11851020 

A  =  S  P  /  (  R  W  *  G  A  *  C  F  *  4  .6  ) 

11851021 

7 

DO  ICO  K  =  1  » 2 

11851022 

T W ( 1 ) =  TWG ( X ) 

1185102? 

DO  90  I  -  1 » NROW 

31851024 

X  =  T  (  1  t 

11851025 

Y=Wi 1 ) 

11851026 

2  =  T  W  (  I  ) 

11651027 

CALL  5 UR FT (A»X  *Y ,7  »  T  S  X » WSX ) 

11851028 

T  !  1  +  1  ) -T ( I )-( T ( I 1-TSX! *CF 

11851029 

l =WSX-W ( 1 ) 

11851030 

I  F  !  Z  )  1  .  ?  .  2 

11851031 

1 

W  (  1  + 1  ) =  W ( I ) - { W ( I 1 -WSX ) *CF 

11851032 

GO  TO  A 

11851033 

7 

W(  It  1  1 =  K( I  i 

11851034 

-J 

TW» r » - T  W ( !)-  '  5”/ !GW*RW)  :*( 

T  CX  •  *»«!), 

118*5108?) 

9  <3 

CONTI NUF 

11851036 

NZ--NRCW.  1 

11851037 

FGSCK  )  =  r W ( N Z  J - T W I 

11851038 

100 

FALU >  =  T ( NZ  » 

11851039 

TEST-FGS! 1 )*FGS(?' 

11851040 

IFITESTi  4.5.6 

11851041 

6 

TWG ( l ! sTWGI ?) 

1105104? 

TWG(?  )«TWGI?)~1. 

11851043 

GO  TO 

11851044 

9 

TWL*TWG  t 

11851045 

T  L  a  T  A  L  <  ?  1 

11851046 

GO  TO  96 

11851047 

A 

C*ABGF(FGG( 1 i/FGSI?)  ) 

11851048 

T  WL*  f  Twr,(  1  •  t  li*G(  2)*C  )/(  1  .4-C) 

11851049 

TL-(Tal(1 )  +  Tal (2! > / 2  * 

11851050 

98 

O,aS--1.OF*GA*(D0S-TL) 

11851051 

QAL=C-W*  (  TWL-Tw  1  l-OAS 

11851052 

return 

11851053 

END 

11851054 

C-43 


c 


SwiROUT  f  Mr-  Stit«FT  (  A  »TA.WA  .Tw.TS.WS) 

SOLVE  FUN  TS  FROM  0 . 243#  i  TA-TS )  +  1 060 . # < WA-WS  )  =  A* ( TS-TW ) 
DIMENSION  XWS  I  150)  150)  *XHS(  150)  . XHW  <  1 50')  . XFS  (  1 50  )  .Yt  2)  «F<  2  ) 

common  xw$*xha»xh;>>xhw»xfs»kws»kha»kfs  .KHW.KHS 

V  (  l  )  *  T  A 

Y  !  ?  )  <j  TA-1 . 

5  DO  1  t  *  1  > .? 

Pl«PV( Y( ) ) .29,02) 

WO.  62?«<*»/l29.  92-PI  i 

1  F  (  l  )  ■  A*M  Y  (  !  1“TW)-0.24?*(TA-Y{  I  )  )  ■*■!  C60  •  *  (  WA-W  ! 
r FST*F( 1 !»"(?) 

IF (TEST)  2  *  5 • A 
4  Y  ( 1  )  =  Y  (  ?  ) 

Y I  2 ) s Y ( ? ) -I » 

GO  TO  5 
3  TS  =  Y( 2  ) 

GO  TO  10 

2  C*ABSF(F( 1)/F( 2) ) 

TS=(Y(1)+Y<2)*C)/(1.+C) 

10  PI =PV ( TS. 29.92 ) 

WS*0.622»P1/ I 29.92-PI ) 

RETURN 

END 


11851001 

11851002 

1185100? 

11851004 

11851005 

11851006 

11851007 

11851008 

11851009 

11851010 

1185101) 

11851012 

11851013 

11851014 

11851015 

11851016 

11851017 

11851018 

11851019 

11 851020 

11851021 

1  1851022 

1185102? 

11851024 


C  -44 


EFFECT ! V £  TEMPERATURE  CALCULATION 
SUBROUTINE  EFTi(CB,WB,EFX) 

D1MENSIFN  XWS(  1 50 ) , XHA ( 1 50 ) , XHS ( 150) ,XHW( 150 ) , XFS ( 150 ) , 

FQB( 100) ,FWB< IOC) ,ET<50,50> 

-OMMON  XHS ,XHA, XH S , XHW , XF S , KHS , KHA , KF S , KHW , KHS , EDB , EWB , ET , KET , NE T 
IFiDB-WB)  1,1,2 

L  EFX=DB 
GO  TO  15 

2  I F ( DB-ECB (  29  )  )  3, A, A 

3  EFX=0. 

GO  TO  15 

4  I F ( WB-EWB (  30  )  )  3,5,5 

5  I F ( DB-ECB ( 1 >  )  6,6,7 

1  EFX  =10C. 

GO  TO  15 

6  i  X  —  0 

UO  8  1=1,29 

I F ( DB-ECB ( I ) )  9  ,10,10 

10  X 1=ED8 (  I ) 

X2=E0B ( l-l ) 

GO  TO  11 


9  IX=IX+1 
8  CONTINUE 
11  JY=0 

UO  12  J  =  1 , 30 
IF ( HB-EHBl J) )  13,14,14 

14  Y 1  =  E  WB ( J ) 

Y2  =  FW6<  J-l  ) 

GO  TQ  K 


^Y  = JY+1 
1  ^CONTINUE 
16  211 =ET(  1X4 1 , JY+1  ) 


2  1 2  =  E  T  IX,  JY+1) 

Z22  =  ET (  IX, JY) 

2  2 l  =  E  T (  IX+i, JY) 
IF(ZU»Z12«222*Z21)  3,3,17 
17  22=222+ i DB-X2 ) * ( Z21-Z22 ) / ( XI- X2 ) 
2 1  =  Z 12+ ( DB-X2 ) * ( Z11-Z12)/ (X1-X2I 
t  F  X  =Z2+(WB-Y2)»(Z1-Z2)/(Y1-Y2) 
IS  RETURN 


NO 


11651  44 


11651  45 


I 


NfHALPY  CALCULATION  BY  CB  AND  WB 

1 1051GO2 

SUBROUTINE 

DBWBH(DB,WB,H) 

11B51G03 

01  MENS I  ON  XWS< 150) tXHA( 1501 ,XHS< ISO > , XHW{ 1 50 ) , XFS ( l 50 ) 

1 1851G0A 

COMMON 

1  xws 

»XHATXHS»XHWfXFS»KWS» KHA»K  FS*  KHM»KHS 

U351G05 

;all 

T3LU 

(WBtKWS.XWS.WS) 

1 1 851G06 

CALL 

TBLU 

(WB,KHWfXHW,HW) 

11351G07 

CALL 

T6LU 

(DB»KWS»XWS»WSS) 

1 1851G0B 

CALL 

TBLU 

( OB»KHA  »XHA»  HA ) 

11851G09 

CALL 

TBLU 

( WBf KHS  «  XHS»HS  ) 

11851G10 

CALL 

TBLU 

IDB,KHS,XHS»HSS) 

11851G1  1 

HAS=HSS-HA 

1 185 1G  1  2 

X=(HS- 

Hk*WS)*HAS-HA*HW*WSS 

1185101  1 

Y=HAS- 

Hh*HSS 

118  5 1 G  1 A 

H  =  X/Y 

11831G1  5 

RETURN 

h-  _ 

11851G16 

END 


C-47 


RELATIVE  M'iHlniT>  C ALCUl. AT  i ON  8r  06  ANO  WR  1)851811 

SUBROUTINE  GfihRRH«03,fcB,P8,8Hl  11851512 

DIMENSION  XWS  (  150)  ,KMAU50>,  XHS  i  i  50  )  ,  XHW  {  t  50 )  ,  XF  SI  1  50  I 

COMMON  XWS,XHA,XES,XHW,XFS,KH5,KhA,KFS,KHW,KHS  1I851SIA 

„6LL  08 1* BOP  (  DB  »  R  5  » UD  r  W  ) 

CALL  TPLLML'B,KhS.XWS,MS}  11851BI7 

ZM=W/KS  1 185  IB  L  6 

DB0=D3/10. 

CALL  T8UI(C8BtXFS,XK,,FSi 

Z  =  Py ( D8 . P8  )  /P 8  1  1851320 

RH*100. *2H/! l.-( i.-’Ml»FS*n  11851821 

RETURN  11351B22 

END 


C-48 


H95LF02 


WET  8Ut8  CALCULATION  8Y  OS  AND  W 
SUBROUTINE  06hK8N<08.W,WB,Ht 

DIMENSION  X*$ ( 150 )  » XHA< 150 ! » XHS  ( 1 50> * XWW { 150 > , XF$ (1 50 1 
COHHON  XWS,XHA,.XHS»XHIi,XF$tXW$,KHA,KF$,KH#,RHS  U851F05 

TRLU{DB,KHA,XHA,HA)  11**’*** 

CALt  T8UH0B»KHS,KHStFJS}  11851F0? 

CALI  TBtU(CBf*WS,XHS»WS> 

!F<N-#S)  R,5,5 
5  *8=06 
H=HS 
GO  TO  6 


A  H*NA*«* (HS-HA »/WS 


XI  *08 

11851F10 

X2*Xl- 

0.5 

nesiFU 

C&U 

T8L0iXl.KHS.XHS,HSli 

11B51F12 

CALL 

TBLUl X2  »KHS • XHS  *HS2) 

U351f  13 

CALL 

TBLlHXlfK*5«XNSf*Sl> 

11051F14 

CALL 

TBLUU2,KWS,X*S.MS2) 

11851F15 

CALL 

T8LU(Xl»XH*f  XHNfHW 1 ! 

1 1851F16 

CALL 

r8LU<X2.KHSt,KHW.HN2) 

11851F1T 

l l*HSl 

11851F18 

22=HS2 

-h*2»(NS2-KJ 

U851F19 

t  F 1 Z2-H )  l j 1 f  2 

Xl*Xl- 

0.5 

1 1851F2  1 

GO  TO 

3 

11851F22 

l  *8= X 2+0. 5* (H-Z2)/(?l-22> 
6  RETURN 
END 


C-49 


VAPOR  PRESSURE  CALCULATION  t 

^UNCTION  PV(X,PB> 
r=X+A59,680 
TS=671.68S 

fMl*-<7.<J02<?8J«CTS/T-l»  > 

TM2*( 5. 02808 ) • ( L0G10F (TS/T  J ) 

TM3=—  ( 1.3816/(10.** <7. )))•( UO.«*t 11. 344*1 l.-T/TS >)>-!. ) 
TM4*  <  8.  1328/1000. J*<10. **<-3. 49149*<TS/T-1. ))-!.» 
AN51*TM1*TM2+TM3*TH4 
AN$2*10. ** I  ANSI ) 

PV-P8*ANS2 

RETURN 

ENO 


11051H1' 
U851NH 
11851HH 
1 1851H1  *■ 

l  kb* itiii 

I1851H2C 
U851H2I 
1 16ilH2< 
118SIH22 
11851H2* 
11851H2( 


C-50 


O 


DEW-POINT  CALCUL AT  I  CM  BY  DB  AND  W8  llflSIOO; 

SUBROUTINE  r>nwonp(D8»HB»DP»W)  1185100.' 

DIMENSICN  awS 1 150 ) »XHA ( 150) t  XHSC 150 ) ,XHWI 150 ) ,XFSf 1501  1185100' 

"  OMNON  >WS , XHA , XHS,XHW,XFS,KMS|KHA,KFS,KHW,KHS  1185100! 

iFIOB-WB)  1,1,2  11851004 

1  DP-08  1165100'. 

CALL  T8LU(0B,KNS,XHS,W)  11851004 

GO  TO  100  1185100' 

2  CALL  T8LU(0B»KHA»XHA*HA)  11851014 


1  IS  FULL.  GIVE  ME  ANOTHER  A3  ANO  START. 

PERMANENT  TAPE  REDUNDANCY  HAS  DETECTED  WHILE  READING  THE  RECORD  THAT  PRODUCED 


CALL  TSl.U(OB,KHS,XHS,HS  1 

_CALL  TBLUIOBfKWS,  XWSj  WS)  _ 

CALL  T  8  LU I WB , KH  S , XHS , HSWB I 
CALL  TBLU<WB,XHW,XHW,HW) 

CALL  T6LUIWB,KWS ,XWS,WSW8) 

W-WS»(HSWB-HA-HW*WSMB)/(HS-HA-HW*WS> 

Yl-WB 

Y2-Y1-0.5 


CALL  TBLUI Vl,KKS«XWS,Zl) 
_CALL  TBLU(Y2,KWS_,XWl,i21 
IF (Z2-W) 3,3,4 
4  Yt-Yl-0^5  _  _____ 

GO  fi)  5  " 

3  0P-Y2+0.5* (W-Z2)/ IZ1-Z2 1 
100  RETURN 
END 


1185101) 

1185101; 


1185101: 

1185101' 

1185101! 

11851014 

1185101 

11851011 


ile5101' 

11851024 

1185102 

1185102: 

1185102 

1185102< 

1185102' 


H851S0 
i 18510C 


WET  BULB  C ALCULAT ION  BY  OB  AND  OP 
SUBROUTINE  DBDPWB ( DB,DP , WB,W ) 

DIMENSICN  XWS ( ISO ) * XHAI 1501  *  XHS( 150 J ,XHW( 150} .XFSI 1501 
COMMON  XWSfXHAfXHSf XHW»XFSfKWS«KHAeKES»XHW»KHS 
CALL  TBLUIDP.KNS.XWS.WJ 
CALL  OBWWBHIOBfMfWBfH) 

RETURN 

END 


11851B0 

11851B0 

1 1851B0 


►  ■ 


LINIER  INTERPOLATION 

U851HO 

SUBROUTINE  TBLUIX.KX.VtYS 

DINENSICN  VI  ISO) 

ItSSIHO 

J*XINTF (X) 

I1S51H0 

*.»J -KX*2 

imiHo 

LL-J-KX+1 

U8SIHO 

VU*V{L) 

I 185 LUO 

VL=V(LU 

11851H0 

Y*<VU~VL}«lX-INTF<X)>+VL 

I1B5LHL 

RETURN 

END 

U851HI 

SUBROUTINE  XTBLO < X,DX ,KX , IX, V, V) 

omensiCN  visoi 

J*XINTF  if  X) 

N“KX+X!NTFIDX>»<LX-il 
IF(J-KX)  1,2,2 

1  V»VU) 

GO  TO  1C 

2  IF( J-H  )  3,4,4 

4  Y-V(LX) 

GO  TO  1C 

3  00  6  1*1, LX 
2*1 

P*KX 

0*P+Z*DX 

.  _IF(X-Q»  5,5,6 

5  Ql*V(I+n 
Q2*V(I) 

V'*QI-(Q1-G21»(Q-X}/DX 
GO  TO  1C 

6  CONTINUE 
10  RETURN 

END 


11651C0 

1 1851C0 
11851C0 

1185IC0 
11851CO 
11851C0 
1 1 85 1C  1 
11851C  1 
11851C1 
11851C1 
1 1851C1 
1 1851C1 
11851C1 
1 1851C1 
1185 1C  1 
11851C1 

11851C2 

U851C2 


S«Bpl«  Input  an A  Oitput  Data  on  M-4  Pro*r*n 


UNDERGROUND  FALLOUT  SHELTER  DATA 


INITIAL  eART  H 
TGI l i  T*<3) 
70.5  72.0 

70.0  72.0 

69.0  71.0 

67.0  70.0 

67.0  70.0 

67.0  70.0 

67.0  70.0 

67.0  69.0 

67.0  68.0 

67.0  67.0 


TEMPERATURES 


TGI  3 » 

TGU) 

72. C 

72.0 

72.0 

72.0 

71.0 

71.0 

70.0 

70.0 

7C.0 

70.0 

70. 0 

70.0 

70. C 

70.0 

69. C 

69.0 

68.0 

68.0 

67.0 

67.0 

INITIAL  ''ONCRETE 
TCin  TH2) 

73.0  73. u 

73.0  73.0 

73.0  73.0 

73.0  73.0 

73.0  73.0 


WALL  TEMPERATURE 


TCI  3  > 

TCU) 

73.0 

73.0 

73.0 

73.0 

73.0 

73.0 

73.0 

73.0 

73. C 

73.0 

TGI5) 

73.0 

72.0 

72.0 

71.0 

7i.0 

70.0 

69.0 

68.0 

67.0 

67.0 


TCI5) 

73.0 

73.0 

73.0 

73.0 

73.0 


0-1 


TGI6) 


TC 16) 
73.0 
73.0 
73.0 
73.0 
73.0 


fine  VARIABLES 


KN> 

T'M  N5 

70.0 

76.0 

72.0 

72.0 

74.0 

74.0 

75.0 

66.  0 

77.0 

68.0 

77.0 

66.  0 

76.0 

64.0 

76.0 

63.0 

75.0 

62.0 

75,0 

81.0 

75.0 

91.0 

75.0 

92.  G 

77.0 

93.  0 

78.0 

92.0 

79.0 

88.0 

80.0 

85.  0 

80.0 

81.0 

90.0 

79.0 

80. 0 

77.0 

80.0 

76.  0 

60.0 

79.  0 

81 .0 

84.0 

fil.O 

87.  0 

60.0 

90.  0 

82.0 

92.0 

82.0 

92,0 

83.0 

88.0 

83.0 

8^.0 

83.0 

82.0 

83.0 

80.0 

83.0 

79.0 

82.0 

78.0 

62.0 

80.  0 

82.0 

86.0 

82.0 

91.0 

62.0 

92.0 

6  3.0 

92.  0 

83.0 

97.0 

64.0 

88.0 

84.0 

83.  0 

64.0 

81.0 

84.0 

79,0 

84.0 

80.0 

84.0 

79.0 

84.0 

81.0 

83.0 

86.0 

83.0 

90.0 

63.0 

92.0 

94*0 

92.0 

84.0 

92,  0 

85.0 

86.0 

85.0 

63.3 

65.9 

80.0 

85.0 

80.0 

84.0 

79.0 

64.0 

79.0 

64.0 

79.0 

CPVINl  USUN(N) 

57.2 

47.9 

61.9 

46.6 
55. C 

49. 6 

48.2 
47. C 

43.5 

6  0  c  2 

71.8 

71.4 
71.  C 

68.2 

69.9 

71.1 

71.3 

69.1 

68.4 

71.6 
7C.6 
73. C 

71.9 

72.2 

71.4 
73. C 
73. C 

71.6 

72.3 
7C.2 

69.1 

69.5 
7C.2 

69.2 

71.8 

71.4 
69. e 

71.4 

71.4 

70.4 

69.8 

69.1 

71.7 
7C.6 

71.3 

72.2 
7C.6 

69.8 

71.4 

71.4 

71.4 
7C.4 

70.2 

70. 2 

72.1 

72.1 

72.1 

0-2 


83.0 
63.0 
83.0 
84.0 
84.0 
8*5.0 
85.0 
85.0 
05.0 
85.0 
85.0 
U5.0 
85. C 
84.0 
84.0 
85.0 
05.0 
05.0 
85.0 
85.0 
85.0 
85.0 
85.0 
85.0 
84.0 
84.0 
03.0 
84.0 
83.0 
85.0 
86.0 
86.0 
86.0 
85.0 
85.0 
85.0 
85.0 
85.0 
65.0 
85.0 
85.0 
85.0 
85.0 
85.0 
65,0 
•5.0 
84. 0 
84,0 
04.0 
84.0 
83.0 
84.0 
85.0 
85.3 
86.0 
85.0 

e?,o 

83.0 

85.0 

84,0 


87.0 
90.0 
92.0 
92.0 
92.0 
88.0 
83.0 
82.0 
80.0 
79.0 
79.0 
80.0 
86.0 
89.  0 
92.0 
93.0 
91.0 
87.0 
83.0 
81.0 
80.0 
79.0 
78.0 
80.0 
85.0 
89.0 
91.0 
93.  0 
92.0 
88.0 
83.0 
81.0 
79.0 
79.0 
79.0 
80.0 
82.0 
90.0 
92.0 
92.0 
91.0 
07.0 
83.0 
80.0 
79.0 
78.0 
77.0 
79.0 
83.0 
89.0 
91.0 
92,0 
91.0 
•8.0 
83.0 
80.0 

79.5 

78.0 

79.0 


TU? 

72.2 
69.8 

71.4 

71.4 
73.0 
72. C 
7C.8 

70.2 

70.6 

70.6 

73.1 

72.2 

72.6 

71.4 
71. C 

70.2 

70.3 

70.4 
69.8 

70.2 
7C.6 
71.0 

71.7 

72.6 

72.6 

71.8 
71.  C 

71.4 

71.4 

70.4 

71.3 

70.6 

70.6 
7C.6 

71.7 

73.8 

72.2 

71.4 

71.4 

71.8 

70.3 
72.0 

71.7 
72.1 
71.0 

71.4 
72.1 
72.6 
72.6 

73.4 

71.4 
71.0 

69.9 

72.5 
71.? 

75.6 
71.0 
71. C 

73.5 


84*0  85.0  72.6 

84.0  89.0  72.6 

84.0  92.0  73.  C 

84.0  92.0  71.4 

88.0  92.0  73.0 

85.0  88.0  73. C 

85.0  84.0  7C.C 

85.0  81.0  71.3 

85.0  30.0  70.2 

85.0  78.0  71. C 

84.0  78„0  71. C 

85.0  80.0  71.7 

84.0  85.0  72.6 

84.0  90. (  72.2 

84.0  92. J  71.4 

84.0  92.0  71.4 

85.0  91.0  71.8 

85.0  88.0  69.9 

86.0  03.0  72. C 

66.0  81.0  ■»».3 

85.0  79.0  70.6 

86.0  79.0  7C.6 

85.0  78.0  7i.C 

85.0  79.0  72.1 

86.0  85.0  72.6 

85.0  89.0  72.6 

85.0  90.0  72.2 

85.0  92.0  73.0 

86.0  91.0  71.8 

86.0  88.0  69.9 

86.0  84.0  7C.C 

86.0  81.0  71.3 

86.0  78.0  71. C 

86.0  78.0  72.5 

86.0  78.0  71. C 

85.0  80.0  71.7 

85.0  85.0  72.6 

85.0  88.0  73. C 

84,0  92.0  71.4 

85. C  93.0  72.6 

85.9  92.0  71.4 

86.0  88.0  69.9 

*6.0  84.0  71.6 

86.0  81.0  71.3 

86.0  79.0  69.1 

86.0  79.0  72.1 

85.0  79.0  72.1 

65.0  *0.0  71,7 

85.0  85.0  71.4 

85.0  90.0  72.2 

85.0  92.0  71.4 

•5.0  92,0  91.4 

86.0  92.0  71*4 

•6.0  #9.0  71. C 

*6.0  #8,0  72. C 

«6.0  80.0  71.7 

86.0  “  TI.B  71.0 

•6.0  78.0  71.0 

88.0  T7.0  71.4 

85.0  #0.0  7J.1 


0-4 


85.0 
85.0 
85. 0 
86.0 
86.0 
86.0 
85.0 
85.0 
84.0 
84.0 
84.0 
83.0 
83.0 
82.0 
83.0 
83.0 
b3.G 
03.0 
03. 0 
83.0 
83.0 
83.0 
8  3.0 
83.0 
82.0 
82.0 


36.  0 
90.  0 

92.  0 

93.  0 
93.  0 
88.  0 
84.  0 
81.0 
79.0 

79.  0 

79.  0 

80.  0 
85.0 
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